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I.  Introduction 

The  research  grant  F33615-84-K-1557  entitled  "Microwave 
Amplifiers"  was  awarded  to  the  University  of  Illinois  at  Urbana-Champaign 
by  the  Department  of  die  Air  Force,  VVright-Patterson  AFB,  Ohio.  Dr. 
James  B  .  Tsui  of  Wright-Patterson  AFB  is  the  Scientific  Program  Officer. 
Mr.  John  W.  Michalski,  Resident  Representative  of  the  ONR  Chicago  Office, 
has  been  assigned  the  responsibility  for  the  administration  of  the  contract. 
The  total  research  project  cost  is  $135,000  to  cover  the  period  from 
September  7,  1984  to  September  30, 1987 

1.  Period:  September  7,  1984  to  September  30, 1987 

2.  Reporting  Date:  September  30, 1987 

3.  Technical  Personnel: 

S.  L.  Chuang  Assistant  Professor  of  Electrical  and  Computer 
Engineering 

S.  W.  Lee  Professor  of  Electrical  and  Computer  Engineering 

C.  T.  Sah  Professor  of  Electrical  and  Computer  Engineering 

(consultant) 

Two  research  assistants 


II.  Technical  Discussions: 


(a)  Modeling  of  the  microwave  GaAs  MESFET 

In  recent  years,  the  gallium  arsenide  metal-semiconductor  field  effect 
transistor  (GaAs  MESFET),  Fig.  1 ,  has  emerged  as  a  very  useful  and 
practical  device  in  a  microwave  low  noise  and  power  amplifier,  oscillators, 
and  several  other  types  of  circuits.  The  main  advantages  of  MESFET  over 
previously  used  solid  state  devices  are  its  high  direct  current  (dc)  to 
microwave  power  conversion  efficiency  and  high  power  gain.  The  GaAs 
MESFET's  electrical  characteristics  and  physical  structure  allow  the  device 
to  be  conveniently  integrated  into  hybrid  microstrip  circuits  as  well  as 
microwave  monolithic  integrated  circuits  (MMICs).  In  order  to  facilitate 
the  optimum  circuit  design,  models  of  the  electrical  behavior  at  microwave 
frequencies  performance  based  on  theoretical  physical  principles  have  been 
developed  here.  First,  an  analytical  small  signal  equivalent  circuit  model  is 
obtained.  Then  the  calculated  results  are  compared  to  that  of  the 
experiments.  A  more  accurate  two-dimensional  model  which  is  based  on  an 
approximate  solution  of  the  semiconductor  equations  applicable  to  large 
signal  circuit  analysis,  has  been  investigated.  1-V  characteristic  curves  of  a 
MESFET  have  been  calculated  (Fig.  2)  and  compared  with  the  experimental 
result  shown  in  Fig.  3  [1]  .  The  agreement  is  fairly  good  except  for  larger 
drain  conductance  at  low  drain  voltage,  which  results  in  a  higher  drain-to- 
source  saturation  voltage.  This  is  due  to  the  parasitic  source  and  drain 
resistance  which  are  not  taken  into  account  in  the  current  voltage  relationship 
calculated  for  the  intrinsic  FET.  The  scattering  parameters  were  calculated 
for  the  complete  GaAs  MESFET  in  Fig.  4,  which  consists  of  an  intrinsic  FET 
with  the  parasitic  circuit  elements.  A  comparison  of  the  S-parameters  with 
measurement  of  a  similar  device  are  plotted  in  Fig.  5  and  Fig.  6,  which  show 
reasonably  good  agreement.  The  difference  is  probably  due  to  the  difference 
between  the  estimated  and  actual  parasitic  element  values,  particularly  at  the 
higher  frequencies.  The  details  of  the  models  were  presented  in  our 
Technical  Report  [2]  submitted  in  September,  1985. 
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Fig.  5  Comparison  of  measured  (x)  and  calculated  (.)  scattering 
parameters  S 1 1  and  S22  as  a  function  of  frequency,  in  GMz 
VGS  =  -2  V  and  VdS  =  V. 


(b)  Nonlinear  performance  of  the  microwave  GaAs  MESFET 
Microwave  linear  and  nonlinear  performances  for  a  MESFET  have 
been  further  investigated  based  on  our  large  signal  model  for  the  intrinsic 
FET(Fig.  7).  The  nonlinear  dependence  of  the  capacitor  coefficients  on  the 
terminal  voltages,  i.e.,  the  gate  voltage  VGS(t)  and  the  drain  voltage  VDS(t), 
is  calculated  and  includes  the  displacement  current  in  the  intrinsic  model. 

The  small  signal  model  is  not  appropriate  in  this  analysis  since  it  is  based  on 
the  the  assumption  of  small  ac  signal;  a  linearization  procedure  is  used.  The 
signal  response  in  the  time  domain  for  nonlinear  performance  is  simulated 
with  the  two  sinusoidal  signal  of  different  frequencies  shown  in  Fig.  8  but  the 
same  power  shown  in  Fig.  9  as  input  signals.  The  resulting  output  voltage 
will  then  be  obtained  as  in  Fig.  10.  To  study  the  power  gain  and  the  third- 
order  spurs,  the  output  waveform  is  then  transformed  into  the  frequency 
domain  by  Fourier  analysis,  where  the  input  and  output  powers  at  the 
fundamental,  first  order,  second  order,  and  the  third-order  harmonics  are 
calculated  (Fig.  11).  The  power  gain  for  the  fundamental  frequency  and  the 
third-order  spurs  versus  input  power  are  shown  in  Fig.  12  and  Fig.  13, 
respectively.  The  large  signal  response  for  asymmetrical  power  input  is  also 
studied.  For  an  asymmetrical  input  power  spectrum  as  shown  in  Fig.  14,  the 
output  voltages  in  the  time  and  frequency  domains  are  obtained  (Fig.  15  and 
Fig.  16).  The  bias  point,  point  A  in  Fig.  17  in  our  study,  is  chosen  to  have  the 
largest  voltage  swing.  We  have  also  shown  that  the  output  voltages  are 
clipped  off  when  the  bias  points  are  either  near  saturation(point  C)  or  near 
cut-off(point  B)  in  the  I-V  curve.  Fig.  17.  A  detailed  report  on  the  nonlinear 
performance  of  the  MESFET  is  presented  in  our  Technical  Report  [3]. 
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(c)  Modeling  of  the  High  Electron  Mobility  Transistor  (HEMT) 

Since  the  high  mobility  electrons  have  been  observed  at  the  interface  of 
GaAs-AlGaAs  layers,  there  is  considerable  interest  in  the  applications  of 
high-speed  devices  ultilizing  such  interface  electrons.  Among  them,  the 
fastest  transistor  in  the  world  today,  the  high  electron  mobility  transistor 
(HEMT),  holds  the  record  of  having  cut-off  frequencies  up  to  70  GHz, 
switching  speed  less  then  ten  trillionth  of  a  second  (10  psec.).  Other 
advantages  of  HEMTs  are  low-power  consumption  and  a  potentially  simple 
fabrication  process.  We  employed  three  analytical  models,  the  classical 
Fermi-Dirac,  the  two-level  two-dimensional  electron  gas  (2DEG)  and  the 
extended  '^DEG  models,  to  describe  the  interface  electron  concentration  in 
the  AlGaAs  interface  (Fig.  18)  in  an  HEMT  (Fig.  19).  Quantization  effect, 
background  doping  and  triangular-well  approximation  are  investigated.  A 
comparison  of  the  surface  electron  concentrations  due  to  the  three  models  is 
shown  in  Fig.  20.  The  classical  model  takes  all  continuum  energy  levels  into 
account  but  neglects  the  quantization  effect.  Due  to  the  very  high  electron 
concentration  confinement  in  the  very  narrow  heterointerface,  the  two-level 
2DEG  model  is  introduced  to  describe  the  quantization  effect.  However,  this 
model  neglects  all  the  allowed  states  but  the  lowest  two,  which  may  introduce 
some  error.  The  extended  2DEG  model  takes  both  the  quantizing  effect  and 
the  electrons  in  the  upper  levels  into  account  and  is  propably  the  best  model 
among  the  three.  A  simple  numerical  procedure  has  been  developed  to 
obtain  the  surface  electron  concentration  versus  band  bending  and  the  I-V 
characteristic  cur/es  of  HEMTs.  Using  available  parameters,  the  calculated 
I-V  curves  agree  very  well  with  experimental  data(Fig.  21).  The  analytical 
model  presented  here  is  very  useful  for  computer-aided  design  of  ultrahigh- 
speed  GaAs  integrated  circuits  as  well  as  a  microwave  amplifier  using  HEMT 
devices.  A  more  detailed  description  of  our  modeling  is  presented  in  our 
Technical  Report  [4]. 

The  numerical  model  for  HEMT  using  the  Finite-element  Method 
(FEM)  is  also  investigated.  Due  to  the  complexity  of  the  HEMT  device  and 
the  need  of  miniaturization  for  high-.speed  operation,  the  gradual  channel 
approximations  used  in  our  analytical  models  no  longer  appropriate  for 
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Fig.  20  A  comparison  of' the  surface  electron  concentration  versus  the 

normalized  potential  using  three  models;  the  Femii-Dirac  integral 
approach  (solid  line),  the  extended  2DEG  model  (short  dash  line) 
and  the  two-dimensional  electron  gas  (2DEG)  model  (long  dash  line). 
The  unintentional  doping  concentralion  is  1014  cnr-^. 
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submicron  devices.  Furthermore,  the  FEM  model  has  the  flexibility  for  an 
arbitrary  geometry,  which  is  an  extremely  useful  tool  in  simulating  any 
practical  device  performance.  The  major  difficulty  in  HEMT  occurs  in  the 
heterojunction  part  in  which  the  electron  and  hole  concentrations 
accumulated  at  the  heterointerface  as  well  as  the  potential  and  electric  field 
distributions  have  to  be  known  precisely  to  describe  the  device  operation. 
Such  information  has  to  be  obtained  by  solving  the  semiconductor  device 
equations  self-consistently.  The  numerical  model  also  has  the  capability  of 
doing  transient  analysis  which  is  essential  for  high  frequency  analysis  as  well 
as  large  signal  analysis.  Our  first  step  to  attack  this  problem  is  to  solve  the 
semiconductor  equation  self-consistently  for  the  heterojunction  which  is  the 
most  essential  part  of  HEMT  devices.  A  simple  heterojunction  diode  has 
been  chosen  as  our  testing  case  since  it  is  the  simplest  heterostmcture  and  the 
result  is  relatively  well  known.  The  diode  is  divided  into  small  triangular 
elements  as  shown  in  Fig.  22  and  the  band  diagram  is  calculated  as  shown  in 
Fig.  23. 
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Fig.  23  Calculated  band  diagram  of  the  double  heterojunction  diode 
obtained  from  our  FEM  code  with  applied  voltage  Va  =  0.7  V 
foivvard  bias 
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(d)  Quantum  well  Stark  effects  and  their  microwave  applications 


1 .  Quantum  well  Stark  Effects 

Subband  eigenstates  in  an  infinite  quantum  well  with  an  applied 
electric  field  are  calulated  using  both  the  exact  analytic  solutions  (linear 
combination  of  Airy  functions)  and  the  variational  trial  functions 
orthogonalized  by  the  Gram-Schmidt  approach.  The  shifts  of  the  subband 
energies  due  to  the  electric  field  are  different  between  the  ground  state  and 
higher  energy  states.  For  increasing  electric  fields,  the  ground  state  shows  a 
large  negative  energy  shift,  while  higher  states  have  small  positive  shifts  for  a 
moderate  range  of  electric  field.  Quasi-bound  states  of  a  finite  quantum  well 
with  an  electric  field  have  been  calculated  by  solving  analytically  one-particle 
quantum  well  Stark  resonance  .  Resonance  position  and  life  time  of  a  quasi¬ 
bound  states  have  been  estimated  from  tlie  real  and  the  imaginary  parts  of  the 
complex  energy  eigenvalue,  respectively. 

2.  Intersubband  optical  absorption  with  an  applied  electric  field 

Tlie  existence  of  discrete  energy  levels  by  quantum  confinement  of 
carriers  in  a  semiconductor  quantum  well  results  in  resonance  optical 
transitions  of  infrared  radiation  of  suitable  polarization  between  the 
quantized  subbands.  The  electric  field  dependence  of  the  intersubband  optical 
absorption  has  been  calculated  for  parabolic  band  and  infinite  well  model. 
Effects  of  the  finite  well  are  taken  into  account  by  choosing  an  effective  well 
width  which  gives  the  same  ground  state  energy.  The  Stark  shifts  of 
intersubband  optical  absorption  spectra  observed  experimentally  agree  with 
our  theoretical  calculations  of  the  intersubband  electro-optical  absorption 
(Ahn  and  Chang,  Phys.  Rev.B.35,4149,1987).  From  our  early  calculations  of 
quantum  well  Stark  effects,  it  was  predicted  that  the  intersubband  peak 
absorption  would  show  blue  shifts  in  energy  with  an  increasing  electric  field. 
For  the  peak  absorption  amplitude,  two  different  experimental  results  were 
reported,  however.  Harwit  and  Harris(Appl.  Phys.  Lett.  50,685.1987) 
observed  that  the  peak  amplitude  decreased  with  an  increasing  electric  field 


and  Bajema  et  al.(phys.Rev.B.36,l  300,1987)  showed  that  that  the  peak 
amplitude  increased  with  an  increasing  electric  field  , which  agree  with  our 
early  calculation  in  which  the  Fermi  level  is  assumed  to  be  constant  (  or 
variation  of  the  Fermi  level  is  smaller  than  that  of  the  energy  level ).  Both 
experiments  show  the  blue  shifts  of  the  peak  energy,  which  agree  with  our 
results.  It  is  well  known  that  the  absorption  coefficient  depends  on  the  dipole 
moment,  the  density  of  states,  and  the  linewidth.  For  intersubband 
absorption,  the  dipole  matrix  element  varies  very  slowly  with  an  electric 
field.  Therefore  a  major  change  of  the  absorption  peak  could  be  attributed  to 
the  change  of  the  poulation  and  the  linewidth  broadening,  which  may 
strongly  depend  on  the  condition  of  the  experiment.  Preliminary  results  for 
the  intersubband  transitions  are  obtained  for  two  different  hypothetical 
conditions:  (1)  the  variation  of  Fermi  level  is  smaller  than  the  energy  level 
change,  and  (2)  the  Fermi  level  change  is  comparable  to  or  greater  than  the 
energy  level  change  so  that  total  electron  density  in  the  well  remains 
constant.  We  found  that  the  peak  absorption  increases  with  field  for  case  (1) 
and  decreases  for  case  (2).  The  linewidth  broadening  if  included  will  further 
reduce  the  absorption(Ahn  and  Chuang, Superlattices  and  Microstructure, in 
press). 

3.  Nonlinear  optical  properties  of  intersubband  transition  with  an 

electric  field 

Recent  expermental  studies  by  West  and  EgIash(Appl.  Phys. 

Lett.46, 1156,1984)  show  that  the  intersubband  optical  transitions  have  a  very 
large  dipole  moment.  This  suggests  that  the  intersubband  transitions  in  a 
quantum  well  may  have  very  large  optical  nonlinearities.  We  have  calculated 
the  third-order  nonlinear  absorption  coefficient.  The  saturation  begins  to 
occur  at  I  =  1 .0  MW/cm^  for  the  well  width  L  =  126.5  A,  for  a  moderate 
range  of  electric  field.  We  also  have  studied  the  possibility  of  tlie  electric- 
field  control  of  the  second  harmonic  generation  in  a  quantum  well .  The 
asymmetry  of  the  quantum  well  due  to  the  electric  field  accounts 
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for  the  non-vanishing  of  the  second-order  susceptibilities.  For  a  moderate 
value  of  electric  field  of  10  kV/cm  to  70  kV/cm,  the  second  harmonic 
susceptibility  is  generally  10  to  100  times  larger  than  that  of  bulk  GaAs. 

4.  Electric-field  effects  on  the  polar  optical  phonon  scattering  in  a 

quantum  well 

To  study  the  electric-field  effects  on  the  linewidth  broadening 
mechanisms,  we  have  estimated  the  field  dependence  of  two  important 
linewidth  broadening  mechanisms:  (1)  polar  optical  phonon  scattering  and 
(2)  tunneling  of  quasi-bound  states.  Our  calculation  of  polar  optical  phonon 
scattering  using  a  simple  model  shows  that  the  scattering  rates  are  enhanced 
in  general  with  an  increasing  electric  field.  Electric-field  effects  are  more 
prominent  for  the  heavy  hole  than  for  the  electron.  Higher  scattering  rates 
and  stronger  field-dependences  of  the  heavy  hole  can  be  attributed  to  its 
heavier  effective  mass  compared  to  that  of  the  electron. 

5.  Gain  control  of  GaAs-AlGaAs  lateral  current  injection  quantum  well 

lasers  with  electric  field 

Quantum  well  lasers  are  of  growing  interest  partly  because  of  their 
superior  characteristics,  such  as  low  threshold  current,  narrow  gain 
spectrum,  and  less  temperature  dependence,  compared  with  those  of 
conventional  lasers.  Recently,  a  new  GaAs-AlGaAs  lateral  current  injection 
quantum  well  laser  was  proposed  in  which  the  stmcture  is  planar  suitable  for 
integration  with  other  optical  devices  and  exhibits  built-in  index  guiding  and 
a  very  low  stray  capacitance.  We  have  calculated  the  electric-field  effects  on 
the  gain  spectra  of  a  lateral  current  injection  quantum  well  laser.  Simple 
parabolic  bands  are  used  in  the  calculation.  An  applied  electric-field  will 
push  the  electrons  and  holes  to  opposite  sides  of  the  well  in  the  direction 
perpendicular  to  the  well,  tlierefore  the  overlap  integral  for  the  optical 
transition  will  decrease.  Both  the  electron  and  the  hole  ground-state  energy 
levels  become  lower  with  increasing  electric  field,  thus  the  gain  peak  is 
expected  to  be  shifted  to  lower  energy.  Our  calculation  shows  that  the  peak 
gain  exhibits  a  red  shift  and  the  peak  gain  amplitude  decreases  with  an 


increasing  electric  field.  It  is  found  that  the  TE  mode  is  more  affected  by  the 
applied  field  than  the  TM  mode.  These  results  are  significant  for  applications 
to  the  direct  microwave  modulations(publication  no.  12). 
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field:  Gram-Schmidt  orthogonalization  approach 
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We  present  variational  calculations  of  subband  eigenstates  in  an  infinite  quantum  well  with  an 
applied  electric  field  using  Gram-Schmidt  orthogonalized  trial  wave  functions.  The  results 
agree  very  well  with  the  exact  numerical  solutions  even  up  to  1200  kV/cm.  We  also  show  that 
for  increasing  electric  fields  the  energy  of  the  ground  state  decreases,  while  that  of  higher 
subband  states  increases  slightly  up  to  1000  kV/cm  and  then  decreases  for  a  well  size  of  1()0  A. 


Electronic  as  well  as  optical  properties  of  quantum  wells 
subject  to  external  electric  fields  have  received  much  atten¬ 
tion.'  ''  These  areas  are  interesting  both  from  a  fundamental 
and  a  practical  point  of  view.  The  possible  device  applica¬ 
tions  include  the  use  of  quantum  confined  Stark  effect’  in 
optical  modulators"  and  optical  switching  devices.”  As  yet. 
most  of  the  theoretical  works  has  been  confined  to  the  calcu¬ 
lations  of  the  ground  state.  Very  recently.  .Matsuuraand  Ka- 
mizato'  reported  an  exact  numerical  calculation  of  subbands 
in  an  infinite  well  and  concluded  that  the  higher  subbands 
behave  very  diffei»*nlly  from  that  for  the  ground  state  when 
the  electric  field  strength  is  increased.  Exact  solutions  em¬ 
ploying  two  independent  Airy  functions''  are  sometimes 
too  complicated  to  use  in  real  problems.  The  variational  ap¬ 
proach  has  the  advantage  of  providing  analytical  expres¬ 
sions  for  the  eigenstate  energies  and  the  wave  functions,  and 
numerical  results  with  reasonable  accuracy  can  be  obtained. 
Analytic  form  of  the  trial  wave  function  for  the  ground  states 
has  been  known.'  '  hut  no  useful  forms  of  the  trial  solutions 
for  the  higher  subbands  which  yield  accurate  results  com¬ 
pared  with  the  exact  ones  are  given  yet.  Recently,  a  trial 
wave  function  of  the  form  y  „  (z)exp(  ~  P,  z)  for  the  «th 
subband  has  been  suggested.'  where  f  „  iz)  is  the  zero-field 
nth  quantum  well  bound  state  wave  function  and  P,  is  a 
variational  parameter.  However,  it  is  pointed  out  that  the 
solutions  may  yield  very  differeni  results  from  the  exact  ones 
( for  n  >  I  )  because  these  trial  functions  arc  not  orthogonal 
to  each  other  That  observation  is  also  confirmed  in  this 
paper  and  the  numerical  results  are  illustraled.  Thus,  it  is 
important  to  find  an  orthogonalized  set  of  Inal  wave  func¬ 
tions  for  the  variational  calculations  ofsubband  energies  and 
wave  functions. 

We  report  in  this  letter  variational  calculations  on  sub¬ 
band  states  in  an  external  uniform  electric  field  based  on  the 
infinite  well  model.  We  find  analytic  forms  of  orthogona- 
lized  trial  wave  functions  by  the  Green-Schmidt  orlhogona- 
lization  procedure.  Our  calculations  agree  very  well  with 
exact  numerical  results  up  to  120<)  kV/cm  with  an  error  less 
than  90  for  L  -  l()<)  .A  f  or  comparison,  we  also  show  the 
numerical  result'  for  which  the  trial  wave  functions  are  not 
orthogonalized 

It  IS  well  known  ihal  ihe  vanalion.il  nielhod  can  also  he 


used  to  obtain  one  of  the  higher  energy  levels  if  the  function 
is  orthogonal  to  the  eigenfunctions  of  all  the  lower  states.'" 
The  most  common  method  of  obtaining  an  orthogonal  set  of 
functions  is  the  Gram-Schmidt  orthogonalization  proce¬ 
dure."  which  is  the  construction  of  an  orthonormal  set 
{^,.<3.....}  from  a  finite  or  an  infinite  independent  set 
{i/i.u ,,...}  wliich  is  not  necessarily  orthonormal. 

Let  us  consider  an  electron  with  charge  —  |e|  and  effec¬ 
tive  mass  m*.  in  an  infinite  quantum  well  width  L  in  the 
presence  of  a  constant  electric  field  F  along  the  positive  di¬ 
rection  of  the  well  z.  We  choose  the  origin  to  be  at  the  center 
of  the  well.  The  Hamiltonian  of  the  system  in  the  effective 
mass  approximation  is  given  by''’ 

//  =  //„  4- Wfz.  (1) 

where  H„  is  the  unperturbed  Hamiltonian  whose  eigenvalues 
are  given  by 

=  (fftr/lm*L-)n-  zt=1.2....  (2) 

For  our  specific  problem,  we  chose  the  nth  vector  u,  to  be 


(z)  =  sin 

xe.xp 

which  IS  not  an  orthogonal  set.  where/?,  is  the  nth  variation¬ 
al  parameter  One  can  easily  .see  that  for  n  =  1,  u,  is  the  trial 
solution  introduced  by  Bastard  et  al.'  We  also  define  the 
inner  product  between  two  functions  f  and  g,  (f'g).  by 

</^^)sJ  ^  /*lz)g(z)dz,  (4) 

where  the  superscript  *  denotes  a  complex  conjugate. 

The  prixiedure  we  use  is  as  follows; 

n  -  I.  ( I )  Let  i/'i  -  M I  .  ( 5a ) 

liil  Minimize  £,(/?, )  -  {f.  H  d'|).(5b) 


and  find  P. 
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and  find  ^2’ 

(iii)  <^2  =  0;/<«f'2|l/'.>'''  • 

/  I 

n  =  l,  (i)  Ltlil>,  =  u,  -  ^  {ui\<^,)<t),  . 


(ii)  Minimize  £,(/?,)  =  <<f’,)// |i/», )/<(/', |(/^;)  ,  (7b) 
and  find /?/. 

(7a)  (m)  d),  =  il/,/(itii\il/i)''- .  (7c) 


This  procedure  requires  mainly  three  inner  products  which  are  in  analytical  forms: 
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ir  +  {P„+P„)^  {n  +  m)-ir  +  (P„  +  P„)- 

l-(-l)'  "'exp[  -(^,  +^„)] 
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X{I  -  (  -  1)-'  '"exp[  -  (P„  ^p^)J}  +  (P^  ^p^)(  -  I)"’  "expl  -(/?,  +P„)]\ 
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The  minimization  of  an  analytical  function  and  finding 
)7"'"  can  be  directly  done  in  the  computer  by  calling  some 
subroutines,  e.g  ,  using  International  Mathematical  Statisti¬ 
cal  Libraries  (IMSL)  subroutines. 

For  the  expectation  value  of  the  ground-state  energy 
£i(/Si),  we  have 

£i()3,)  = 


where  the  parameter  £is  defined  by 

£=  :c:££ /£','”.  tl2) 

which  turns  out  to  be  the  normalized  electric  field.’ 

In  the  low  field  limit,  F <  1, 


I 

lU, 

II 

E- 

(  13) 

and 

A£,  ~-F,(p'r>  -  £'i"'  — 

n  V  3 

2  y  >rt\  'F  'L  ‘ 
r'J  ”  fi 

1  141 

In  the  high  field  limit.  /•  '•.>  1 

PT"'  -  (,r'F)'  ' 

(15) 

and 


Ai,  .  -  -,,ia  +  (-)  (_)  . 


(16) 


Equations  (11)-(16)  have  already  been  obtained  by 
Bastard  ei  al. '  except  they  missed  the  factor  of  1  /2  in  front  of 
P  in  Eq.  ( 7 )  and  the  expressions  for  P  7'"of  Ref  1  in  the  low 
and  high  field  limits  are  not  correct.  This  does  not  affect  their 
finai  results  for  A£|,  which  are  correct.  Calculations  of  high¬ 
er  energy  levels  are  straightforward  following  the  proce¬ 
dures  in  Eqs.  ( 5  )-(7)  The  results  of  the  normalized  energy 
£„  =  £„/£','”  for  the  first  three  states  are  plotted  versus  the 
normalized  electric  field  £  in  Fig.  1.  The  plot  in  terms  of 
these  two  normalized  parameters  is  universal  which  can  be 
seen  from  the  original  Schrddinger  wave  equation.  It  is 
readily  seen  that  the  shift  of  the  subband  energy  due  to  the 
electric  field  is  different  between  the  n  =  1  state  and  the 
higher  energy  states  For  increasing  electric  fields,  the 
ground  state  shows  a  large  negative  shift,  while  higher  states 
have  small  positive  shifts  for  fields  up  to  F^IQ  and  then 
negative  shifts  We  have  also  plotted  the  subhand  energies 
obtained  from  the  e.xact  numerical  solution  We  find  that 
both  methods  ol  calculation  give  verv  close  results,  even  up 
to  F~  25.  The  parameters  we  use  are  m*  ~  OOb.s/n,,  and 
/-  -  1(X)  A  where  in,,  (leiiotes  ihe  free-electron  mass  F  .  Id 
corresponds  to  s's  s  |,,V  emfi'i  :his/.  The  results  ol  calcu¬ 
lations  employing  unorthogoriali/ed  iri.il  functions'  defuievl 
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FIG.  1.  Normalized  subband  energies  (£.  =  £,/£'“’)  are  plotted  vs  the 
normalized  electnc  field  (F  =  lel/X  ( 1 )  The  enact  solution —  (sol¬ 

id  lines).  (2)  the  variational  method  with  Gram-Schmidt  orthogonaliza- 

tion  procedure - (dashed  lines),  (3)  the  vanational  method  (see Ref.  5) 

without  orthogonalization . (dotted  lines). 
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FIG.  2.  Magnitudes  of  the  wave  functions,  L  |tf,  (r)l’,  are  plotted  vs  the 
normalized  distance  z/L  for  the  first  three  subbands.  The  dotted  lines  are 
for  the  zero  electric  field.  The  solid  lines  are  the  exact  solutions  and  the 
dashed  lines  are  those  for  the  variational  method  with  the  Gram— Schmidt 
onbogonalization  procedure  for  a  quantum  well  with  an  applied  electric 
field  (£=:=  20). 


in  Eq.  (3)  are  plotted  in  Fig.  1  (dotted  lines).  One  readily 
notices  that  these  solutions  begin  deviating  significantly 
from  the  exact  ones  at  around /■  =  1.5  (/■=  86.8  kV/cm  for 
L  =  100  A).  In  Fig.  2  we  plot  the  square  of  the  wave  func¬ 
tions  1^1^  for  the  first  three  states  when  the  electric  field  is 
F"  =  20.  is  normalized  to  L.  We  see  in  Fig.  2  that  the 
ground-state  wave  function  is  shifted  to  z  <  0  region  signifi¬ 
cantly,  While  higher  subband  wave  functions  still  have  a 
large  amplitude  in  the  z  >  0  region  even  at  F  =  20. 

In  conclusion,  we  have  derived  orthogonalized  trial 
wave  functions  which  yieid  very  close  results  to  those  of  the 
exact  numerical  solutions.  It  is  shown  that  for  decreasing 
electric  fields  the  energy  of  the  lowest  subband  decreases, 
while  that  of  the  higher  subbands  increases  slightly  up  to 
F=  20  (FZ  lOOOkV/cm  for  L  =  100  A). 
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Exact  calculations  of  quasibound  states  of  an  isolated  quantum  well  with  uniform  electric  field; 

Quantum-well  Stark  resonance 
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We  present  universal  plots  from  exact  numerical  calculations  for  the  energy  level  and  the  reso¬ 
nance  width  of  quasibound  states  in  a  quantum  well  with  an  applied  electric  field  (quantum-well 
Stark  resonance)  by  solving  the  Schrddinger  equation  directly.  This  calculation  gives  both  the 
resonance  positions  and  widths  for  the  complex  eigenvalue  £'0  — iT/2  of  the  system.  Our  theory 
also  shows  that  the  energy  shifts  of  the  ground  states  for  the  electrons  and  holes  have  the  same 
behaviors  in  high  fields  without  any  turnaround  phenomenon,  contrary  to  the  results  of  Austin  and 
Jar  os. 


Electronic  and  optical  properties  of  quantum  wells  with 
applied  external  electric  fields  are  of  increasing  interest. 
Studies  of  these  areas  are  important  both  from  a  funda¬ 
mental  and  a  practical  point  of  view.  Optical  modulators' 
and  optical  switching  devices^  based  on  the  quantum  con¬ 
fined  Stark  effect  have  been  suggested.  Possible  device  ap¬ 
plications  of  the  field-induced  tunneling  in  quantum-well 
and  quantum-barrier  hetcfostruclurcs  include  high-speed 
resonant  tunneling  devices. 

More  recent  theoretical  studies’"’  of  the  effects  of  exter¬ 
nal  electric  fields  on  the  quantum-well  systems  have 
predicted  both  the  field-induced  level  shifts  and  the  field 
dependence  of  the  carrier  lifetime.  In  this  paper,  we  report 
exact  numerical  calculations  on  quasibound  states  of  a 
quantum  well  in  an  external  electric  field  (quantum-well 
Stark  resonance)  by  solving  the  Schrddinger  equation  for 
Stark  resonance  directly.  It  is  found  that  the  previous  re¬ 
sults  based  on  phase-shift  analysis’’  and  the  stabilization 
method*  agree  very  well  with  our  results  over  a  wide  range 


PoieniiQi  Energy 
V(2) 


Slope  =Ie|F 


of  the  electric  field.  At  an  extremely  high  electric  field, 
there  is  no  turnaround  behavior  in  the  energy  shift  for  both 
the  electrons  and  the  holes,  contrary  to  the  results  in  Ref. 
7,  where  no  explanation  can  be  provided  for  that 
phenomenon.  We  believe  that  our  direct  numerical  ap¬ 
proach  is  very  reliable  even  at  a  very  high  electric  field, 
while  the  results  using  the  phase-shift  analysis  may  have 
drawbacks  in  the  high-field  limit.  Our  approach  has  an 
advantage  over  the  previous  results’"’  in  that  both  the 
Stark  resonance  position  (quasibound-state  level)  and  the 
width  can  be  obtained  from  the  single  complex  energy 
eigenvalue  of  the  quantum-well  Stark  resonance  problem. 
The  disadvantage  is  that  numerical  subroutines  of  the  Airy 
functions  with  complex  arguments  are  required. 

Consider  an  electron  with  charge  —  \  e\  and  effective 
mass  m*,  in  a  finite  quantum  well  of  width  L  and  depth  Kq 
in  the  presence  of  a  constant  electric  field  F  along  the  posi¬ 
tive  direction  of  the  well  2  (Fig.  l).  We  choose  the  origin 
to  be  at  the  center  of  the  well.  The  Schrddinger  equation 
of  the  system  in  the  effective-mass  approximation  is  given 
by’-’ 


2m*  dz^ 


\p(z  )  -I-  I  e  I  Fztpiz  )  ’“Eipiz  ), 


</,/2  , 


(1) 


2m*  dz^ 


tl/(,z  )  -h  (Fq-I-  1  c  I  Fc  Itplz  )  =E\i/{z  ), 


I2  I  >£/2  . 


FIG.  1.  Potential-energy  profile  V(z)  for  a  single  quantum 
well  with  depth  Fo  and  width  L  subject  to  an  external  electric 
field  F. 


Since  the  potential  energy  term  in  Eq.  ( 1 )  lends  to  —  <»  as 
2  goes  to  —  00,  the  system  does  not,  strictly  speaking,  have 
true  bound  states. In  other  words,  the  particle  initially- 
confined  in  a  weil  can  always  lower  its  potential  energy  by 
tunneling  out  of  the  well  when  the  field  is  not  zero.  It  may 
happen,  however,  that  the  tunneling  probability  is  very 
small.  In  such  a  case,  we  can  regard  the  system  as  having 
ouasibound  stales,  in  which  the  particles  move  “inside  the 
well''  for  a  considerable  period  of  lime  and  leave  through 
tunneling  only  when  a  fa'rlv  long  time  interval  r  has 
elapsed.  In  discussing  the  quasibound  slates,  we  m  \y  use 
the  following  formal  method.  Instead  of  -.-onsidering  the 
solutions  of  the  Schrddinger  equation  with  a  boundary 
condition  requiring  the  finiicness  of  the  wave  func'ion  at 
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infinity,  we  shall  look  for  solutions  which  represent  outgo¬ 
ing  waves  at  infinity;'®  this  implies  that  the  particle  finally 
leaves  the  well  by  tunneling.  Since  such  a  boundary  condi¬ 
tion  is  complex,  we  cannot  assert  that  the  eigenvalues  (en¬ 
ergy)  must  be  real.  By  solving  the  Schrddinger  equation, 
we  obtain  a  set  of  complex  eigenvalues,  which  we  write  in 
the  form 

f-Eo-tT/2,  (2) 

where  f  is  found  to  be  pwsitive.  Eq  and  f  correspond  to 
the  quasibound-state  energy  level  and  the  resonance  width, 
respectively.  The  tunneling  probability  per  unit  time  is  de¬ 
fined  by 


The  solutions  to  Eq.  (1)  with  the  outgoing-wave  condition 
are  linear  combinations  of  two  independent  Airy  func- 


flilBi(i72)  +  i  Ai(t72)l.  z  <  “E/2. 
aoAi(r;i)  +  6oBi(rj(),  |z  |  :<E/2, 
fl2Ai(rj2))  -z  >  “E/2  , 


2m*  I 

(ehFy  J 


<5b) 

The  wave  function  for  z  <  “E/2  represents  an  electron 
traveling  to  z  “  —  »  after  tunneling.  The  complex  energy 
E  can  be  found  by  solving  the  secular  equation  obtained  by 
matching  the  value  of  ig  and  its  first  derivative  at  the 
points,  z  ~  ±E/2.  The  resulting  determinantal  equation 
is 


Ai(f)i''')  Bifqi*')  “Ai(fj2^) 
Ai'(r„+)  Bi'frj,-'')  “Ai'frjj*') 
Aifrjr)  Ci(qr)  0 
Ai’fqf)  Bi'(qf)  0 


“  [Bi(rj2  )  +f  Aifqi  )1 
—  [Bi'(rj2~  )  +  <  Ai'(r;2” 


where  ric  and  T]-r  are  the  values  of  rji  and  rji  evaluated  at 
z  “■E/2  and  “E/2,  respectively.  If  we  introduce  a  new 
parameter  E'®’ defined  by 

f  rr  I 


(v.hich  happens  to  be  the  ground-state  energy  of  an  infi¬ 
nite  quantum  well  with  width  E ),  and  define  the  normal¬ 
ized  energy  £”E/E'®’,  the  normalized  electric  field 
E  =*  I  e  I  EE/E'®',  and  the  normalized  well  depth  Lq 
”Ko/E'®’,  we  may  express  r/r  and  rir  by  these  three 
normalized  quantities;  E,  E,  and  Vq. 

rir--\-^]  (E  +  tE)  ,  (8a) 


This  means  that  the  solution  of  E  from  Eq.  (6)  is  universal 
and  can  be  used  for  both  electrons  and  holes  with  the  re¬ 
placement  of  the  parameter  E'®’  with  their  corresponding 
elfective  masses  (Here  the  effective  masses  inside  and 
outside  the  well  are  assumed  to  be  equal.)  The  nuiT.Mlized 
energy  E  can  be  expressed  in  terms  of  only  two  normalized 
parameters,  I'o  and  E.  Thus  it  is  clear  that  both  electrons 
and  holes  should  have  the  same  behaviors  in  their  energy 
shift  and  the  resonance  width.  To  obtain  the  results  of  Eq 


TABLE  I.  Comparison  of  the  numerical  results  for  Eo“iT/2 
using  the  exact  numerical  method  of  this  paper,  the  phase-shift 
analysis  (Refs.  7  and  9),  and  the  stabilization  method  (Ref.  8). 
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FIG.  3.  The  real  part  of  the  normalized  energy 
for  various  normalized  well  depths  is  plotted  vs  the 

normalized  electric  field  F  “  |  e  |  FL/E^°\ 


and  r  for  holes,  one  need  only  multiply  E  by  Eq  using  the 
effective  mass  of  the  hole.  We  have  solved  Eq.  (6)  numer¬ 
ically  to  the  desired  accuracy  using  the  series  and  asymp¬ 
totic  expansions  of  the  Airy  functions  with  complex  argu- 


To  check  the  validity  of  our  approach,  we  com¬ 


pared  our  results  with  those  of  the  previous  methods’*  in 
Table  1.  The  values  of  Vq,  L,  and  m*  for  the  heavy  holes 
used  in  the  calculations  are,  respectively, 

Ko”  100  meV,  L  -37  A,  m*  >=0.45mo  ,  (9) 


“0.45mo  , 


where  nto  is  the  free-electron  mass.  It  is  readily  seen  that 
our  results  agree  very  well  with  those  of  the  phase-shift 
analysis’'^  and  the  slibili/.ation  method.* 

In  Fig.  2,  the  real  part  of  the  energy  Eo  (resonance 
position- solid  line)  for  the  ground-state  energy  with  the 
values  of  Fg,  L,  and  w*  for  electrons  given  by  Kg “340 
meV,  L  “  100  A,  and  m*  “0.0665mg  is  compared  with  the 
results  of  infinite-well  variational  calculations**”'^  (dashed 
lines),  where  we  have  used  an  effective  well  width 
Z.eff“  126.5  A,  cho.sen  to  give  the  same  £g  at  zero  field  for 
the  variational  calculations.  It  can  be  easily  noticed  that 
both  calculations  gave  very  similar  results  even  up  to 
2x  10*  V/cm.  However,  the  variational  calculation  for  the 
infinite-well  model  cannot  give  the  resonance  width  since 
no  tunneling  exists  for  the  infinite  well.  The  results  of  the 


are  plotted  versus  E  in  Fig.  .3.  In  contrast  to  the  previous 
results’  which  are  still  controversial,  the  resonance  position 


normalized  field  F.ieiFL/E” 


FIG.  4.  The  normalized  resonance  width  f  for  vari¬ 

ous  Fg  is  plotted  vs  normalized  elecric  field  F  —  1  e  |  F/,/£“”. 


is  found  to  be  in  the  well  even  at  very  high  field.  The 
behaviors  of  the  resonant  position  are  the  same  for  both 
electrons  and  holes  with  proper  f*®*  used  together  with 
Fig.  3  as  discussed  before.  Thus  the  turnaround  behavior 
for  the  holes  and  electrons  in  the  energy  shift  shown  in 
Ref.  7  is  probably  a  drawback  of  that  method  itself.  IJ.sing 
the  same  numerical  values  for  holes  as  those  in  Ref.  7, 
£“30  A,  Kg“70  meV,  m*  “0.45mg,  we  obtain 
“92.26  meV,  Fg—0.76.  We  do  not  have  any  turnaround 
behavior  even  up  to  F“10,  or  the  electric  field  £—3075 
kV/cm,  wiiiLii  CO. CIS  a  mucii  wiJci  lu.i^e  of  electric  field 
than  that  of  Ref.  7.  In  Fig.  4,  we  plot  the  normalized  reso¬ 


nance  width  F' 


for  various  Fg  vs  £.  Since  the  life¬ 


time  r  is  defined  by  r”/t/r,  the  results  plotted  in  Fig.  4 
predict  a  rapid  decrease  of  the  carrier  lifetime  with  in¬ 
creasing  applied  field  by  field  enhanced  tunneling. 

In  conclusion,  we  have  solved  the  Schrodinger  equation 
for  a  quantum  well  with  uniform  electric  field  directly. 
Complex  eigenvalues  for  quantum-well  Stark  resonance 
are  obtained.  Our  approach  has  an  advantage  over  previ¬ 
ous  analyses’"’  in  that  both  the  resonance  position  and 
width  can  be  obtained  from  a  single  complex  eigenvalue  of 
the  problem. 


normalized  resonance  energy  £g“£g/£  for  various  Vq 
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Intersubband  optical  absorption  in  a  quantum  well  with  an  applied  electric  field 
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We  present  new  results  for  the  electric  field  dependence  of  the  intersubband  optical  absorption 
within  the  conduction  band  of  a  quantum  well.  We  show  that  for  increasing  electric  field  the  ab¬ 
sorption  peak  corresponding  to  the  transition  of  states  I  —  2  is  shifted  higher  in  energy  and  the 
peak  amplitude  is  increased.  These  features  are  different  from  those  of  the  exciton  absorption.  It 
is  also  found  that  the  transition  1  — »  3,  forbidden  when  F~0,  is  possible  when  F  is  nonzero. 


Quantum  confinement  of  carriers  in  a  semiconductor 
quantum  well  leads  to  the  formation  of  discrete  energy  lev¬ 
els  and  the  drastic  change  of  optical-absorption  spectra.' 
The  interband  absorptions  near  the  band  gap  have  been 
extensively  studied,  and  it  has  been  shown  that  their  ab¬ 
sorption  and  luminescence  spectra  are  dominated  by  exci- 
tonic  effects. Some  more  recent  studies  have  concen¬ 
trated  on  the  electric  field  dependence  of  energy  levels'*"’ 
and  band-edge  optical  absorption  including  the  exciton 
effect.*"'®  Very  recently,  experimental  studies  of  the  in¬ 
tersubband  absorption  within  the  conduction  band  of  a 
GaAs  quantum  well  without  an  applied  electric  field  have 
been  reported."  A  very  large  dipole  strength  and  a  nar¬ 
row  bandwidth  were  observed.  In  this  paper,  we  present 
theoretical  calculations  for  the  electric  field  dependence  of 
the  optical  absorption  between  the  discrete  subbands 
within  the  conduction  band  of  a  quantum  well  based  on 
the  infinite-potential-barrier  model.  One  of  the  reasons  for 
increased  interest  in  this  area  is  the  possibility  of  practical 
device  application.  For  example,  in  1970,  Kazarinov  and 
Suris'’  proposed  a  new  type  of  infrared  laser  amplifier  us¬ 
ing  the  intersubband  transition  and  resonant  tunneling.  A 
far-infrared  photixietecior  with  high  wavelength  selectivity 
based  on  the  intersubband  absorption  and  the  sequential 
resonant  tunneling  has  also  been  suggested.  ” 

The  Hamiltonian  of  the  system  (a  single  quantum  v/ell) 
subject  to  a  uniform  electric  field  perpendicular  to  the 
quantum  well  (the  z  direction)  in  the  presence  of  optical 
radiation  (Fig.  1)  is  written  as 

+  ,  (I) 

where  I/q  is  the  unperturbed  Hamiltonian  for  an  electron 
in  the  quantum  well  in  the  piesence  of  perpendicular  elec¬ 
tric  field,  and  the  interaction  Hamiltonian  is  given 
hy'^ 

//;,p--—  A  p=----  .-lo[e'‘'*‘'“"  +  cc.l6  p  .  C) 

Wifi  2wi() 

where  A  i.s  the  vector  potential,  e  is  the  polari/.ition  vector, 
q  IS  the  wave  vcctiir  for  incoming  t'ptical  radiation,  c  is  the 
magnitude  of  the  charge  of  the  electron,  wta  is  the  free- 
spaee  electron  mtiss.  and  p  is  the  momentum  vector  ol  the 
electron  in  the  crystal  The  first  term  in  (21  gives  the  ab¬ 
sorption  and  the  second  term  gives  the  emission  of 

photon 


Then,  for  a  given  interaction  potential  Hop,  the  transi¬ 
tion  rate  from  the  initial  state  i/r,  to  the  final  state  ip/  for 
absorption  is  given  by'^ 

FF/i  “  I  <ig/  I  Wof  I  I  ^S(E/  -Ei~ha>)  ,  (3) 

n 

where  £,  and  £/  are  the  energies  of  the  electron  in  the  ini¬ 
tial  state  and  the  final  state,  respectively,  and  tu  is  the  an¬ 
gular  frequency  of  the  incident  photon.  If  we  neglect  the 
interaction  between  the  electrons  in  the  well,  the  wave 
functions  for  the  initial  state  ig,  and  the  final  state  i///  after 
absorption  can  be  written  as'* 

1//,  “Ui.(r)4,(r) 

"“A  \z\<kf  ■  (“fa) 

V'/“Wf'(r).J/(r) 

“zl  “‘^’uf'(r)e''‘' ''i5/(c ),  Ul<y  ,  (4b) 

where  A  is  the  area  of  the  well.  L  is  the  width  of  the  well, 
k,,  kf'  are  the  wave  vectors  of  the  electron  in  the  .x-p  plane 
for  the  initial  and  the  final  states,  respectively,  r,  is  the  po- 

PotentiQi  t 


FIG.  1.  Polcnti.tl-energy  profile  (or  ,tn  infinuc  qu.intuni  well 
with  width  1.  suhjtcl  to  an  cvlern.il  cicciric  field  /■'  in  lire  pres¬ 
ence  of  incoming  radi.illon  with  .mgufir  Irequencv  hio. 
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sition  vector  in  the  .v-y  plane,  and  Uc  and  Uc'  are  the  cell 
periodic  functions  near  the  conduction-band  extremum. 
The  envelope  functions  0,  and  0/  satisfy  the  following 
Schrddinger's  equation  in  the  effective-mass  approxima- 


where  the  cell  periodic  function  part  has  been  taken  care 
of  as  in  Ref.  16,  and  we  have  used  the  dipole  approxima¬ 
tion.  Since 


tion; 


6=  cf' 


2m*  dz- 


7  0(r ) -1- 1  e  I  Fzpfz  )  =£0(2 ),  jzj 


_£ 

2 


(5) 


^t/2 

/  <>/•(-)(!>/(- )^/z=0/,  , 


and  are  given  by  the  linear  combination  of  two  indepen¬ 
dent  .Airy  functions  Aifrj)  and  Bi(r;),  where  rj  is  defined 
bv 

(/r-IHf.-).  (6) 


In  Eqs.  (5)  and  (6),  m*  and  £  denote  the  effective  mass  of 
an  electron  and  the  electric  field,  respectively. 

For  intersubband  transitions,  the  matrix  element 


we  find  that  the  major  contribution  to  the  optical  matrix 
element  is  the  z  component  of  the  r  vector,  and  so  the  ab¬ 
sorption  is  strongly  polarization  dependent.  The  absorp¬ 
tion  constant  a  in  the  well  is  defined  as'^  ho)  times  the 
number  of  transitions  per  unit  volume  per  unit  time  divid¬ 
ed  by  the  incident  power  per  unit  area 


i/r/ 1  /fop’*  I  Vd  t-'afi  he  approximated  by ' 


^  >•  /  k,  k,' 


(8) 


V/|//of  1^,) 


'?/ !  "op 

eAo  . 


2mo 

eAp 

2ih 


(£,  -£/)fi'<^/|rU,>  , 


(7) 


where  the  summations  over  i  and  /  are  for  the  quantized 
initial  and  final  energies,  respectively,  for  the  2  com¬ 
ponents  of  the  momenta.  If  we  calculate  the  total  transi¬ 
tion  rate  and  take  into  account  the  line  broadening,'^  we 
obtain 


a** 'Ll.  V' (cos  ^  e )  I  A//,- 1  ^In 
/  j  nn  m^Uird) 


+  exp 


ksT 


j  1  +exp 


Ef-E^^ 


ksT 


(r/2) 


(/iw-£/i)-  +  (r/2)2 


(9) 


with  the  matrix  element 


Mfi 


niiAEjf'-Ef^’) 


where  £/, 


ih 


(•LH 


(10) 


and  £/^'  and  £/^*  denote  the 
quantized  energy  levels  for  the  initial  state  and  final  state, 
respectively,  q  is  the  permeability,  c  is  the  speed  of  light  in 
free  space,  ks  is  Boltzmann’s  constant,  T  is  the  tempera¬ 
ture,  B  is  the  angle  between  the  polarization  vector  and  the 
normal  to  the  quantum  well,  n,  is  the  refractive  index,  Ef 
is  the  Fermi  energy  which  depends  on  the  density  of  elec¬ 
trons  in  the  well,  and  F  is  the  linewidth. 

The  oscillator  strength  /  is  given  by ' ' 


I - 

(solid  line)  for  an  effective  well  width  101,27  A,  which 
gives  the  same  ground-state  energy  for  £“0  with  the  true 
well  width  of  65  A  and  the  barrier  height  A£c“245 
meV."  We  use  £f— 6.49  meV  which  corresponds  to 
about  l.6x  10”  cm  electrons  and  T-'IO  meV  from  the 
experimental  results"  and  is  assumed  to  be  independent 
for  the  variation  of  F.  One  can  easily  see  that  the  transi- 


2wo(£/'’-£/") 


f- 


_ 2  \Mfi  I  ^ 

mo(£/'’-£/"T 


(2)  2 


(II) 


S  ?oook 


In  the  zero-field  limit,  /'=  14,45  for  the  1  *  2  transition, 
which  is  independent  of  the  width  of  the  well  for  the 
infinilc-poiential-well  model.  The  experimental  result"  of 
/  for  a  well  width  of  65  A  (or  effective  /,  =  101.27  A)  is 
12.2  and  slightly  depends  on  the  well  width, 

Wc  calculated  u  for  the  first  three  states  with  field 
dependence  niinicrically  for  £  =  .100  K.  In  Fig.  2,  wc  plot 
the  absorption  coelficient  a  for  the  incident  photon  with 
polarization  perpendicular  to  the  well  (d  =  0),  taking  into 
account  the  first  three  states  as  a  function  of  the  energy  of 
the  photon  with  £=0  (d.ished  line)  and  /•'”25()  kVAm 


H(i.  2.  Comparison  of  the  intersuhhand  absorption  eocllicicnl 
n  for  an  infinite  well  width  /.  =  101 .27  \  for  the  first  three  states 
with  electron  density  l.fix  10'  /cm  ’  electrons  tor  the  zero  elec¬ 
tric  field  (dashed  line)  and  for  the  electric  field  of  2.^0  k\7cm 
(solid  line) 
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tion  1  2  is  dominant  for  both  electric  fields  F* **0  and 

250  kV/cm.  f  or  F“250  kV/cm  the  absorption  peak  is'- 
shifted  by  1 6  ineV  from  1 65  to  1 8 1  meV,  and  the  peak  am¬ 
plitude  is  increased  from  3153  to  4619  cm“'.  There  are 
two  distinct  features  for  the  case  of  the  intersubband  ab¬ 
sorption  compa'^ed  with  the  exciton  absorption.’ 

(i)  The  absorption  peak  for  intersubband  optical  absorp¬ 
tion  is  increased  in  energy  with  increasing  electric  field 
over  a  wide  range  of  the  electric  field,  because  for  increas¬ 
ing  electric  fields  the  energy  of  the  ground  state  decreases 
rapidly,  while  those  of  the  higher  subband  states  increase 
slightly  then  decrease  slowly  as  cited  in  Ref.  6.  On  the 
other  hand,  for  the  exciton  absorption,  both  the  ground 
states  of  the  electrons  and  the  holes  decrease.  Thus  the 
absorption  peak  is  decreased  in  energy  with  increasing 
electric  field. 

(ii)  The  absorption  peak  for  intersubband  optical  ab¬ 
sorption  is  increased  in  magnitude  with  increasing  electric 
field  because  the  eiectrons  are  shifted  to  the  same  side  of 
the  well  for  both  the  initial  and  the  final  states  with  in¬ 
creasing  electric  field,  and  the  energy  difference  Ei  —  E\ 
also  increases  for  the  reason  mentioned  in  (i).  As  a  result, 
the  absolute  value  of  the  overlap  integral  <V//,  for  1 2 
transition  increases.  For  the  exciton  absorption,  increasing 
electric  field  causes  further  separation  of  electrons  and 
holes  in  the  well  as  well  as  the  decrease  of  the  energy 
difference  between  the  electron  and  the  hole  ground  states, 
thus,  the  decrease  of  the  absolute  value  of  the  overlap  in¬ 
tegral.'' 

It  is  ai.so  remarkable  that  the  forbidden  transition  1  — *  3 
for  /■  =0  becomes  possible  when  F  is  nonzero  because  the 
parity  which  prohibits  the  transition  1  -*  3  no  longer  exists 
when  F  is  nonzero.  In  Fig.  3  we  plot  as  a 

function  of  F  for  the  1  2  transition,  where  is  the 

value  of  Mfi  for  the  1  2  transition  with  F="0.  One  can 

easily  see.  as  expected,  that  the  ratio  increases  slightly 
from  I  as  F  increases.  In  our  calculation,  we  assume  T  is 
constant;  however,  to  account  for  the  effect  of  the  electric 


FIG.  3.  The  normalized  overlap  integral  1 3/2i/.3/|fM  ^  where 
MkV  is  for  the  zero  electric  field,  is  plotted  vs  electric  field  F. 


field  on  the  absorption  completely,  further  analysis  of  the 
electric  field  dependence  of  the  line  broadening  is  desired. 

In  conclusion,  we  have  calculated  the  electric  field 
dependence  of  the  intersubband  absorption  within  a  con¬ 
duction  band  of  a  quantum  well.  It  is  found  that  the  ab¬ 
sorption  peak  is  shifted  in  energy  and  is  also  increased  in 
magnitude  with  increasing  electric  field.  The  forbidden 
transition  1  — -  3  when  F=0  becomes  allowable  for  the 
nonzero  electric  field. 

This  work  was  partially  supported  by  NASA  Grant  No. 
NAG  1-500  and  Air  Force  Contract  No.  F33615-84-K- 
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Nonlinear  intersubband  optical  absorption  in  a  semiconductor  quantum  well 

D.  Ahn  and  S.  L  Chuang 

Department  of  Electrical  and  Computer  Engineering,  University  of  Illinois  at  Urbana-Champaign,  Urbana, 

Illinois  6I80I 

(Received  6  February  1987;  accepted  for  publication  8  June  1987) 

The  third-order  nonlinear  intersubband  absorption  in  a  semiconductor  quantum  well  is  studied 
theoretically  using  the  density  matrix  formalism  including  intrasubband  relaxation.  It  is  shown 
that  the  peak  absorption  is  reduced  by  half  for  an  optical  intensity  1  MW/cm^  for  the  well  size 
L  =  126.5  A  with  3. Ox  lO'Vcm' electrons. 


The  optical  properties  of  the  two-dimensional  electron 
gas  in  quantum  wells  are  of  growing  interest.  The  linear  opti¬ 
cal  properties  including  the  clectrooptical  absorption  for  the 
interband  transitions  have  been  extensively  stiidicd  '  'Non¬ 
linear  optical  properties  and  potential  device  applications  of 


bulk  semiconductors  and  quantum-well  structures  have 
been  studied  recently  (or  interhand  absorption.'  Recent  ex¬ 
perimental  studies  show  that  the  optical  transitions  between 
the  subbands  in  a  conduction  band  of  quantum  well  have  a 
large  dipole  moment.'"  "  This  suggests  that  intersubband 
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optical  transitions  in  a  quantum  well  may  have  very  large 
optical  nonlinearities.  Stark  effects  on  the  intersubband  opti¬ 
cal  absorption  have  been  studied  experimentally'^  and  theo¬ 
retically.'^  Yuen'"'  has  calculated  the  intersubband  optical 
nonlinearities  due  to  carrier  confinement  neglecting  the  col¬ 
lision  broadening  of  the  subbands  and  concluded  that  the 
nonlinear  optical  properties  are  caused  by  the  electron  redis¬ 
tribution  among  the  subbands  by  the  optical  field.  But  neg¬ 
lecting  the  collision  broadening  leads  to  the  divergence  of  his 
results  near  the  subband  edge.  In  this  communication,  we 
present  theoretical  calculations  for  the  nonlinear  optical  ab¬ 
sorption  associated  with  intersubband  transitions  taking  into 
account  the  intrasubband  relaxation.  Our  analysis  is  based 
on  the  density  matrix  formulation  with  the  relaxation-time 
approximation. '' 

Let  us  consider  an  electron  with  charge  -|e|  and  effective 
mass  m*  in  an  infinite  quantum  well  with  a  width  L  in  the 
presence  of  an  optical  radiation  with  the  polarization  along 
the  direction  of  the  well  z  (TM  polarization).  The  matrix 
elements  for  the  other  polarization  in  the  intersubband  ab¬ 
sorption  is  much  weaker  or  nearly  zero.  We  choose  the 
origin  to  be  at  the  center  of  the  well.  Letp  be  the  one-electron 
density  matrix  for  this  system,  Hq  be  the  unperturbed  Hamil¬ 
tonian  for  the  infinite  quantum  well,  M  be  a  dipole  operator 
which  has  only  off-diagonal  matrix  elements,  Af,,* 
due  to  the  symmetry  of  the  system,  and  E(t)  be  the  electric 
field  strength  of  the  optical  field  with  angular  frequency  ta. 
When  the  quantum  well  is  not  symmetric,  e.g.,  a  symmetric 
well  with  an  applied  static  electric  field,  the  dipole  operator 
may  have  diagonal  elements  which  require  special  consider¬ 
ations.'*’'’  The  one-electron  density  matrix  equation  with 
intrasubband  relaxations  becomes '  ’ 

+  (1) 

at  in 


in  powers  of  the  electric  field  strength  E  are  given  as 
P(.t)=eoX(nE(t) 

=  ^'Ll[Ma„P^(t)+M^p^,{t)],  (5) 

^  K  k; 

where  V  is  the  volume,  and  the  factor  of  2  accounts  for  the 
summation  of  the  spin  variables.  The  susceptibility  x  is  relat¬ 
ed  to  the  absorption  coefficient  a  by 

a  =iy.y/(/r/e^  )Im[eQj^((y)]  ,  (6) 

where//  is  the  permeability,  is  the  real  part  of  the  permit¬ 
tivity,  and  Im  denotes  the  imaginary  part. 

We  solve  Eq.  ( 1 )  pwrturbatively  by  expanding^  in  pow¬ 
ers  of  £  as 

P  =  Ip‘’”  (7) 

n 

with  the  unperturbed  density  matrix  p‘°’  having  only  diag¬ 
onal  terms.  We  neglect  higher  harmonic  terms  which  corre¬ 
spond  to  the  successive  absorption  or  emission  of  photons 
and  consider  only  the  steady-state  responses.  Under  these 
assumptions,  the  nth-order  perturbation  term/j*"’  (r)  can  be 
written  as 

=  (8) 

Furthermore,  p  has  the  symmetric  property, 
FaftCO  =pta(>)- Using Eqs.  (1),  (3),  (4),  (7), and  (8),  we 
obtain  pU’((i>)  and  pi^’(w)  after  some  mathematical  mani¬ 
pulations.  The  results  are  as  follows: 


pU’(a>) 


2^(  VEi^-co-  iYab ) 


(9) 


where  [  ,  ]  is  the  quantum-mechanical  commutator,  ^  is  the 
Planck  constant  divided  by  Itt,  and  T  is  the  phenomenologi¬ 
cal  operator  responsible  for  the  damping  due  to  the  electron- 
phonon  scatterings,  collisions  with  electrons,  etc.  We  as¬ 
sume  that  r  is  a  diagonal  matrix  and  its  element  r„„  is  the 
inverse  of  the  relaxation  time  for  the  state  |  n ) . '  ’  For  simpli¬ 
city,  we  will  consider  only  the  two  lowest  subbands  of  the 
well  in  the  conduction  band.  We  introduce  two  shorthand 
notations  a  and  b  such  that 

|a)  =  |l,k,),|6)  =  |2,k;>  ,  (2) 

where  k,  and  k;  are  the  wave  vectors  of  the  electron  in  the  x-y 
plane,  and  1  and  2  refer  to  the  quantized  z  components  of  the 
energy  of  the  first  and  the  second  states,  respectively.  H(,  is 
diagonal  in  \a)  and  |6>  with  energies  £■„  and  ,  respective¬ 
ly.  The  monochromatic  incident  field  is 


F(/)=Re(F,/? 


(3) 


=  J  £■„<?  -  -I-  \E  ?  e"-" 

and  the  diagonal  elements  of  the  operator  T  are 

(b\r\b)  =  \/t^,  and  (air|/7)  =  l/r„  ,  (4) 

where  r„  and  t,,  are  the  relaxation  times  for  the  states  lo) 
and  \b  ),  respectively. 

The  electronic  polarization  P(t)  and  susceptibility  \{t) 


and 


piV) 


^  Fo|Fol’(r„  +  n  Wi.  I  iplV  . 
4/P(lF^ -co-in,)[iW^-wy-  +  yi,] 

1 


(10) 


where  =  {b  \M  |a>,  =  (F(,  -  F„  )/fi,  and 

Xofc  =  l/r„  -I- 1/rj )  .  For  simplicity,  we  assume  that 
Ta  =Tt,  ~  Tin  ,  the  intrasubband  relaxation  time.  Then, 
from  Eqs.  (5),  (6),  (9),  and  ( 10),  we  obtain  the  linear  and 
the  third-order  absorption  coefficients  a" ' (<u )  and  a"' (a;,/ ) , 
respectively. 


o"'(®) 


1  -t-  exp(Ff  -  F|)/^gr\ 
\  1  -I-  exp(Ff  —  E2)/k„  T ) 


X 


X 


1  -I-  exp(Ff  - 


(F,  -  E,  -  fko]'  +  (^/r,„  )’ 


and 
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where  1  is  the  optical  power  per  unit  area,  n,  is  the  refractive 
index  (for  GaAs  —  3.2,)  Eq  is  the  permittivity  of  free 
space,  c  is  the  speed  of  light  in  free  space,  m*  is  the  effective 
mass  of  electron,  T  is  the  temperature,  kg  is  the  Boltzmann 
constant,  and  In  denotes  the  natural  logarithm.  The  matrix 
element  for  a  quantum  well  has  been  evaluated  analyti¬ 
cally  in  Eqs.  (11)  and  ( 12).  In  deriving  Eqs.  (11)  and  (12), 
we  substitute /„  and  fg ,  the  Fermi  distribution  function  with 
a  given  Fermi  level  Eg  of  the  states  a  and  6,  for  the  unper¬ 
turbed  diagonal  density  matrix  elements  and  and 
use  the  identity 

l-\(f  _  f  -  I  (  *  +  t^pjEg  -  E^)/kBT\ 

\\+txx>{Eg-E:,)/kgT)' 


I  +  exp(£^.  —  E^)/kgT ) 


The  total  absorption  coefficient  a(o),f)  is  given  by 

a (<y,/)  =  a' "(<y) -I- «'’•(&>,/).  (14) 

We  have  calculated  a  with  intensity  dependence  nu¬ 
merically.  The  parameters  we  used  are  as  follows: 

m*  =  0.0665w„,  7=11  K,  L  =  126.5  A, 

=  25.38  me V,  ri„=0.14ps.  (15) 

We  used  an  effective  well  width  2,  =  126.5  A,  chosen  to  give 
the  same  ground-state  energy  of  a  finite  well  with  a  width 
Lg  =  100  A  and  the  barrier  height  F„  =  340  meV  in  the 
GaAs  —  AljGa,  As  system.  The  Fermi  level 
Ef  =  25.38  meV  implies  that  there  are  about  3.0  X  lO'Vcm’ 
electrons  in  this  well  at  T  =11  K.  Our  choice  of 
=  0. 14  ps  should  be  reasonable  since  the  experimental 


FIG  I .  The  peak  absorption  cocfncicni  a  is  plolied  vs  intensity  for  an  infi* 
nitc  well  with  a  width  L  ■-  1 26.5  A  and  the  electron  density  3  0  «:  lO'Vcm' 
for  r-  77  K. 
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Photon  energy  tm^(ev) 

F!G  2  Comparison  of  the  iniersubband  absorption  coefficient  a  for  an  infi¬ 
nite  well  with  a  width  L  =  126.5  A  and  the  electron  density  3. Ox  lO'Vcm^ 
at  T*  =  77  K  for  three  different  optical  intensities:  (i)  /  —  0,  (ii)  /  =  0.5 
MW /cm’,  and  (iii)  7=  1.0  MW/cm’. 


results  show  subpicosecond  relaxation  time.'®  In  Fig.  1,  we 
plot  the  peak  of  the  absorption  coefficient  a  (cm  ” ' )  which 
occurs  atfe  =  0. 106eV  for  various  intensities  /  (kW/cm^) 
with  light  propagating  along  the  well  with  its  polarization 
perpendicular  to  the  well.  It  can  be  seen  from  the  figure  that 
the  saturation  intensity  is  1  MW/cm^  where  the  peak  ab¬ 
sorption  coefficient  is  reduced  by  half  In  Fig.  2,  we  plot  the 
absorption  coefficient  a  versus  the  incident  photon  energy 
^(eV)  for  the  same  well  for  three  different  optical  intensi¬ 
ties:  (i)  /  =  0,  (ii)  /=  0.5  MW/cm^,  and  (iii)  /  =  1,0 MW/ 
cm*.  We  see  that  the  strong  absorption  saturation  begins  to 
occur  at  around  /=  1.0  MW/cm^  The  large  third-order 
nonlinearities  are  the  result  of  modulation  of  phase  coher¬ 
ence  by  the  strong  oscillating  fields. 

In  conclusion,  we  have  studied  the  third-order  nonlin¬ 
ear  intersubband  optical  absorption  taking  into  account  the 
intrasubband  relaxation.  It  is  shown  that  the  magnitude  of 
the  peak  absorption  coefficient  is  reduced  by  half  when  the 
intensity  is  1  MW/cm^  for  the  well  width,  L  =  126.5  A. 

This  work  was  partially  supported  by  the  NASA  Grant 
No.  NAG  1-500  and  the  Air  Force  Contract  No.  F33615-84- 
K-I557.  D.  Ahn  would  al.so  like  to  thank  General  Telephone 
and  Electronics  Laboratories,  Inc.  for  financial  support. 
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Calculation  of  Linear  and  Nonlinear  Intersubband 
Optical  Absorptions  in  a  Quantum  Well  Model  with 

an  Applied  Electric  Field 

DOYEOL  AHN  and  SHUN-LIEN  CHUANG,  member,  ieee 


AnaiyCic  forms  of  th«  linoar  and  the  third>order  nonlinear 
optical  intersubband  absorption  coefficients  are  obtained  for  general 
asymmetric  quantum  well  systems  using  the  density  matrix  formalism, 
taking  into  account  the  intrasubband  relaxation.  Based  on  this  model, 
we  calculate  the  electric  field  dependence  of  the  linear  and  the  thirds 
order  nonlinear  intersubband  optical  absorption  coefficients  of  a  semi¬ 
conductor  quantum  well.  The  energy  of  the  peak  optical  intersubband 
absorption  is  around  100  meV  (wavelength  is  12.4  )im).  Thus,  elec- 
trooptical  modulators  and  photodetectors  in  the  infrared  regime  can 
be  built  based  on  the  physical  mechanisms  discussed  here.  The  con¬ 
tributors  to  the  nonlinear  absorption  coefficient  due  to  (he  electric  held 
include  1)  the  matrix  element  variation  and  2)  the  energy  shifts.  Nu¬ 
merical  results  are  illustrated. 


I.  Introduction 

Quantum  confinement  of  carriers  in  a  semiconduc¬ 
tor  quantum  well  leads  to  the  formation  of  discrete 
energy  levels  and  the  drastic  change  of  optical  absorption 
spectra  [  1 ) .  One  of  the  most  remarkable  properties  of  these 
quasi-two-dimensional  electronic  stems  is  that  the  op¬ 
tical  transitions  between  the  size-quantized  subbands  are 
feasible.  Recently,  the  linear  intersubband  optical  absorp¬ 
tion  within  the  conduction  band  of  a  GaAs  quantum  well 
has  been  studied  experimentally  without  an  electric  field 
[2],  [3]  and  with  an  electric  field  [4],  A  very  large  dipole 
strength  and  a  narrow  band  width  were  observed.  This 
suggests  that  the  intersubband  optical  transitions  in  a 
quantum  well  may  have  very  large  optical  nonlinearities. 
The  linear  intersubband  absorption  with  an  applied  elec¬ 
tric  field  was  calculated  by  the  authors  [5]  and  nonlinear 
intersubband  absorption  without  an  electric  field  has  been 
calculated  in  several  papers  [6],  [7],  These  linear  and 
non  near  size-quantized  transitions  are  interesting  both 
from  a  fundamental  and  a  practical  point  of  view.  These 
have  the  potential  for  device  applications  in  far-infrared 
(wavelength  -  10  /xm)  laser  amplifiers  [8],  [9],  photo¬ 
detectors  [9],  [10],  and  high-speed  elecirooptical  modu¬ 
lators  Fin. 
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(0)  F  =  0 


lb)  F  >0 


Fig.  1.  Schematic  diagrams  for  energy  levels  and  wave  functions  of  uu 
inhniie  quantum  well  lor  (a)  F  =  0  and  tb)  F  >  0. 


Consider  an  infinite  quantum  well  with  a  width  L  with¬ 
out  .m  applied  electric  field  [Fig.  l(a)]  and  with  an  elec¬ 
tric  field  (Fig.  I(b)j.  Previous  calculations  [12],  [13]  show 
that  1)  the  energy  level  of  the  ground  state  £|  lowers  and 
the  first  excited  state  f.  raises  slightly  when  there  is  an 
intermediate  applied  electric  field.  2)  The  wave  functions 
are  pushed  to  one  side  of  the  quantum  well.  Detailed  nu¬ 
merical  results  for  the  energy  levels  E„  in  terms  of  the 
normalized  parameters  £„  and  £  are  shown  in  Fig.  2(a) 
where  £„  =  E„/E\°\  /{2m*L^)  and  £  = 

lelFL/£',“’.  The  amplitudes  of  the  wave  functions  are 
shown  in  Fig.  2(b).  Thus  it  is  interesting  to  see  that  I)  the 
energy  difference  £.  -  £,  corresponding  to  the  peak  en¬ 
ergy  for  optical  absorption  when  there  is  an  incident  light 
is  raised  when  the  electric  field  £is  increased,  and  2)  the 
dipole  matrix  element  |  < '/'.  1  ec  |  iZ-i  >  1  decreases  slightly, 
and  the  magnitude  of  the  absorption  coetficieni  assuming 
a  C’  -nstant  line  width  increases  with  an  increasing  electric 
fie;J  if  the  Fermi  level  is  assumed  to  be  fixed  (£f  -  £, 
thus  increases).  If  the  line  width  broadening  effects  due 
to  tunneling,  scattering  processes,  and  inhomogeneity 
with  the  electric  field  are  included,  the  peak  value  will 
drop  as  seen  by  a  recent  experiment  [4]. 

In  this  paper,  we  derive  analytical  forms  for  both  linear 
and  nonlinear  optical  absorption  coefficients  for  a  general 
asymmetric  quantum  well  system  which  may  consist  of  a 
symmetric  quantum  well  with  an  applied  uniform  electric 
field.  We  consider  a  two-level  system.  Our  analvsis  is 
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Fig.  2.  fa)  Normalized  subband  energies  f  arc  plotted  ver¬ 

sus  the  normalized  electric  field  kP  ^  lei  FL!  E'l" )  ( I0|.  f  1 )  The  exact 
s«)lui»on  —  (solid  lines).  (2)  (he  variational  method  with  Grum-Schmidt 
onhogonalization  procedure  —  (dashed  lines),  (b)  Magnitude  ol  the 
wave  functions.  L  |  (/'.,( 2)t\  are  plotted  versus  the  normalized  distance 
:/Z.  for  the  first  three  subbands.  The  dotted  lines  are  for  the  zero  olectnc 
field.  The  solid  lines  are  ihe  exact  solutions  and  the  dashed  tines  are 
those  for  the  variational  method  with  the  Gram-Schmidt  onhogonaliza- 
Mon  procedure  for  a  quantum  well  with  an  applied  electric  field  iP  ~ 
20)  110). 


bu.sed  on  the  one-particle  density  matrix  formalism  taking 
into  account  the  intrasubband  relaxation.  We  then  calcu¬ 
late  the  linear  and  the  third-order  nonlinear  intersubband 
absorptions  in  a  semiconductor  quantum  well  with  an  ap¬ 


plied  electric  field  as  a  special  case  of  the  general  for¬ 
mulation  for  asymmetric  quantum  well  svstems  assuming 
constant  intrasubband  relaxation  time.  Our  calculation  is 
based  on  the  infinite  quantum  well  model  with  an  effective 
well  width  which  gives  the  same  ground  state  energy  of  a 
given  finite  quantum  well.  When  the  barrier  height  of  the 
quantum  well  is  finite,  it  is  weil  known  that  the  wave 
functions  of  the  electrons  leak  into  the  barrier  region.  Thus 
the  effective  well  width  L.,f  is  usually  larger  than  the  true 
well  width.  The  ground  state  properties  of  realistic  quan¬ 
tum  well  systems  have  been  discussed  extensively  in  the 
literature  [17].  Our  previous  calculations  [18]  show  that 
the  ground  state  energy  of  an  infinite  quantum  well  with 
an  effective  well  width  L  =  126.5  .A  agrees  very  well 
with  that  of  a  finite  quantum  well  with  width  Lj  =  100  A 
and  the  barrier  height  k'o  =  340  meV.  up  to  the  electric 
field  F  =  200  kV /cm.  The  absorption  peak  for  the  states 
1  —  2  transition  is  shifted  from  106  meV  when  F  =  0  to 
112  meV  when  F  =  100  kV /cm.  These  energies  corre¬ 
spond  to  optical  wavelengths  of  11.7  ixm  (at  106  meV ) 
and  11.1  ^m  (at  112  meV).  We  assumed  that  the  intra¬ 
subband  relaxation  time  t,„  is  constant.  In  the  real  quan¬ 
tum  well,  relaxation  times  of  the  electron-phonon  scat¬ 
terings  and  electron-electron  scatterings  may  strongly 
depend  on  the  applied  electric  field.  When  the  quantum 
well  has  a  finite  barrier,  the  tunneling  effect  tends  to  in¬ 
troduce  another  time  constant  for  broadening  of  the  line 
width  (i.e..  the  energy  levels  become  complex  [18]  £„  - 
(■  r/2).  Both  the  scattering  processes  and  the  tunneling 
effects  together  with  inhomogeneity  may  enhance  the  line 
width  factor. 

Nonlinear  optical  properties  in  bulk  semiconductors  and 
quantum  well  structures  have  been  studied  recently  [19]- 
[23].  However,  most  of  these  works  consider  the  inter¬ 
band  process;  the  electric  field  effect  is  not  considered.  In 
this  paper,  we  consider  the  intersubband  process  with  an 
electric  field  applied  to  the  quantum  well  stuctures. 

The  one-clectron  density  matrix  equations  are  devel¬ 
oped  and  solved  for  the  sinusoidal  steady  state  taking  into 
account  the  intrasubband  relaxation  for  asymmetric  quan¬ 
tum  well  systems  in  Section  II-A.  Analytical  forms  of  the 
linear  and  the  third-order  nonlinear  absorption  coefficients 
of  asymmetric  quantum  well  systems  are  derived  in  Sec¬ 
tion  Il-B,  Section  Ill-A  gives  expressions  for  the  dipole 
moments  of  the  infinite  quantum  well  with  an  applied 
electric  field.  Results  obtained  in  Section  II-B  are  used  (in 
Section  III-B)  to  calculate  the  linear  and  the  third-order 
nonlinear  intersubband  coefficients  of  a  semiconductor 
quantum  well  with  an  applied  electric  field.  Finally,  we 
present  numerical  results  and  discussions  in  Section  IV. 

II.  iVlATHrviATir '.L  .'^ORMLLATIO.N  Ol  Ll.NhAR  AND 

Nonlinear  Optical  Absorptions  for  Assaimetric 
Quantlm-Well  Systems 

In  this  section,  we  present  the  theoretical  formulation 
for  both  the  linear  and  nonlinear  optical  absorption  coef¬ 
ficients  for  asymmetric  quantum  well  svstems  The  cal- 
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culations,  though  mathematical,  are  quite  elementary.  The 
application  of  this  formulation  to  the  intersubband  optical 
absomtion  in  a  quantum  well  with  an  applied  electric  field 
is  then  considered  in  Section  III. 

A.  Density  Matrix  Equations 
We  consider  an  asymmetric  quantum  well  in  the  pres¬ 
ence  of  an  optical  radiation  of  angular  frequency  u  with 
the  polarization  along  the  direction  of  the  well  c.  Cases 
for  the  arbitrary  polarization  will  introduce  an  additional 
factor  cos  d  in  the  matrix  element  where  d  is  the  angle 
between  the  z-axis  and  the  direction  of  polarization  [5]  for 
intersubband  transitions.  Let  p  be  the  one  electron  density 
matrix  for  this  system,  Wq  be  the  unperturbed  Hamiltonian 
for  this  system  with  an  asymmetric  potential  energy,  M 
be  a  dipole  operator,  and  £(r)  be  the  electric  field  strength 
of  the  optical  radiation  with  frequency  cj.  Since  we  are 
considering  an  asynimetric  quantum  well  system,  the  di¬ 
pole  operator  may  have  nonvanishing  diagonal  elements, 
that  is  different  from  the  symmetric  system  which  has  only 
off-diagonal  elements  [14].  The  one-electron  density  ma¬ 
trix  equation  with  intraband  relaxation  becomes  [15] 

-i[r(p-p'"')  +  (p-p'°')r]  (1) 

where  [.]  is  the  quantum  mechanical  commutator,  h  is  the 
Planck’s  constant  divided  by  lit,  p'°'  is  the  unperturbed 
density  matrix,  and  F  is  the  phenomenological  operator 
responsible  for  the  dampling  due  to  the  electron-phonon 
interaction,  collisions  among  electrons,  etc.  We  assume 
that  r  is  a  diagonal  matrix  and  its  element  is  the  in¬ 
verse  of  the  relaxation  time  for  the  state  |ni>  [15].  For 
simplicity,  we  shall  confine  our  attention  to  two-level  sys¬ 
tems  only.  We  introduced  two  short-hand  notations  |a> 
and  I  ) ,  such  that 


|a>  =  |l.  k,) 


1^)  =  |2.  k\) 


£„  =  £, 


2/nr 


c*  =  £2  k',  (3b) 

where  m* ,  m*  are  the  effective  masses  for  the  initial  and 
the  final  slates,  respectively,  for  the  intersubband  transi¬ 
tion,  We  assume  m*  =  m*  =  m* ,  the  effective  mass  of 
electrons  near  the  conduction  band  edge.  The  effect  of 
nonparabolicity  of  the  band  structure  is  neglected  here.  It 


can  be  included  by  using  different  effective  masses  m*  and 
m*.  or  more  correct  expressions  for  f„  and  E^.  The  mono¬ 
chromatic  incident  fielH  is  defined  as 

E(t)  =  Re  (E'’?-""') 

=  i£'’e'‘‘'  (4) 

=  £e’'“"  4-  E*e‘-‘ 


and  the  diagonal  elements  of  the  operator  F  are 

<f>|F|£>  -  1/r*,  and  <alF|a)  =  l/r,  (5) 

where  and  T/,  are  the  relaxation  times  for  the  states  |  a  > 
and  \  b),  respectively. 

We  solve  (1)  perturbatively  by  expanding  p  in  powers 
of  £  as 


with  the  unperturbed  density  matrix  p'“'  assumed  to  have 
only  diagonal  terms  [15].  Let  pj^’  =  <  a  |  p""  |  a  >.  pj^' 
=  <a|p''-'|£>.  pZ'  =  (^Ip'^’k),  and  p\V  = 
(.b\  p'"'\by  Then,  p  has  the  symmetry  property  that 
Pab(l)  =  pti(i)-  We  obtain  four  equations  for  each  ele¬ 
ment  of  the  density  matrix  from  (1)  and  (6). 

Using  (l)-(6),  we  obtain  the  (n  +  1  )th  perturbation 
terms  for  /t  >  0 


3/  ■  lift 


pZ^'^  -]-AEb-  EJ  -  y^ipZ^'^ 


j^{pZ’  -  pZ')MboE{t) 


-  -  -Vf„J£(r)p'"„’ 


^  ^  7bii  ( 


-  pZ')Kt,E{t) 


where  k,  and  k'  are  the  wave  vectors  of  the  electron  in 
the  xy  plane,  and  1  and  2  denote  the  size-quantized  z- 
components  of  the  momenta  for  the  initial  and  the  final 
states,  respectively.  Hq  is  diagonal  in  |a>  and  |fi>  with 
energies  £^  and  £(,,  respectively,  which  are  given  as 


-  M^)EU)p':Z 


^hhPhb  .  *  i  P  ab 


-  m^pZ' 


-M^pZ')E{t)  flO) 

where  Ma„  =  (  a  |  iW  |  p  > .  .W,,*  =  (  a\M  \by  M^,  = 
<  I  W  I  a> ,  M^^  =  <  £  I  -Vf  ]  =  1  /  T^.  =  1  /t(,. 

and  =  7s<,  =  :(  1  /t„  +  I  /t*).  Equations  (7)-(10i  are 
readily  solved  by  expanding  the  density  matnx  elements 
as  sums  of  terms  proportional  to  exp  (  ±iu:i)  and  equating 
terms  on  both  sides  having  the  same  lime  dependence.  Wc 
will  neglect  higher-harmonic  terms  which  correspond  to 
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the  successive  absorption  or  emission  of  photons  and  we 
consider  only  the  steady-state  response.  Under  these  as¬ 
sumptions.  the  nth  order  perturbation  term  p'"\t)  can  be 
wiiiicn  as 

p^"'(n  -  ^""(01)^-'-'  ^  (11) 

for  odd  n.  When  n  is  even,  only  dc  terms  are  dominant. 
Using  (4),  (7)-(ll),  and  neglecting  the  olT-resonancc 
terms,  we  obtain  and  after  some  mathe¬ 

matical  manipulations  (Appendi.x): 

,  EMHoip'l'  ~  P'^1) 

(“)  =  ITTir— 7— — T  (12) 


Pii'(w)  = 


A  (  -  to  -  /7.J.) 

-E\E\-M^ip‘°'  -  pJD 
h\iV,^  -  -  iy^) 

,  p((l/7„,.)  + 

L  ( -  to)'  -f-  7.-,* 

_ -  A/„J~ 

(  -  iyo>,)(  -  ta  -  h.,,). 


2iE  E  'IM. 


^■*(<0  +  hhh) 


(A/,„  -  -  pjg*) 

(  ^tsi  *“  t7,,/, ){  “  ca  —  iy 


P2’(ta)  =  - 


2iE\E\-\M,, 

A^UI  -P  /7,„ 


^ (0)  _  f 

Phh  -  Jb 


1  +  exp  [(f^  -  E,)/Kf,T\ 


where  Ef  is  the  Fermi  level  of  the  system,  T  is  the  tem¬ 
perature,  and  Kg  is  the  Boltzmann  constant. 

D.  Linear  u/iu  .  .jniinear  Opiuat  Absorptions  of  the 
Asymmetric  Quantum  Systems 

The  electronic  polarization  P{t)  and  susceptibility  \  (t) 
caused  by  the  optical  field  E(  t)  can  be  expressed  through 
the  dipole  operator  \l  and  the  density  matrix  as 

P{t)  =  e„x(ta)£:c-'“'  +  e,x(-ta)EV“' 

=  -pTr(pW)  (15) 

where  Kis  the  volume  of  the  system,  Eq  is  the  permi"ivity 
of  free  space  and  Tr  denotes  the  trace  or  summation  over 
the  diagonal  elements  of  the  matrix  pM.  The  susceptibil¬ 
ity  X  is  related  to  the  absorption  coefficient  cir(<a)  by 

^(o))  =  O)  I— Im  (eox(‘a))  (16) 

V  £/i 

where  p  is  the  permeability  of  the  system,  e«  is  the  real 
part  of  the  permittivity,  and  xfta)  is  the  Fourier  compo¬ 
nent  of  x(t)  with  e“'“'  dependence.  From  (12),  (13),  (15), 
and  (16),  we  obtain  the  linear  and  the  third-order  absorp¬ 
tion  coefficients  a"’(w)  and  o(‘^’(w,  / ),  respectively, 

a"*(w)  =  ‘^J~  p  S  S  I  ^bo\ 

(pL!!*  -  pZ^)fiy.., 

(E,  -  -  hwf  +  {hyj 


a'^Vw,  /)  =  -ca,  -^^Im  SZ 

\  fR  r  i,  L 


.  r  (M,,  -  m„j(p:°' -  pj,?,’) 

X^ba  ^ytib')i^bfi  vJ  iy^ig) 

(13c) 

where  =  (£*  -  E,, ) /A,  and  Im  denotes  the  imaginary 
part.  However,  the  contributions  of  p2'(w)  and  pU’lw) 
to  a'^'(ta,  / )  later  turn  out  to  be  negligible.  The  diagonal 
elements  of  the  unperturbed  density  matrix  p“”  are  the 
equilibrium  electron  state  occupations  determined  by  the 
Fermi  level.  So  we  can  represent  pl,^’  and  Pw’  by  the 
Fermi-Dirac  distribution  functions /„  and /(,,  respectively: 


- ^ - :  {14a) 

1  +  exp  [(E,  -  Er)/KgT] 


\E\'\Mba\  {piV  -  p!,^) 

^'(l^h,,  -  W  -  i7„a) 

^y lib  (  y isu  4-  7l>fc  )  I  ^ab  I 
7u<i  7/»fr  {  (  ^ba  W  )  b  7  (jt»  I 

(  ^ha  -  hahX  -  W  -  (7„i) 


=  -  ^^m  f.  Z  Z 


\Mba\-{p:V  -  P'^IX 

-  u,  -  ,7us) 

^  Zee ) !  I 

yuayhb{(^bc  -  <^)'  +  7^} 

'  \M,.b  -  mJ' 

(  ^'h,,  -  ‘y..h){  -  w  -  ty^) 


k k  r.jt  r. 
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where  I  is  the  optical  power  per  unit  area,  n,  is  the  re¬ 
fractive  index  (in  =  «;«(,),  c  is  the  speed  of  light  in  free 
space,  ana  we  have  neglected  ti\e  off-resona.ict  ^.-nns  in 
deriving  (18).  A  factor  of  2  has  been  added  in  (17)  to  (18) 
to  account  for  the  summation  over  two  spins  and  the  se¬ 
lection  rule  for  the  spin  orthogonality  between  the  initial 
and  the  final  states. 

III.  Line..\r  a.nd  Nonline.ar  Intersubbanp  Optical 

Absorptions  of  a  Semiconductor  Quantu.m  Well 

WITH  AN  .Applied  Electric  Field 

A.  Dipole  Moment 

Here  we  consider  an  infinite  quantum  well  with  a  width 
L  in  the  presence  of  a  constant  electric  field  F  along  the 
positive  direction  of  the  well  :.  We  choose  the  origin  to 
be  at  the  center  of  the  well.  In  the  independent  electron 
approximation,  the  wave  functions  for  the  initial  state 
and  the  final  s'.nte  in  the  conduction  band  can  be  written 
as  [16] 

=  <  r  |a> 

=  A  ■'■'-(;,(?)  exp  (<•;,  •  ?,)  4>f(z), 

\z\  <  L/2  (19a) 

'PA  r)  =  (rib) 

= /I  ■''■(/(.(  r)  exp  (i?,  •  ?,)<<);(:). 

\z\  <  L/2  (19b) 

where  A  is  the  cross  section  area  of  the  well,  r,  is  the 
position  vector  in  the  .rv  plane,  and  and  are  the  cell- 
periodic  functions  near  the  conduction  band  extremum.  In 
the  effective  mass  approximation  the  envelope  functions 
<6 1  and  0;  satisfy  the  following  Schrddinger  equation  [5], 
[17],  [18] 

!z|  <  L/2  (20) 

and  are  given  by  the  linear  combination  of  two  indepen¬ 

dent  Airy  functions  Ai(r))  and  fli(T)),  where  ij  is  defined 
by 


In  (20)  and  (21).  m*  and  e  denote  the  effective  mass  and 
charge  of  the  electron,  respectively.  For  intersubband 
transitions,  only  the  optical  fields  propagating  along  the 
well  with  polarization  along  the  direction  of  the  well  c 
have  nonvanishing  matrix  elements  [5|.  The  dipole  matnx 
eiements  can  be  approximated  by  [5],  [24] 

=  <^’i  kU!"  ) 

s  f  4>f(2)|e!2<fi,(Od2  5;,;.,  (22) 

*2  -  L  2 


M^  =  <h|!e|;|h> 


=  ] 

J  -L/: 


)|ei:  Oi{2)  c/z  (24) 


If  the  polarization  of  the  electric  field  makes  an  angle  0 
with  the  z  axis,  an  additional  constant  factor  cos  d  shoul  t 
appear  in  (22)-(24). 

B.  Linear  and  Nonlinear  Intersubhand  Absorptions  wiiL 
an  Applied  Electric  Field 

We  consider  the  intersubband  optical  transition  of  a 
GaAs-AljGa,  _,.As  quantum  well.  For  siruplirity .  we  ti'-- 
sume  that  =  Tj  =  T,n,  the  intrasubband  relaxation  time 
From  ( 14),  (17),  and  (22),  we  obtain  the  linear  absorption 
coefficient  ct' ' \l>). 


a "  ’( 0) )  =  oj  —  I  M-i  r  — 

■  2  (/a  -fh)(fi/r„) 

»•,  (£:  -  E,  -  +  (h/rj 

where  M-^  is  given  by 

p*t./2 

/Vf:i  =  gi.;*(;)|e|c  ci|(:)  Jc.  (26) 

J  ~L/2 

We  also  assume  that  r,„  is  a  constant.  The  two-dimen¬ 
sional  integration  in  (25)  can  be  done  easily.  We  obtain 

.  ,  ri  +  exp[(E^-  £,)/Ar,7-)) 

0  +  e.xp((£f  -  E.)/KsT\y 


a‘"(uj)  =  w |,V/:i| 


1  +  exp  [(£,  -  EA/K,T 


+  exp  [(£^  -  E,)/KsT\ 


( £2  -  £,  -  ho))'  -t-  [h/T^„y 

The  above  result  for  the  linear  absorption  coefficient  has 
al-so  been  derived  u.s)ng  the  Fermi  Golden  Rule  |5|.  It  can 
easily  be  seen  that  the  peak  absorption  coefficient  occurs 
at  hw  =  El  E,  from  (28).  Since  £,  -  £,  increases  with 
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an  applied  electric  field,  the  peak  absorption  energy  is  also 
shifted  hiaher.  The  maanitiiH  ■  of  th-*  peak  onf-c-ji  absorp¬ 
tion  coefficient  depends  on  1)  the  dipole  matrix  element 
j  which  decreases  slightly  with  the  electric  field, 

and  21  the  line  width  factor  which  is  assumed  to  be 

constant  here  In  fact,  the  line  width  may  be  broadened 
by  the  electric  field,  as  will  be  discussed  later. 

Now  we  calcul.ite  the  third-order  nonlinear  intersub¬ 
band  absorption  coeffici-ent  /  )  where  /  is  the  op¬ 

tical  intensity.  Using  ( IS).  (22)-(24).  and  f27),  we  obtain 
I) 


where  is  the  free-electron  mass.  We  use  an  effective 
well  width  L  =  126.5  A  chosen  to  give  the  same  ground 
state  energy  of  a  finite  well  w  ith  a  width  L,  -  100  A  and 
the  barrier  height  U,,  =  J40  me\’  n.  the  GaAs-AlG.iAs 
sysum.  The  Fermi  level  E^  =  25.3.s  meV  implies  that 
there  are  about  3.0  x  10"'/cm'  electrons  In  the  well  at  T 
=  77  K  when  F  =  0.  Our  choice  of  t,„  =  0  14  ps  is  from 
the  experimental  results  [2|.  Here.  t„  is  assumed  to  be  a 
constant.  More  work  needs  to  be  done  on  the  variation  of 


I  +  exp[(£,  - 
I  -I-  exp  [(£f  -  E:)/KijT\  j 


X 

X 


_ I'^ei  r(^/rin) _ 

(£,  -  -  huif  +  (h/Tjf 

I' 


1 

J 


where  .V/n  and  are  defined,  respectively,  by 


Thus,  the  effects  of  the  applied  static  electric  field  £  on 
the  nonlinear  optical  absorption,  similar  to  the  linear  ab¬ 
sorption  coefficient,  include  1)  the  matrix  element  |iV/,,  | 
decreases  slightly  with  £.  2)  the  peak  absorption  at  £.  - 
£,  increases  with  £  roughly,  and  3)  the  diagonal  terms 
M2:  4nd  /V/||  are  not  zero  when  F  *  0.  We  now  return  to 
the  total  absorption  coefficient  a(w,  /  1  of  the  intersub¬ 
band  transition.  From  (28)  and  (29),  the  total  absorption 
coefficient  a(ui,  I )  is  given  by 

cx(w. /)  =  a"’(aj)  -E  /).  (32) 

Note  that  when  F  =  0.  the  diagonal  elements  M,,  and  .V/,, 
become  zero  due  to  the  symmetry  property  and  a(w.  /) 
becomes  that  of  the  symmetric  quantum  well  (7|.  The  nu¬ 
merical  results  of  a  for  several  cases  will  be  presented  in 
the  next  section. 

IV.  Numerical  Results  and  Discl.s.sions 
In  this  section,  we  calculate  the  absorption  coefficient 
a(u.  /)  obtained  in  (36)  numerically  for  various  inten- 
.sities  /  and  applied  electric  fields  £.  The  parameters  are 
as  follows: 

m*  =  0.0665  mo.  £  =  77  K,  f,  =  126,5  A, 

Ef  =  25.38  meV,  n.  =  3.2.  and  r,„  =  0.14  ps  (33) 


the  relaxation  time  due  to  the  field  effect.  The  electric 
field  may  decrease  the  relaxation  time  r,„  due  to  the  elec¬ 
tron-phonon  scatterings,  the  electron-electron  scatter¬ 
ings.  the  impurity  scatterings,  and  the  inhomogeneity  ef¬ 
fects  because  of  the  electron  wave  function  change.  When 
the  quantum  well  has  a  finite  barncr  tne  tunneling  effect 
tends  to  introduce  another  time  constant  for  broadening  of 
the  line  width  (i.e.,  the  energy  levels  become  complex 
(18]  £„  -  iT/2).  Both  the  electron-phonon  scattering 
and  the  tunneling  effects  enhance  the  line  width  factor. 

In  Fig.  3.  we  plot  the  linear  absorption  coefficient  a(  uj. 
/  =  0)  =  a'"(ai)  with  £  =  0  (dotted  line)  and  F  =  100 
kV/cm  (solid  line).  We  can  see  that  the  peak  of  the  ab¬ 
sorption  coefficient  is  shifted  upward  both  in  energy  and 
in  magnitude.  Upward  shifts  of  energy  can  be  explained 
by  the  quantum  confined  Stark  effects  (5],  112].  [171  and 
were  observed  by  Harwit  and  Harris  [4|.  The  increase  of 
the  magnitude  of  the  peak  absorption  is  mainly  because 
we  assume  (hat  the  Fermi  level  Ey  is  fixed  with  an  applied 
electric  field  and  panially  due  to  the  fact  that  the  overlap 
integral  decreases  only  slightly  as  electrons  in  states  1  and 
2  are  shifted  to  the  same  side  of  the  well  by  the  applied 
electric  field.  If  the  broadening  effects  due  to  electron- 
phonon  scatterings  and  tunnelings  are  included,  the  solid 
curve  should  be  broader  in  line  width  and  the  peak  value 
will  drop  in  agreement  with  the  experimental  results  in 
|4j.  In  Fig.  4.  we  plot  the  absorption  coefficient  a(w,  /) 
in  the  absence  of  the  electric  field  (  £  =  0)  for  three  dif¬ 
ferent  optical  intensities  (7|:  (1)  /  =  0,  (iil  /  =  0.5 
,VIW/cm^  and  (iii)  /  =  1.0  MW/cnr'.  We  sec  that  the 
strong  absorption  saturation  occurs  at  /  =  1  0  MW  cnr. 
In  Fig.  5,  we  compare  the  absorption  coefficient  oc  for 
three  different  cases:  (1)  /  =  0  and  £  =  0  doited  linei. 
til)  /  =  0  and  £  =  l()0  kV  cm  (dashed  linei.  and  (ml  / 
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Fig.  3.  Compari.son  of  the  linear  intersubband  absorption  coethcicni  ctt/ 
=  0 )  for  an  infinite  well  of  width  L  -  1 2b.  5  A  with  3.0  x  10"‘/cm‘ 
electrons  for  the  zero  eler'tric  field  (dotted  line)  and  for  the  electric  Held 
of  (00  kV/cm  (solid  line)  for  7"  =  77  K. 


Fig.  5.  Intcrsubhand  absorption  coefHcieni  a  lor  an  mnnitc  well  of  a  wui"i 
L  =  126  5  .A  with  3  0  X  !()'''/cm'  electrt'ns  at  T  =  77  K  is  plotted  to: 
various  electric  fields  F  and  optical  intensities  ( /  i.  (i)  F  =  0.  /  =  C 
(dolled  line),  tit)  F  *  100  kV/cm.  /  =  0  idashcd  line;,  and  (iii)  F  = 
lOOkV  '  in.  /  -  I  t)  MW  cm*  (solid  line* 


ELECTRIC  FiElO  (xvrem) 


Fig  6-  -  iW,,  )  /.V/j,  I*  IS  plotted  versus  the  electric  licid  F. 


Fig.  4.  Intcrsubhjnd  absorption  coeHicient  or  for  an  infinite  well  of  a  width 
L  *  126.5  A  with  3.0  X  lO'^/cm*  electrons  at  f  »  77  K  with  zero 
electric  held  is  plotted  for  three  different  optical  intensities:  (i)  /  =  0. 
Ill)  /  =  0  5  MW/cm*.  and  nil)  /  »  I  0  MW/cm^  (From  151.) 


=  1.0  MW/cm‘  and  F  =  100  kV/cm  (solid  line).  We 
can  see  that  the  absorption  is  reduced  by  half  at  /  =  1.0 
MW/cm".  Again,  the  etfect  of  the  applied  electric  field 
is  to  shift  the  total  absorption  peak  to  a  higher  energy.  In 
Fig.  6,  we  plot  l(W.n  -  /W, ,  )/,V/,|  |-  for  various  electric 
fields  F.  Note  that  when  F  =  200  kV/cm.  !(.V/i.  — 
Wii)/''^:ir  becomes  larger  than  1  for  Z.  =  126.5  A. 
However,  the  net  contribution  of  the  diagonal  terms  to  the 
third-order  nonlinear  coefficient  a'^'  in  the  last  espre.s- 
sion  of  (29)  is  less  than  one  percent  when  F  =  100 
kV/cm.  This  can  be  seen  by  evaluating  the  last  term  with 
I  Mil  ~  W||  !"  in  (29)  at  £i  -  £'i  =  hu,  and  noting  that 
hui  »  h / T,„  (by  approximately  a  factor  of  10).  .So  these 


terms  can  be  neglected  for  the  moderate  range  of  electric 
field  f . 

In  Fig.  7,  we  plot  the  calculated  linear  absorption  coef¬ 
ficient  with  phenomenological  line  width  broadening  with 
F  =  0  (dashed  line)  and  £  =  36  kV/cm  (solid  line)  for 
an  infinite  quantum  well  v-ith  an  effective  width  151.2 
(or  a  finite  quantum  well  with  L,  =  120  A  and  a  barrier 
height  =  390  meV)  and  5.83  x  10'  /cm'  electrons  in 
the  well  at  T  =  90  K  to  compare  with  ihe  recent  experi¬ 
mental  result  by  Hai\sii  and  Harris  [-11  They  showed  that 
the  line  width  is  increased  by  30  percent  when  £  =  36 
kV/cm  over  that  of  the  zero  Held  in  Fig  2  of  (4).  Our 
model  IS  in  good  agreement  with  the  measured  .'^tark 
shifts.  Slight  deviation  of  intersubb.and  ibsorpnon  peak 
energy  at  zero  held  (15  meV)  is  within  the  experimental 
error  as  discussed  in  [4| 
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Fig.  7.  The  linear  intersubband  absorption  coefficicni  a{ /*  0)  for  an  m* 
finite  v.ell  with  an  effective  width  L  -  151.2  .A  wuh  5.83  x 
electrons  at  T  =  90  K  with  a  field  dependent  line  width  is  plotted  for  the 
zero  electric  field  case  (dashed  tine)  and  for  an  applied  electnc  field  of 
36  kV/cm  (solid  line)  These  (heoreiical  results  are  in  good  agreement 
with  the  experimental  results  in  |4). 


(A4) 

If  we  substitute  (A2)  and  (A3)  into  (Al),  we  obtain 
p2’(u)) 


V,  Conclusions 

An  analytic  and  useful  form  of  the  linear  and  the  third- 
'^'-der  nonlinear  absorption  coefficients  is  obtained  for 
general  asymmetric  two-level  quantum  well  systems.  We 
have  calculated  the  electric  field  dependence  of  the  linear 
and  nonlinear  intersubband  optical  absorptions  using  these 
results.  We  found  that  absorption  is  reduced  by  half  when 
/  =  1 .0  MW / cm^  for  a  quantum  well  with  a  width  L  = 
126.5  A.  The  contributions  to  the  nonlinear  absorption 
coefficient  due  to  the  electric  field  F  include  1)  the  matrix 
element  variations,  and  2)  the  energy  shifts  (quantum 
confined  Stark  effects). 


-  pZ') 

h\W^- 


L  -  w)'  +  y'ub 

_ (M,,/,  -  M^S' _ "j 

(Ka  -  -  0)  -  iyjb)}' 

We  now  derive  pl^'fui)  and  p2'(w).  From  (9)  and  (11), 
we  find 

-'wpU’(w)  =  -TssPm' 


Appendix 

Derivation  of  pii'(w).  Pm'(w).  and  p,'i'(w) 


-  ^  (^N/Pli’(0)  - 


(0)]E 


(A6) 


From  (4),  (7),  and  (11)  we  obtain 


where  we  have  neglected  the  off-resonance  terms.  If  we 
substitute  (A3)  and  (A4)  into  (A6),  we  obtain 


E 

h{W^  -  iP  -  iy^) 


•  [(pL:’(o)  -  p;,;’(o)).v/^  (At) 


2(£|£|'|iVU|' 
+  hub) 


•  \m 


A^ba  - 


-  -  p'bV) 

hub ){  W'n,  -  W  -  I7„^ )  ' 


+  (M,,  -  Muu)p'b:'{0)] 


where  p;^'(0),  pj,;'(0).  pZAO)  are  the  dc  terms  and  we 
have  neglected  higher-harmonic  terms.  Using  (7)  to  (12) 
and  neglecting  the  off-resonance  terms  such  as  /o,'J,'(w) 
and  p^;'(  ~w).  we  obtain 


(A7) 

Likewise  p.'i'l  u  )  can  be  derived  easily  following  the  same 
procedure.  Note  that  the  symmetry  propeny  p^^h)  = 
pZ(0  and  (15)  also  leads  to  /j|,^'(  o  )  =  [  p  Z'^  -  u;  )  ]*  and 
pl^'fw)  =  [  -w  ))*,  etc.  These  two  terms,  p],^'i  u  ) 
and  p'Z'i'^  ).  turn  out  to  be  negligible  when  evaluating  the 
third  order  absorption  coefficient  a’  "( w.  /  )  in  i  IS). 


52 


.y  .jfvy  i^y.  /-  if-  -r  v.  v.  ..i .  .  >f  >  y.  .■r.  v.. 


r-W  ' '■ 


:204 


IEEE  journal  of  QUANTUM  FLF.CTROMCS.  VOL  QF.  SO  l?.  DECEMBF'  :V‘s7 


References 

R.  Dingle.  W  Wiegman.  and  C  H  Hcnr>'.  ‘'Quantum  stales  of  con¬ 
fined  camcr'i  in  ^er\'  ihin  Ai,  Ga,  .  .  A\  -  GaAs  -  Al.  Ga,  ,  As  hei- 
i">is!ruciures. '■  Pin;  Rt‘\,  Lt^n  .  '■ol  .''T,  pp  S27-830.  |074 
L.  C  West  and  S.  J .  Ec'ash .  ‘ ' First  oosers  jtion  *'l  an  cxirenicly  large- 
dipoic  infrared  iran-ition  within  the  conduction  band  ot  a  GaA^  quan¬ 
tum  well. "  z-f/jn/.  phwi.  Lett.,  so),  db.  pp  1156-I15K.  Id83. 

B  F.  Levine,  R.  J  Mihk.  J  Walker.  K  K.  Choi.  C.  G  Bethea,  D. 
A.  Kleinman.  an..  J  M  Vandenbere.  “Strong  8  2  nm  infrared  inter- 
subb.ind  absorption  in  doped  CaAi  AAlAs  quantum  well  waveguide.** 
Afipi.  Phys.  Lett  .  vol  50.  pp.  273-275.  108? 

A.  Harw  itt  and  J  S.  Hans.  Jr. .  “Obser\'3tion  of  Stark  shifts  in  quan¬ 
tum  well  mtersubband  tran.sJtJor,s.  “  Appl.  Ph\s.  L*rtt..  vol.  50.  pp. 
685-687.  1^87. 

D  Ahn  and  S.  L  Chuang.  “Iniersubband  optical  absorption  in  a 
quantum  well  with  an  applied  electric  held.”  Ph\s.  B.  vo!  35. 
pp  4i49-i!5i.  1^87. 

S  Y  Yuen.  “Fast  relaxing  absorption  nonlinear  refraction  in  super- 
Ijtticcs. ' ' /J;?/'/  Phys  ten.,  vol,  43,  pp  813-815.  1983. 

D  .Ahn  and  S.  L.  Chuang,  “Nonlinear  ntersubband  optical  absorp¬ 
tion  in  a  semiconductor  quantum  well,**  J.  ^ppl.  Phys.,  Sept.  15. 
198'’. 

R.  F  Kazannov  and  R.  A.  Sun.s.  “Possibility  of  the  amplification  .>f 
eicctromapnetic  w-ave.s  in  a  semiconductor  with  a  superlatiicc.”  Sov 
Phvs.-Semiconc/. ,  vol.  5,  pp  707-709.  (971. 

F  Capasso.  K.  Mohammed,  and  A  Y  Cho.  “Resonant  tunneling 
through  double  barriers,  perpendicular  quantum  transport  phenomena 
.n  >uperiat(icc.  and  their  device  applications. '  /£££  J.  Quanaun 
E/r-rrort  .  vol  OE-22.  pp.  1853-1869.  1986 
B  F  Levine.  K  K  Choi.  C.  G.  Bethea.  J  Walker,  and  R.  J  Malik, 
■‘.New  iO  ;im  infrared  detector  using  mtersubband  absorption  m  res¬ 
onant  tunneling  GaAlA.s  superlatticcs."  Appi  Phvs  Lett.,  vol.  50. 
pp  1092- !094'  1987 

D.  Ahn  and  S.  L.  Chuang.  “Eleciro-optical  modulator  using  inter* 
subband  absorption  m  a  quantum  well  structure."  in  Proc  .Vk'  Int 
Quantum  Electron  Conf.  (IQEC'P?).  Baltimore.  MD.  1937, 

D  Ahn  and  S  L  Chuang.  “Variational  calculations  of  subbands  in 
a  quantum  well  with  uniform  electric  held;  Gram-Schmidt  onhogo* 
nalizjtion  approach."  Appl.  Phvr  Lett.,  -.ol.  49.  pp.  1450-1452. 
1986.  Note;  the  terms  t  in)  in  (9)  should  read  + 

m  I, 

M  MaisuuM  and  T  Kamizato.  “Subbands  and  exciiions  in  a  quan¬ 
tum  well  in  an  electric  field,"  Ph\.^  Rev.  B.  vol  33.  pp,  8385-8389. 
1986 

h  V  Demidov  and  Yu.  A,  Romanov.  “Noniincar  susceptibility  of 
nonsymmciric  quantum  systems."  fcv.  k'yij/i.  Uchebn  Zxived..  Ra- 
Utojiz..  vol.  28.  pp  43-50.  1985. 

N  Bluembereen.  .V<  ilineur  Optics  New  York.  Benjamin.  1965. 
ch  2. 

C  Kitlel.  Quantum  Theor,'  Sulids  New  York;  Wilev.  1963.  pp. 
286-290 


D.  A-  B.  Miller.  D  S.  Chcmla  T  C  Damcn.  A  C.  Gossaru,  W 
Wiegman.  T.  H.  Wood,  and  C  V.  Burrus.  “F.lecinc  6c!d  di.p«*niiencc 
of  optical  absorption  near  the  band  gap  nt  quantun.-woil  .'.tnji.turcs." 
Phys.  Rev.  8.  vol  32.  pp.  H)43M()60,  |OS5. 

D  Ahn  and  S.  L  Chuang.  ‘ '  Et act  calculatiotis  ^  jI  qii.isiboc^id  uatc.s 
Ol  an  isolated  quantum  well  with  unit.irm  eiectnc  ndJ  Quanium-weli 
Stark.  rer.«*nancc."  Phvs.  Rev.  B.  \cA  34.  np  9034-9<)3"'.  l')K6 
.A.  .\tiller.  DAB  .Miller,  and  S  D  Smiih.'  Ovnamic  non  1  near 
optical  processes  in  semiconductors. “  AJ\.  P/iyi..  vol  3t).  pp  697- 
800.  198! 

N.  Peygha/nb.irian  and  H.  M.  Gibbs.  “Opticjl  n-)ni:r.canfy.  hi.sijbil- 
ity.jnd  sicnal  processing  m  semiconductors.  J  Opt  Soc.  iiner  B. 
vol.  2.  pp”  i:i5-i:27.  1985. 

D.  S.  ChemJa  and  D  A.  B.  .Miller.  “Room-temperature  exr.s  .uc 
nonlinear-opticai  c'.Tccis  in  semiconjuctui  quanium-wdi  siiuciui..*' 
J.  Opt.  Site.  Arncr.  B.  vol.  2.  pp  1I55-I1T3.  1985 
S-  Schmirr-Rmk.  D.  S  Chemla,  and  D  .A  S  .Miller.  “Ther.  of 
transient  exciionie  optical  nonlincantici  in  semiconductor  qua'.  ..  n- 
well  structure,  ‘  Phvs  Re'..  B.  vol.  32.  pp  6601-6609,  1985. 

V.  C.  Chang.  "Nonlinear  optical  properties  of  semieondueioi-  o- 
pcriatiices."  J.  Appl.  Phys.,  vol.  58.  pp.  499-509.  1985. 

D.  D  Co«>n  and  R.  P.  G.  Kurunasin.  '  New  mode  of  [R  deiec’.on 
using  quantum  weds.  "  Appl.  Phys.  Lett  .  vol  45.  pp  649-651  i  i 


Doyeoi  Ahn  wjs  bom  In  .Seoul.  Korea,  on  iu’ 

tl96t.'  He  received  :hc  B  S  and  .M.S  degree 
electronic  engineering  from  Seoul  National  I 
vtTMty.  .Seoul.  Korea,  in  1983  and  1985.  respee 

He  IS  now  working  toward.s  the  Ph  D.  degree 
in  electrical  engineering  at  the  University  of  lib’ 
nois  at  Urbana-Champatgn  While  m  Korea,  he 
held  a  Korea  Science  and  Engineenng  Foundatioti 
Fellowship  and  Teaching  Assistantship  He  is  no. 
a  Re.scarch  A.s.sistant  and  a  GTE  Fellow.  He  is 
currently  conducting  research  on  the  electronic  and  optical  propcrijes  of 
quantum  well  structures  including  linear  and  nonlinear  eicctrooptical  ab 
•.orption  and  res<*nant  tunneling. 

Mr.  Ahn  i.s  a  member  of  T.iu  Beta  Pi.  Phi  Kappa  Phi.  and  (he  Ameiican 
Phs  steal  S<k:iciv 


Shun-Lien  Chuanu  (S'78-M'32).  lor  a  photograph  and  biography,  see  p 
509  of  the  May  198?  is.sue  of  this  Journal 


53 


APPENDIX  F 

Electron  states  in  two  coupled  quantum  wells — A  strong  coupling-of-modes 
approach 

S.  L.  Chuang  and  B.  Do 

Univenily  of  IlUiiois  at  Urhana-Chanipaign,  Department  of  Electrical  and  Computer  Engineering,  Urhana, 
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In  this  paper,  coupled-mode  equations  for  two  strongly  coupled  quantum  wells  are  derived  via 
two  methods:  (i)  a  “reciprocity”  formulation  similar  to  the  reciprocity  theorem  in 
electromagnetics,  and  (ii)  a  variational  formulation.  It  is  shown  that  both  formulations  give 
the  same  coupled-mode  equations.  Numerical  results  for  the  electron  energ_,  levels  and  wave 
functions  for  two  coupled  asymmetric  or  symmetric  quantum  wells  are  presented,  using  both 
the  strong  coupling-of-modes  theory  and  the  conventional  weak  coupling-of-modes  theory,  and 
are  compared  with  those  of  the  exact  calculations.  The  two-coupled-quantum-wells  model  is 
useful  in  the  study  of  electron  transport  in  superlattice  via  hopping  or  tunneling. 


I.  INTRODUCTION 

Interesting  properties  of  semiconductor  superlattice 
and  multiple  quantum-well  devices  have  been  under  exten¬ 
sive  investigation'  since  the  proposal  of  Esaki  and  Tsu.’  A 
two-coupled-quantum- wells  model  has  been  used  in  the 
study  of  hopping  conduction’  and  phonon-assisted  tunnel¬ 
ing  in  superlatiice  structures.''  A  coupled-mode  theory  is 
useful  for  the  calculations  of  the  energy  levels  and  w'ave  func¬ 
tions  when  the  electron  states  in  the  quantum  wells  are  cou¬ 
pled.  The  coupled-mode  theory  has  been  used  in  various  re¬ 
search  fields  such  as  the  studies  of  optical  waveguide 
couplers  and  filters,'  distributor  feedback  lasers  or  nonlinear 
optics,'’  and  bond  orbitals  in  solids.’ 

A  conventional  coupled-mode  theory  usually  assumes 
that  the  coupling  is  weak  and  a  perturbation  approach  can  be 
taken  to  find  the  coupling  coefficients  and  the  coupled-mode 
equation.  When  the  coupling  is  strong,  since  the  wave  func¬ 
tion  in  each  individual  waveguide  in  optical  couplers  or  in 
each  orbital  in  solids  is  not  orthogonal  to  each  other,  one 
needs  to  modify  the  coupled-mode  equation  by  including  the 
overlap  integrals  due  to  overlap  of  waveforms  from  each  part 
of  the  coupled  physical  system.  The  nonorthogonality  of  the 
wave  functions  contributing  to  the  nonzero  overlap  integral 
has  been  considered  in  the  bond-orbital  model  in  Refs.  7  and 
8  in  the  study  of  the  energy  levels  using  the  hybrid  wave 
functions  of  the  atoms.  However,  since  most  of  the  theoreti¬ 
cal  parameters  are  usually  adjusted  to  fitting  the  experiment 
data,  all  of  the  corrections  due  to  the  nonzero  overlap  inte¬ 
grals  are  thereby  absorbed.  On  the  other  hand,  the  improve¬ 
ments  of  the  coupled-mode  theory  in  the  optical  waveguides 
have  come  only  recently'  ''  where  the  vector  nature  of  elec¬ 
tromagnetics  and  the  overlap  integrals  are  taken  into  ac¬ 
count  It  is  shown  '  ' '  that  power  conservation  may  be  violat¬ 
ed  significantly  if  the  overlap  integral  is  not  taken  into 
account  when  coupling  is  strong 

In  this  paper,  a  strong  cotipliiig-nf-modes  derivation  for 
the  electron  states  in  two  coupled  qiianluin  wells  similar  to 
that  of  the  oplical  waveguides  is  presented  using  .i  general 
“reciprocity”  rclalii.n  and  a  variational  formulation-  Both 


formulations  give  the  same  coupled-mode  equations.  A 
model  of  two  coupled  asymmetric  quantum  wells  is  then 
considered  numerically  This  model  has  been  used  in  the 
study  of  perpendicular  conductions  in  quantum-well  struc¬ 
tures  with  an  applied  electric  field  ’  ■*  Numerical  results  us¬ 
ing  the  coupled-mode  theory  for  the  electron  energies  and 
wave  functions  are  illustrated  for  both  symmetric  and  asym¬ 
metric  wells  with  various  separations  between  the  wells.  The 
results  are  also  compared  with  those  from  the  direct  numeri¬ 
cal  method  and  those  of  the  w’eak  coupling  of  modes  theory 
neglecting  the  overlap  integrals. 

II.  STRONGLY  COUPLED-MODE  THEORY  FOR 
QUANTUM  WELLS 

In  this  section,  we  derive  a  general  “reciprocity"  rela¬ 
tion  from  the  Schrodinger  wave  equation  in  Sec.  II  A  for  two 
wave  functions  tA|(r,r)  and  t(';(r,t)  corresponding  to  two 
Hamiltonians  //,  and  //j.  This  relation  is  very  similar  to  the 
general  reciprocity  theorem  in  electromagnetics.  'We  then 
derive  the  coupled-mode  equations  using  this  quantum  me- 
chani'-al  reciprocal  relation  in  Sec.  II  B.  A  variational  for¬ 
mulation  is  also  presented  in  Sec.  II  C,  which  shows  that 
identical  coupled-mode  equations  are  obtained  using  the  two 
formulations. 

A.  A  quantum  mechanical  “reciprocity"  relation 

Consider  two  wave  functions  i/',  (r,t)  and  il’fr.t)  which 
arc  solutions  to  the  Schrodinger  equations  with  Hamilto¬ 
nians  //,  and  //,,  respectively; 

//,!/■,=  (1) 
(It 


Multiplying  (  1 )  by  the  complex  conjugate  of  I'c  and  sub¬ 
tracting  the  complex  coniugatc  of  (  2  i  multiplied  by  I'y,  and 
integrating  ovei  the  whole  space,  we  have 


n  I,'-,)  tit-  {t!\  i',). 
ft 
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'Ji  •>  -J-  ■>  ■->■ 


where  the  conventional  notations  for  the  ket,  bra,  and  the 
inner  product  are  used,  e.g.. 


,>  =  JiA?(r, 


t)tp^(T,t)d^r. 


(4) 


Equation  (3)  is  the  general  “reciprocity”  relation  which  is 
very  similar  to  the  reciprocity  theorem  in  electromagne¬ 


tics. ‘Mf 


//,  =  -  {fr/2m*)V^+  E,(r,r), 
Hj  =  -  («V2m*)V2 
Eq.  (3)  reduces  to 


(5) 

(6) 


V2(t,t)\lb^)  =  /fr— 
dt 


In  the  one-dimensional  case,  we  have 

J  -  x 


(7) 


id- 

dt 


=r 


(8) 


which  is  an  exact  relation  if  0,  and  xl!2  satisfy  the  Schrddinger 
equations  in  ( 1 )  and  (2)  and  the  corresponding  boundary 
conditions  with  the  Hamiltonians  and  //2.  Equations  (7) 
and  ( 8 )  can  also  be  used  when  is  a  slight  perturbation 
from  // 1 ,  e.g. ,  due  to  the  electron-photon  or  electron-phonon 
interactions.  One  may  have  a  stationary  potential 
y^(r,i)  =  K,(r)  with  a  set  of  solutions  tpy^ird)  to  (1) 
known  and  treat  the  difference  k'2('‘>^)  ~  as  a  time- 

dependent  perturbation.  The  unknown  wave  function 
(/’jCr.t)  can  then  be  expanded  in  terms  of  a  linear  combina¬ 
tion  of  the  set  of  solutions  with  time-dependent  co¬ 

efficient  Cj  it)-  The  time  evolution  of  Cjit)  can  then  be  in¬ 
vestigated.  The  Fermi  golden  rule  can  also  be  derived 
following  this  procedure. 

B.  Derivation  of  coupled-mode  equations  using  the 
reciprocity  relation 

Since  the  choice  of  ip,  and  1P2  can  be  arbitrary  as  long  as 
they  satisfy  the  corresponding  Schrddinger  equations  and 
boundary  conditions,  we  consider  the  following  cases. 

Case  (1): 


and 


We  choose 

Ip,  =(Paix)e 

(9a) 

F,(x)  =  yjx), 

(9b) 

t/': --<3„(x)e 

( lOa) 

y.ix)  =  F,(x), 

(10b) 

(/■:,  ^  Ki, 

or  in  matrix  form 

CE  EC  -  K  ■  K, 


K) 


(lla) 


(11b) 
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(c)  v(i) 


FIG.  1.  (a)  A  single  asymmetric  quantum  well  a  of  width  described  by 

A  single  asymmetric  quantum  well  ^  of  width  described  by 

(jc).  (c)  Two  coupled  quantum  wells  described  by  y(x). 


with 

0 

.  0  Bi,j 

and  the  matrix  elements  are  defined  as 


E- 


(12) 


r 


c/x  =  (/[/■), 


(13) 


K(x)  -  y/x)]c6.  dx^.{i\y-  y^D),  (i4) 


where  F^(.v)  and  F,,  (x)  are  single  quantum  wells  as  shown 
in  Figs.  1(a)  and  1  (b).  The  solutions  (x)  and  (x)  are 
well  known  in  texts  on  quatum  mechanics  (see  the  Appen¬ 
dix).  For  simplicity,  we  assume  here  the  effective  ma.ssof  the 
electron  is  the  same  everywhere  in  space.  Substituting  (d) 
and  ( 10)  into  (8),  we  obtain  an  analytical  relation 


where  i.j  —  a  or  b.  The  superscript  +  denotes  complex  con¬ 
jugate  and  transpose  of  the  matrix.  We  note  that  F„(x), 
F^(x),  £■„,  and  E,,  considered  here  are  all  real,  i.e.,  quasi¬ 
bound  states  such  as  those  in  a  leaky  quantum  well"’  are  not 
considered.  The  potential  energy  F(x)  here  is  chosen  to  be 
that  of  the  two  coupled  quantum  wells  as  shown  in  Fig.  1(c). 
The  differences  yix)  —  y^ix)  and  F(x)-F,,  (x)  are 
shown  in  Fig.  2. 

Case  ( 2 ) : 

Weehoose  1/1,  i.x.t)  to  be  the  solution  of  the  two-coupled 
quantum  wells  F(x): 

i/’,~a(i)(p^,{x)  4  A(/)(3„(x), 
r.fx)  =  ('■(x). 


t/'.  d>„ix)<’ 


S  i  C'  uai'a  .ind  B  Do 


-56 


V 

u 


\v 


(15a) 

f 

(  i5b) 

O' 

(  1  ba ) 

I29t 

1 

(Q)  V(»)-Vo(») 


FIG.  2.  (a)  The  difference  between  potential  energies,  F{jc)  —  V^ix).  (b) 
The  difference  between  potential  energies  F(jc)  —  (x)  used  in  calculate 

ing  the  coefficients  K,j . 


y2(x)  =  yAx).  (16b) 

We  obtain  the  first  coupled-mode  equation: 

^  =  £,  [C,,u(/)  +  C,,6(/)  J 

+  K;,a(t)  +Kt,bU). 

(17) 

Case  ( 3 ) : 

We  choose  i/*!  and  K,  (a:)  the  same  as  ( 15a)  and  (15b), 
and 


=  (18a) 

y,(x}=yAx),  (18b) 


and  obtain  the  second  coupled-mode  equation; 

^  +  C,„  =  E,[C,M<)  +  c,Mn  ] 


+  K:,a(n+K*b(t). 


Equations  (17)  and  (19)  can  be  put  in  matrix  form 


/fiC-- 

di 


a{l) 

b(t) 


Q 


a(t) 

bU) 


(19) 


(20) 


where 

e  =  EC  K  * 


=  CE  +  K 


(21) 


using  (lib).  We  note  that  E  is  real,  and  Q  and  C  arc  Hermi- 
tian.  Equations  (20)  can  also  be  put  in  the  form 


.j.d\ci{t) 

ifi~~  , 
d!  b(l) 


-  M 


all) 

hd) 


(22) 


where 

M  =  C-'Q 


=  1 

E-bC-'K 

(23) 

M  = 

y«  ' 

(24) 

n  . 

Ya  = 

E^  +  (Ar„a 

-  C^bXb. 

)/(l  -C„,C^ 

). 

(25a) 

n  = 

+  (^il) 

- 

)/(l  -Q,C^ 

). 

(25b) 

= 

{K.b  -  C 

ab^bb  )/(  1 

~  ^ab  ('be  )  ’ 

(25c) 

-C 

baK.A/0 

•c 

1 

(25d) 

C.  Derivation  of  coupled-mode  equations  using  a 
variational  approach 

The  variational  method  has  been  very  useful  in  calculat¬ 
ing  the  energy  and  wave  functions  of  the  electrons.’  It  was 
also  applied  to  the  waveguide  couplers  for  lossless”  and 
lossy  systems.'’  Consider  the  Schrodinger  equation 


HiP{x,t)  =  -^^^(x,t), 
i  dt 

(26) 

where 

zi  If  d^ 

2..  ,/x- + 

(27) 

and 

i/>(x,t)  =(/>(x)e-“^'''*. 

(28) 

Here  T(x)  is  the  potential  energy  of  two  coupled  quantum 
wells.  A  variational  formula  for  the  energy  E  is  given  by 


11 

(29) 

Note  that  T(x)  is  stationary.  Let 

ip{x,t)~-a(t)(l>Ax)  +  b[t)(l>i,{x) 

(30) 

or  use  ( 28 ) 

<fi(x)  ^a4>^(x)  +b<pi,ix). 

(31) 

where  a  and  b  are  constants. 

H  =  +  [V{x)  -  f^Jx)] 

(32a) 

+  [y(x)  -  (^.(x)], 

(32b) 

„  ff  d^  „  ,  s 

(32c) 

u  if  d^  ,  s 

(32d) 

We  obtain  ( using  tj ,  ee i/,  a.  =  6 ) 

(33) 

where  C,^  and  Q,,  are  identical  to  those  defined  in  (13)  and 
(21),  and  both  are  Hermitian  matrices.  By  taking 
dE  /da*  =  0  in  (33),  we  obtain 


^  (-^4) 

/  / 

or  in  tornis  of  the  coupled  mode  equation,  using 
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rw  Hvyj  Vw 


iftl.  ^E, 

dt 


we  obtain 


a(/)l  ^  Qp(') 

.A(r)J  L6(/)J’ 


which  is  identical  to  (20). 

III.  SOLUTIONS  TO  STRONGLY  COUPLED-MODE 
EQUATIONS 

A.  General  solutions 

A  general  solution  to  the  coupled-mode  equation  can  be 
obtained  using  (20)  by  setting 

a(r)  =  ae-'^"*,  (37) 

where 


[0(0 

~L6(/)J’ 


'  '  LM/)J’ 

and  solving  the  eigenequation 

£Ca  =  Qa,  (39) 

which  is  the  same  procedure  used  to  diagonalize  the  matrices 
C  and  Q  simultaneously.  Since  C  and  Q  are  Hermitian,  the 
eigenvalues  are  all  real  and  the  eigenvectors  satisfy  the  orth¬ 
ogonality  property” 

a'-''"Ca"'  =  0  {orE2^E„  (40) 

where  a''’  is  the  eigenvector  of  (39)  with  an  eigenvalue  E,. 
In  terms  of  the  matri.x  elements  of  M  =  C  ~  'Q,  we  have  the 
eigenvalues 

E±  =[(rh +yJ/2]±yl^^TK^  (4i) 

where 

^  =  (n-ra)/2  (42) 

and  the  corresponding  eigenvectors  (un-normalized) 

'll  [a  -f-  .y/A^  4-  1 

a"=  ^  (43) 

for  E  and 


[  -  A  -  \ 

for  E  We  keep  minus  signs  in  front  of  and  k,,^  since  in 
our  examples  to  be  shown  later,  k^^  and  k,^  are  negative. 
The  above  two  eigenvectors  can  be  easily  normalized  by  di¬ 
viding  (43)  and  (44)  by 


D"’ = /a”' +Ca'" 

(45) 

1 

1  and 

1  Z>'’'  =  v'af”^Ca'’’. 

(46) 

1  respectively. 

1  LTsing  the  conservation  of  probability  in  the  whole  space 

'  for 

1  —  (  \1’* { x.t)dx  ^  0 

1  dtj 

(47) 

by  setting  I'|(a)  —  VAx)  -  I’ix)  in  (8) 

for  two  coupled 

wells  where 
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ili(x,t)~a(i)(f>„(x)  +  bit)i/)^ix)  (48) 

and  the  coupled-mode  equation  (20),  one  finds  the  condi¬ 
tion  for  (47)  to  be  satisifed  is  that  Q  matrix  is  Hermitian.  Q 
is  indeed  Hermitian  in  this  formulation  as  has  been  shown  in 
Sec.  II.  The  statement  Q,,^  =  Qt,  with  Q  =  CM  also  leads  to 

(rt  -y.)C^  =k^  -  k%.  (49) 

Substituting  (25a)-(25d)  into  (49)  and  using  (11a)  also 
shows  that  (49)  is  satisfied  by  our  coupled-mode  formula¬ 
tion  ( 22 ).  For  most  cases,  ,  ki„,  can  be  chosen  to 

be  real  as  shown  in  the  Appendix.  We  have 

>^ba  -  Kh  =  {yh  ~  y,)C^,  (50) 

which  is  the  condition  for  probability  conservation,  which 
can  also  be  shown  by  substituting  (43)  and  (44)  into  <40). 

B.  Special  case  for  two  identical  quantum  wells 

When  the  two  wells  are  identical  and  symmetric,  we 
have 

C^=C,  (51a) 

K.b=K^sK,  (51b) 

£„=Ei,^E,  (51c) 

and 

r=y.  =  yb=E+  -  cat  )/(i  -O,  (52a) 

k  =  k,,  =k^  =  {/(-  CK,, )/( 1  -  C)  <0.  (52b) 

The  eigenvalues  due  to  the  coupling  of  states  £„  and 
will  be 

E,  =y±\k\=y  +  k^E^-^  +  -^.  (53) 

1  4-  C  1  +  c 


The  energy  splitting  is 

.  ^  ^  ^  ,  ,  -2K  +  2CA„, 

ii.E  =  E  ,  -E  _  =  2  A-  = - - 

'  1  -C 

The  normalized  eigenvectors  are  found  to  be 

1  rn 


/2(l4C)  n 


a'-'  =  — -  ,  ■  (55b) 

v'2(l  -  C)  ^  -  'J 

The  above  results  have  also  been  derived  in  the  method  of 
molecular  orbitals'"  to  study  the  bonding  and  antibonding 
states  of  the  electron  between  the  nuclei  for  two  identical 
atoms. 

Since  both  C  and  are  small  in  the  weak  coupling 
limit, 

Af - 2A:(  1  4-  C  •’ )  4  2CA'„„ ,  ( 56 ) 

while  the  leading  term  --  2K  is  the  well-known  result  for 
energy  splitting  in  a  conventional  perturbation  approach.' 
One  notes  that  this  result  ( 56)  is  also  restricted  to  two  identi¬ 
cal  wells.  The  normalized  eigenvectors  for  the  weak  coupling 
of  modes  theory  are 
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(51a) 

(51b) 
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(51c) 

(52a) 

(52b) 
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FIG.  3.  Two  wave  functions  in  two-coupled  asymmetric  quantum  wells 
with  F,  =  0. 1 5  eV.  =  0.0  eV,  K,  =  0. 1 7  eV.  =  0.02  eV.  F,  =  0. 19  cV. 
Solid  lines:  exact  solution.  Dashed  lines:  strongly  coupled-mode  theory. 


'  =  -L[  ‘ 

H  -  ll’ 


<l>,  (x)  =  [4>a(x)  +  4>k(x)]/yj2, 


(^2M  =  [<f>Jx)  -<l>,(x)]/,l2,  (58b) 

which  are  well-known  results  when  <j>^  (x)  and  (x)  are 
identical  and  orthonormal  to  each  other. 

IV.  NUMERICAL  RESULTS 

In  this  section,  we  compare  the  results  of  the  strongly 
coupled-mode  theory  to  those  of  the  exact  numerical  meth¬ 
od.  In  Fig.  3,  we  show  the  wave  functions  of  two  coupled 
asymmetric  quantum  wells  which  are  used  in  the  model  to 
calculate  the  phonon-assisted  tunneling  between  quantum 
wells.'*  We  can  see  very  good  agreement  between  the  exact 
solutions  (solid  curves)  and  the  coupled-mode  theory 
(dashed  curves).  The  parameters  are  K|=0. 15  eV, 
Fj  =  0.0  eV,  Fj  =  0. 1 7  eV,  K,  =  0.04  eV,  and  K,  =  0. 19  eV, 
In  Figs.  4-7,  we  consider  two  coupled  symmetric  quan¬ 
tum  wells.  In  Fig.  4(a),  we  plot  the  energy  levels  versus 
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FIG  4.  (a)  Energy  levels  of  two-coviplcd  symmetric  quantum  wells  are 

plotted  vs  the  separation  between  two  wells  r.  V,  -  Vs  =0.3  cV.  FIG.  5.  Wave  functions  of  iwo-coupled  quantum  wells  in  Fig.  4(a)  at  a 
—  Fj  -0.0  eV,  d„  -  (If,  -  20  A.  Solid  lines;  exact  s<>)ution.  Dashed  fixed  spacing  (a)  /  20  A.  (b)  i  --  30  A.  Solid  lines,  exact  solutions 

lines:  strongly  coupled-mnde  theory.  Dolled  lines:  weakly  coupled-mode  Dashed  lines:  strongly  couplcd-imxle  theory.  Dotted  lines,  weakly  coupled 


theory,  (b)  Percentage  errors  (»f  energies  for  (a) 


mode  theory. 
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the  separation  t  between  two  identical  wells  of  width  20  A 
and  barrier  height  0.3  eV.  The  results  of  the  strongly  cou¬ 
pled-mode  theory  (dashed  lines)  agree  well  with  those  of  the 
exact  calculation  (solid  lines).  The  results  of  the  conven¬ 
tional  weakly  coupled-mode  theory  (ignoring  the  overlap 
integrals  and  )  are  also  shown  (dotted  lines)  which 
agree  well  with  the  exact  solutions  when  the  separation  t  is 
large,  ( thus  the  coupling  is  weak ) ,  and  they  deviate  from  the 
exact  solutions  at  small  t.  The  percentage  errors  of  the  ener¬ 
gy  levels  and  ^_/E_  are  shown  in  Fig.  4(b)  for 

the  strongly  coupled-mode  theory  (dashed  curves)  and  the 
weakly  coupled-mode  theory  (dotted  curves).  One  can  see 
significant  improvements  of  the  strongly  coupled-mode  the¬ 
ory.  In  Fig.  5  ( a ),  we  plot  the  wave  functions  of  the  two  states 
at  a  separation  of  20  A.  The  wave  functions  are  all  normal¬ 
ized.  The  solid  curves  are  obtained  for  the  exact  solution,  the 
dashed  curves  are  for  the  strongly  coupled-mode  theory  and 
the  dotted  lines  are  for  the  weakly  coupled-mode  theory.  In 
Fig.  5(b),  we  increase  the  separation  to  30  A.  One  sees  better 
agreement  among  the  three  methods. 


■■ifPflRflriON  BfTWtEN  TWO  WfiLS  (Xt 


F-'Kj.  6.  fa)  Energy  ievcK  of  two-coupled  symmetric  quantum  wc)Is  are 
plotted  vs  the  separation  between  two  wells  i  E,  —  E,  -  E,  eV, 
E,  --  E,  -  0.0  eV.  o',,  -  d„  --  30  A.  Solid  lines,  exact  solution.  Dashed 
lines  strongly  coupled-mode  theory  Dotted  lines:  weakly  coupled-mixlc 
theory. 
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FIG.  7.  Wave  functions  of  two-coupled  quantum  wells  in  Fig.  6(a)  at  a 
spacing  (a)  /  =  20  A,  (b)  (  =  .30  A.  Solid  lines:  exact  solutions.  Dashed 
lines:  strongly  coupled-mode  theory.  Dotted  lines:  weakly  coupled-mode 
theory. 


In  Figs.  6  and  7,  we  consider  cases  similar  to  Figs.  4  and 
5,  except  the  widths  of  the  two  quantum  wells  are  increased 
to  30  A.  The  general  agreements  between  the  three  methods 
are  better  because  the  confinement  of  the  wave  functions  in 
each  well  is  better,  thus  the  coupling  is  weaker  for  a  fix  t  than 
that  of  Figs.  4  and  5. 

In  Figs.  8  and  9,  we  consider  two  coupled  asymmetric 
quantum  wells  with  F,  =  0.3  eV,  Fj  =  0.0  eV,  F,  =  0.32  eV, 
F»  =  0.02  eV,  and  F;  =  0.34  eV.  Again,  the  exact  solutions 
(solid  lines),  the  strongly  coupled-mode  theory  (dashed 
lines),  and  the  conventional  coupled-mode  theory  (dotted 
lines)  agree  better  with  each  other  when  the  separation  t 
increases.  In  the  limit  when  two  wells  approach  each  other 
(/— 0),  the  percentage  of  errors  for  the  conventional  theory 
is  20%  and  above  while  the  strongly  coupled-mode  theory 
has  errors  of  5%-8%  .  Large  discrepancies  between  the  wave 
functions  of  the  weakly  coupled-mode  theory  (dotted  lines) 
and  the  exact  solutions  (solid  curves)  are  shown  in  Fig.  ‘^{a) 
when  the  separation  of  two  wells  is  20  A.  while  the  results  of 
the  strongly  coupled-mode  theory  agree  well  with  the  exact 
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FIG,  8.  (a)  Energy  levels  of  two-coupled  asymmetric  quantum  wells  are 
plotted  vs  the  separation  between  two  wells.  F|  =  0.3  eV,  F,  =  0.0  eV, 
y,  =  0.32  eV,  y^  =  0.02  eV,  F,  =  0.34  eV.  30  A.  Solid  lines: 

exact  solution.  Dashed  lines:  strongly  coupled-mode  theory.  Dotted  lines: 
weakly  coupled-mode  theory,  (b)  Percentage  errors  of  energies  for  (a). 
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APPENDIX 

Evaluations  of  the  parameters  and  K,/ 

The  normalized  wave  functions  for  quantum  well  a 
[Fig.  1(a)]  is  given  by 

fe"'-'  cos  x<0 

!  cosCAr^x -t- 0., )  Q<x<d^ 

-  -  d,j  )  x>d^, 

(Al) 


=C„ 


where 

a,  =V(2mV^")(K,  -£•„),  (A2) 

-  K,),  (A3) 


(lOVV'el 


solution.  The  agreements  are  better  for  the  three  theories 
when  the  two  well  widths  and  the  separation  are  increased  to 
30  A  as  shown  in  Fig.  9(b). 

V.  CONCLUSIONS 

A  strongly  coupled-mode  theory  for  two  coupled  quan¬ 
tum  wells  has  been  presented.  Two  formulations  using  a 
quantum  mechanical  “reciprocity”  relation  and  a  variation¬ 
al  formulation  are  illustrated  and  they  give  identical  cou¬ 
pled-mode  equations.  The  presentation  here  is  very  similar 
to  the  electromagnetics  case  for  two  coupled  optical  dielec¬ 
tric  waveguides."”  Numerical  results  for  the  energy  levels 
and  wave  functions  of  the  electrons  in  the  coupled  quantum 
wells  are  illustrated  and  compared  with  those  from  the  exact 
numerical  calculation.  It  shows  clearly  that  the  overlap  inte¬ 
grals  Cab  and  should  be  taken  into  account  when  the 
coupling  is  strong. 


0  ?5  SO  lOO 

horizontal  dimension  (A)  • 


FIG.  9.  Wave  funclions  of  two-coupled  quantum  wells  with  I’,  =  0.3  eV. 
(,'^^0.0  eV.  !'■,  ^0.32  eV.  -  0,{)2  cV.  i;-^  0  34  eV.  (a) 
~t  -  20  A  (b  I  f/,  --  /  d,.  ^  30  A.  Solid  lines,  e.xiict  solution 

Dashed  lines  strongly  ctiuplcd-rnodo  theory  Dotted  lines;  weakly  coupled¬ 
mode  theory 
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^V-V^ 


I 


a,  =  ^/2^^Wy(T7=X).  (A4) 

=  -  tan'^'(a,/fcj,  (A5) 

Q  ^  ^ - - - .  (A6) 

d,  +  (1/a,)  +  (1/aj) 

The  energy  is  obtained  by  solving  the  eigenequation 

=  tan“'(a|/A:„  )  +  t&n  ' {a^/k^  )  4-  nir,  (A7) 

where  n  =  0  here  for  the  ground  state.  The  normalized  wave 
function  for  quantum  well  b  (Fig.  1(b)]  is  similarly  ob¬ 
tained. 


e“'‘'-*>cos(94 
4>^  =  C^'  cos(k^ix-s)  -f 


s<jc<s  -(- 


-  a^(x  —  df,  —  s) 


cosik^di,  -f-  &!, )  x>s  -I-  d^. 


~  yj  -  E^) ,  (A9) 

k,  =  yJ(2m*/fi^)(E,  -y,),  (A  10) 

a,==^(2m*/fi‘){y,-E,).  (All) 

and  similar  equations  to  (A5)-(A7).  Using  the  wave  func¬ 
tions  (Al),  (A8)  and  the  definitions  for  Q , /(.y  in  (13)  and 
( 14),  it  is  straightforward  to  evaluate  these  integrals  either 
analytically  or  numerically.  We  note  that  we  have  used 
Coo  =  ^bb  =  1  lo  determine  the  normalization  constants 
and  C(,  in  ( A 1 )  and  ( A8 ) .  A  good  way  to  ci.cck  tf  -  '.c.-neri- 
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cal  results  is  to  use  the  analytical  relation  (11a),  which 
should  be  satisfied  to  the  desired  accuracy  if  F, ,  C,j ,  and  K,j 
are  evaluated  correctly. 
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THE  ELECTRIC-FIELD  DEPENDENCE  OF 
INTRASUBBAND  POLAR  OPTICAL  PHONON 
SCATTERING  IN  A  QUANTUM  WELL 
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University  of  Illinois  at  Urbana-Champaign 
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ABSTRACT 

The  electric-field  dependences  of  Irtrasubbcnd  polar  optical  phonon  scat¬ 
tering  for  electrons  and  holes  in  a  semiconductor  quantum  well  are  studied  the¬ 
oretically  using  a  simple  Infinite  well  model.  It  Is  found  Chat  the  polar 
optical  phonon  scattering  rates  are  enhanced  with  an  applied  electric  field. 
They  are  higher  for  heavy  holes  for  all  ranges  of  electric  -fields ,  and  the 
electric-field  dependence  of  the  hole-polar  optical  phonon  scattering  Is  much 
stronger  than  that  of  the  electron-polar  optical  phonon  scattering.  The  higher 
subbands  have,  in  general,  weaker  electron-phonon  scattering  rates  and  electric 
field  dependences  than  those  of  the  ground  state.  The  tunneling  effect  in  a 
finite  well  is  also  discussed.  It  is  suggested  that  recent  experimental  result 
of  Che  field-dependent  line  broadening  of  near  band-edge  optical  absorption  can 
be  attributed,  at  least  qualitatively,  to  the  dominance  of  heavy  hole-polar 
optical  phonon  interaction  and  heavy  hole-tunneling. 


(PACS  index:  71.38.*-!,  72.10.Di,  72. 20. Op,  63.20.Kr). 
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I.  INTRODUCTION 


There  are  many  stimulating  works  on  the  semiconductor  quantum-well  struc¬ 
tures  because  remarkable  changes  in  electronic  and  optical  properties  are 
possible. Since  the  polar  optical  phonon  (POP)  scattering  is  the  major  scat¬ 
tering  mechanism  In  a  polar  semiconductor  like  GaAs  over  a  considerable  tem¬ 
perature  range POP  scattering  In  a  two-dimensional  structure  has  been  the 


subject  of  study  In  some  recent  papers. 


There  has  been  discrepancy  between 


some  of  the  early  calculations  of  POP  scattering  in  two-dimensional  layers'* 

Q 

and  the  experimental  results  which  show  a  weakening  of  the  electron-phonon 
Interaction  In  a  two-dimensional  structure  compared  to  that  for  a  three- 
dimensional  structure.  The  discrepancy  Is  partially  due  to  the  crude  approxima¬ 
tions  used  In  the  early  calculations,  e.g.,  q^  L  «  1  (strictly  two-dimensional 
gas),  where  q^  Is  the  transverse  component  of  the  phonon-wave  vector  and  L  Is 
the  quantum-well  width.  ’  Recently,  Leburton  calculated  the  phonon  scat¬ 
tering  rate  for  two-dimensional  electrons  In  the  lowest  subband  of  a  quantum 

well  without  an  applied  electric  field  and  showed  chat  his  results  generally 

9  7 

agree  with  the  experimental  data  from  Chiu  et  al.  and  the  model  of  Price.  In 
this  paper,  we  study  Che  effects  of  an  applied  electric  field  perpendicular  to 
Che  quantum  well  for  both  electrons  and  heavy  holes  in  Che  ground  level  and  the 
excited  state.  Previous  calculations  for  POP  scatterings  show  that  the  scat- 
tering  rate  is  enhanced  as  the  well  width  L  is  reduced.  When  one  applies  an 


electric  field  perpendicular  to  the  quantum  well,  electrons  (or‘;  holes)  are 


pushed  to  one  side  of  the  well  thus  the  effective  well  width  Is  reduced. 


10,  1  1 


This  will  result  In  the  enhancement  of  the  POP  scattering  rate.  Our  results  show 
that  the  POP  scattering  rate  is  enhanced  with  increasing  electric  field.  However, 
there  are  remarkable  differences  between  the  scattering  rates  for  the  electrons 
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and  heavy  holes  on  the  electric-field  dependence  due  to  their  different  effec¬ 
tive  masses.  Not  only  is  the  POP  scattering  rate  dominated  by  heavy  holes,  but 
the  electric-field  dependence  of  POP  scattering  rate  is  also  stronger  for  heavy 
holes.  The  former  can  be  explained  by  Price's  /m  law  for  the  two-dimensional 
collision  rate^  and  the  latter  can  be  accounted  for  by  the  electric-field  depen¬ 
dence  of  the  wave  function  variation.  The  wave  function  variation  of  the  heavy 


hole  Is  more  pronounced  than  that  of  the  electron  wave  function  at  the  same 


electric  field. 


10,11 


Let  us  consider  an  infinite  quantum  well  with  a  width  L  (Fig.  1(a)).  When 

the  quantum  well  has  an  applied  electric  field  (Fig.  1(b)),  it  is  found  that 

(1)  the  energy  level  of  the  ground  state  lowers  and  the  first  excited  state  E2 

raises  slightly,  and  (2)  the  wave  functions  are  pushed  to  one  side  of  the  quantum 

well.  Numerical  results  for  the  energy  levels  in  terms  of  the  normalized 

energies  and  normalized  electric  field  F  are  shown  in  Fig.  2(a)  where  E^  = 

E^/Ej^°\  /(2m  L^)  and  F  »  je]  FL/Ej^°^  where  m  is  the  effective 

mass.  The  amplitudes  of  the  wave  functions  are  shown  in  Fig.  2(b).  Since  Che 

★ 

effective  mass  of  heavy  holes  (m  "  0.45  m  ,  m  is  the  free  electron  mass)  is 

0  0 


*  (0) 

heavier  than  Chat  of  the  electron  effective  mass  (m  =  0.0665  m  ),  E, 

o  1 


for  the 


heavy  hole  is  lower  Chan  that  of  the  electron.  As  a  result,  the  normalized  F 
for  a  heavy  hole  is  larger  than  that  for  an  electron  at  the  same  electric  field 
strength  F.  This  causes  the  stronger  electric-field  dependence  of  the  POP  scat¬ 
tering  in  a  two-dimensional  gas  for  heavy  holes.  Our  calculations  show  that  the 
POP  scattering  rate  for  Che  electron  is  enhanced  by  about  10%  when  the  electric 
field  F  is  200  kV/cm  for  L  »  126. 5A  compared  with  that  for  Che  zero-field  case, 
while  the  POP  scattering  rate  for  Che  heavy  hole  is  almost  doubled  when  F  =  200 
k'//cm  compared  with  the  zero-field  results.  This,  together  with  tunneling  of  the 


4 


heavy  hole,  may  explain  qualitatively  the  recent  experimental  results  where  the 
exclton  line  width  is  almost  doubled  when  ?  ■  220  kV/cm  compared  with  those  for 
the  zero-field  case.^^  We  also  found  that,  in  general,  the  POP  scattering  rates 
and  their  electric-field  dependences  for  higher  subbands  are  smaller  chan  chose 
for  the  ground  state.  Many  body  effects,  however,  are  not  considered  here  and 
will  be  a  subject  of  future  study.  The  electric  field  dependence  of  the  carrler- 
POP  scattering  should  also  be  important  in  Che  calculation  of  the  intersubband 
optical  absorption^ ^  in  a  quantum  well. 
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ir.  ELECTRIC-FIELD  DEPENDENCE  OF  POLAR  OPTICAL  PHONON 

SCATTERING  OF  TWO-DIMENSIONAL  CARRIERS  IN  A  QUANTUM  WELL 

The  Hamiltonian  of  the  system  (a  single  quantum  well)  subject  to  a  uniform 
electric  field  perpendicular  to  the  quantum  well  (the  z-direction)  with  carrier- 
polar  optical  phonon  Interaction  Is  written  as 


H  -  H  +  H'  (1) 

o 

where  is  the  unperturbed  Hamiltonian  for  an  electron  (or  hole)  in  the  quantum 

well  in  the  presence  of  a  perpendicular  electric  field.  For  simplicity,  we  as¬ 
sume  isotropic  end  t.arabolic  bands,  and  consider  only  one  type  of  hole  (heavy 
hole) . 

The  Frohllch  Hamiltonian  H'  In  SI  units  Is^** 


If  (a^ 

*-  q  q 
<l  ^ 


+ 

where  a  and  a  are  the  second  quantized  creation  and  annihilation  operators  of 

q  H 

the  phonon  with  the  wave  vector  q,  respectively,  and 


Here,  V  is  the  volume  of  the  system,  is  the  permittivity  of  free  space, 
and  are  the  optical  and  static  dielectric  constants,  respectively,  and 
Is  the  phonon  frequency,  =  36.202  meV.  Here  iWe  consider  the 

Intrasubband  transition  only,  i.e.,  the  energy  separation  of  subbands  is  suf¬ 


ficiently  larger  Chan  .  If  we  neglect  Che  interaction  between  the  electrons 
(or  holes)  in  the  well,  the  wave  functions  for  Che  initial  state  >1/'^  (r)  and  Che 
final  state  ili'^  (r)  can  be  written  as 
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where  A  Is  the  cross-sectional  area  of  the  well,  L  is  the  width  of  the  well,  a 
denotes  either  electron  or  holes,  and  k""^  are  the  wave  vectors  of  electrons 
(or  holes)  in  the  x-y  plane  for  the  Initial  and  final  states,  respectively,  and 
r^  is  the  position  vector  in  the  x-y  plane.  We  assumed  in  Eqs.  (4)  that  the 
perturbations  vary  slowly  so  that  the  cell  periodic  function  part  will  be 
Integrated  to  one  in  Che  matrix  element  to  be  evaluated  later.  Thus,  the  cell 
periodic  function  is  not  Included  here.  The  envelope  function  (z)  for  the 
nth  subband  satisfies  the  following  Schrodinger  equation  in  the  effective  mass 
approximation : 


— —  — ~  ♦  (z)  -  e  F  z  ♦  (z)  =  E  i>  (z) 

.  *  ,  z  n  n  n  n 

2m  dz 
a 


(5) 


where  F  is  the  electric  field,  the  subscript  n  refers  to  Che  nth  subband  with  the 
energy  E^,  and  m^  is  the  effective  mass  for  specie  a  Note  that  e  =  -  |ej  for 
electrons  (o  =  e)  and  e  =  +  |e[  for  holes  (a  =  h).  The  solutions  of  Eq .  (7)  are 
well  known  and  are  given  by  a  linear  combination  of  the  two  independent  Airy 
functions  Ai(n)  and 


(z  ) 
n 


a“  /'i(Ti“)  +  b**  Bi(n°) 
n  n  n  n 


<  ^ 
-  2 


(5) 


with  the  boundary  conditions 
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r  -) 


where  n  is  defined  by 
n 


-  [- 


2m 


(eflF) 


(e“  +  e  F  z) 
n 


(7) 


(8) 


The  two  constants  a**  and  b**  and  the  energy  e'*  can  be  determined  from  Eqs. 

n  n  n  ^ 

(8)-(10)  and  the  normalization  of  the  wave  function.  From  now  on,  we  will  omit 
the  superscript  a  for  convenience. 


The  transition  rate  for  POP  scattering  is  expressed  in  the  Born  approxima¬ 
tion  by  Fermi's  Golden  Rule: 

Wf5  iMf^l  5(E^  -  Ej  i  (9) 


where  E.  and  E^  are  the  total  initial  and  final  energies  of  the  carrier,  respec¬ 
tively,  and  and  denote  the  transition  rate  for  phonon  emission  and  absorp 
tion,  respectively.  The  matrix  elements  are  written  as 


M 

f  i 


-  (n  +  1)'^^  5  , 

‘‘t  -  '’t  -L/2 


/  dze'^'’"'^  I*  (z)(^ 


(10) 


and 


f  i  q  q 


C  1/2  . 

—  n  6,  ^ 


k  ,  k  +  q  r 
t  t  ^t  - 


I  i2 

/  dze  !♦  (z)  I 


(11) 


L/2 


The  initial  and  the  final  energies  E.  and  are  defined  in  the  parabolic  band 


approximation,  respectively. 


deviation  from  / m*  law  for  the  scattering  race  Is  expected.  The  numerical 


results  of  W  and  W  for  several  cases  will  be  presented  In  the  next  section, 
n  n 
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NUMERIC^.L  RESULTS  AND  DISCUSSIONS 


The  total  transition  rates  W  (E  )  and  W  (E^)  of  the  nth  subband  with  a 

n  t  n  t 

two-dimensional  wave  vector  for  phonon  emission  and  absorption,  respectively, 
have  been  calculated  by  numerically  Integrating  I^  and  I^  in  Eqs.  (14)  and  (16) 
for  the  lowest  two  subbands  with  an  applied  electric  field  F.  The  parameters 
are  as  follows: 

m  (electron)  “  0.0665  m^ ,  m  (hole)  “  0.45  m^ ,  T  -  300  K,  L  •  126. 5A 
e  "  10.92,  e  -  12.90,  and  fim  -  0.036202  (eV), 

00  *8  q 

where  m  is  the  free  electron  mass.  The  phonon  occupation  number  n  is  given  by 
o  q 

the  Bose-Einsteln  distribution. 


n^  -  {exp  [(fi<o^)/KgT]  -  l}"\  (13) 

where  K  is  the  Boltzmann  constant.  We  use  an  effective  well  width  L  “  126. 5A 

O 

chosen  to  give  the  same  ground  state  energy  of  a  finite  well  with  a  width  L^  » 
lOOA  and  the  barrier  height  V  =  340  meV  in  the  GaAs-Al  Ga,  As  system. 

In  Fig.  3,  we  plot  the  scattering  rates  for  the  electrons  and  the  heavy 
holes  with  an  applied  electric  field  (1)  F  ■  200  kV/F  (solid  lines)  and  (2) 

F  “  0  (dashed  lines)  for  the  first  subband.  Although  noticeable,  the  effects  of 
an  applied  electric  field  on  electron-POP  scattering  rates  are  generally  small 
( e lect ron-POP  scattering  rates  are  increased  by  about  102  when  F=  200  kV/cm  with 

L  =  126. 5A  compared  with  those  for  the  zero-field  case).  However,  as  seen  from 

the  figure,  the  electric  field  changes  the  heavy  hole-POP  scattering  rates 

drastically  by  almost  a  factor  of  2  when  F  =  200  kV/cm  compared  with  those  for 
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the  zero-field  case.  Higher  scattering  rates  of  the  heavy  hole  for  all  electric 

fields  can  be  attributed  to  its  heavier  effective  mass  than  that  of  the 

electron.  For  a  finite  quantum  well,  however,  tunneling  of  carriers  is  also 

important  at  high  electric  field  regime  compared  with  the  zero  field  case  for 

which  the  polar  optical  phonon  scattering  Is  the  dominant  scattering  mechanism 
10  ' 

for  T  >  200  k.  lo  compare  the  POP  scattering  rate  with  the  tunneling  rate, 

we  have  also  calculated  the  tunneling  rate  at  F  ■  200  kV/cm  from  the  complex 

energy  eigenvalues  of  quasi-bound  state  for  electron  (L  “  lOOA  and  potential 

well  height  -  340  meV)  and  heavy  hole  (L  ”  lOOA  and  »  100  meV).  Tunneling 

rates  for  electron  and  hole  in  their  ground  states  at  F  »  200  kV/cm  are  5.37  x 
9  -1  12  -1 

10  (sec  )  and  21.22  x  10  (sec  ),  respectively.  Recently,  a  band  offset 

ratio  of  60:40  for  AlGaAs/GaAs  (265  meV  for  electrons,  175  meV  for  holes)  Is 

more  commonly  accepted.  In  this  case,  the  tunneling  rates  for  electron  and  hole 

In  their  ground  states  at  F  -  200  kV/cm  are  3.05  x  10^^  (sec  S  and  7.6  x 
12  -1 

10  (sec  ),  respectively.  So  heavy  hole-POP  scattering  and  heavy  hole¬ 
tunneling  have  comparable  transition  rates.  From  these  results,  we  can  explain, 

1  2 

at  least  qualitatively,  the  recent  experimental  data  of  field-dependent  exci- 

ton  linewidth  broadening  in  a  two-dimensional  structure  by  heavy  hole-POP 

interaction  and  heavy  hole-tunneling.  We  do  not  consider,  however,  ocher 

2  0 

linewidth  broadening  mechanisms  such  as  the  interface  roughness,  Che  field 
inhomogeneity,  and  the  carrier-carrier  interaction.  In  Fig.  4  the  scattering 
races  for  electrons  and  heavy  holes  with  an  applied  electric  field  (1)  F  =  200 
kV/cm  (solid  lines)  and  (2)  F  =  0  (dashed  lines)  for  the  second  subband  are 
plotted.  In  general,  Che  scattering  rates  and  their  electric-field  dependences 
of  the  second  subband  are  smaller  than  those  of  the  ground  level.  Smaller 
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electric-field  dependences  can  be  explained  by  the  fact  that  the  wave  function 
change  due  to  the  applied  electric  field  of  higher  subbands  is  smaller  than  that 
of  the  ground  state^^. 
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IV.  CONCLUSIONS 

The  eleccric-f leld  dependence  of  the  polar  optical  phonon  scattering  rates 
I  for  both  the  electrons  and  the  holes  in  a  quantum  well  are  studied.  For  all 

i 

I  electric  field  ranges,  the  scattering  rates  for  the  holes  are  higher  than  those 

of  the  electrons  due  to  the  effective  mass  difference.  When  the  electric  field 

!  F  •  200  kV/cm,  the  scattering  rates  for  the  holes  are  almost  doubled,  while  the 

! 

scattering  rates  for  the  electrons  are  increased  by  about  lOX  only,  compared 
with  those  for  the  zero-field  case.  It  is  found  that  higher  subbands  have 
smaller  scattering  rates  and  weaker  electric-field  dependences,  for  both  the 
electrons  and  the  holes,  than  those  of  the  ground  state.  In  our  opinion,  the 
field-dependent  line-width  broadening  of  Interband  edge  optical  transitions 
could  be  dominated  by  heavy  hole-polar  optical  phonon  interaction  and  heavy 
hole-tunneling.  Other  mechanisms,  however,  such  as  the  Interface  roughness 
scatterings,  may  also  contribute  to  the  total  llnewidth. 
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APPENDIX:  DERIVATIO:.’  OF  EQS.  (13)  to  (16) 


The  total  transition  rates  W  (k  )  and  V  (k_)  for  the  carrleis  in  the  nth  sub 

n  t  n  t 

band  with  a  two-dimensional  wave  vector  k^  are 


<  ^‘‘t>  -f  ^  l”fll'  ^  -  ^1  " 

‘‘t 


|C|^  (n  +  1)  I  I  ^  /  dz  e  (♦^(2)1^1 

"  k'  q  *t’  *t  ‘’t  -L/2  ••  I 
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f."  2 
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k,  q 


- '  '"t  ♦  ’t*' '  ^7  "t'  ■  SI 

2  (A2 

V|C|  n  .2  . 

“  - 27"  ^  '’t  ^n  ^ - *  ^'^t  ^  k^q^.cos  $ )  -  fiu  ] 

^Ti  n  2a  ^ 


where  is  the  azimuthal  angle  between  k^  and  q^.,  -  ti  £  (f  £  it  •  Here,  (q^,L,F) 
is  defined  by 


G„(q^.L.F)-/  dq^-  -  /  dze  |*^(z)r 

^z  ‘'t  I 


and  can  be  evaluated  by  contour  integration.  The  result  is 

^  /  .  ..V  r .  1  I  .  -l^zz  ,.  , 
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^z  ^t 


2  -qt|z-2 
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The  evaluation  of  the  delta  function  is  straightforward  and  is  given  by 
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2 
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In  Eq.  (A6),  the  condition  qt  ^  0  makes  the  other  delta  function  with  an  argu 
ment  q^  +  cos  ((>  +  /  cos^  <))  +  2b_  vanish. 

2  2 

For  the  emission  process  (Eq.  A5),  the  conditions  q^  ^  0,  k^.  cos 

* 

0)  >  0  and  cos^  (^  >  1  give 

h  q  -  -  “ 


^max  ^  ^max 
E  /lid)  >  1 

t  q 


^  =  cos  ^  )  and  E  »  - — ^ 

^max  V  ^  t  . 


The  condition  E  /fuu  >  1  gives  the  two-dimensional  emission  rate  a  steplike 
t  q 

behavior.  For  Che  absorption  process  (Eq.  (A6)),  all  values  of  ^  between  -  n 

and  7t  are  allowed.  Using  Eqs.  (Al)  to  (A8),  we  obtain  the  total  transition 

rates  (k. )  and  W  (k  )  : 
n  c  n  c 

W'^  (k  )  =  0)  (n  +1)  a  (E_  F,  L)  (AS 

n  t  q  q  epn  t’ 
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where  we  have  reduced  the  ♦  integration  to  0  < 
since  cos  is  an  even  function  of 


with  a  factor  of  2  added 
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and  a  is  defined  in  equation  (17). 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  diagrams  tor  the  energy  levels  and  the  wave  functions  of  an 
Infinite  quantum  well  for  (a)  F  «  0  and  (b)  F  >  0. 

Figure  2.  (a)  The  normalized  subband  ener^gles  (E^  “  E^/Ei^o))  are  plotted  vs. 

the  normalized  electric  field  (F  ■  |e|  FL/Ej^o').  (1)  The  exact  solu¬ 
tion  (solid  lines),  (2)  the  variational  method  with  the  Gram-Schmidt 
orthogonalizatlon  procedure  -  (dashed  lines). 

(b)  The  magnitudes  of  the  wave  functions,  L|i|;jj(z)p,  are  plotted  vs. 
the  normalized  distance  z/L  for  the  first  three  subbands.  The  dotted 
lines  are  for  the  zero  electric  field.  The  solid  lines  are  the  exact 
solutions  and  the  dashed  lines  are  those  for  the  variational  method 
with  the  Gram-Schmldt  orthogonalizatlon  procedure  for  a  quantum  well 
with  an  applied  electric  field  (F  «  20).  (From  Ref.  11) 

Figure  3.  Plots  of  the  polar  optical  phonon  scattering  rates  of  electrons  and 
holes  for  an  applied  electric  field  (1)  F  =  200  kV/cm  (solid  lines) 
and  (2)  F  ■  0  (dashed  lines)  for  the  ground  state  with  a  well  width 
L  »  126. 5A  at  T  *  300  K. 


Figure  A.  Plots  of  the  polar  optical  phonon  scattering  r^tes  of  the  electrons 
and  the  holes  for  an  applied  electric  field  (1)  F  =*  200  kV/cm  (solid 
lines)  and  (2)  F  »  0  (dashed  lines)  for  the  second  subband  with  a 
well  width  L  =  126. 5A  at  T  =  300  K. 
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ELECmiC-FIELD  DEPENDENCE  OF  THE  INTERSUBBAND  OPTICAL 
ABSORPTION  IN  A  SEMICONDUCTOR  QUANTUM  WELL 

D.  Ahn  and  S.  L.  Chuang 

Department  of  Electrical  and  Computer  Engineering 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  Illinois  61801 


We  present  theoretical  results  of  intersubband  linear  optical  absorption  in  the 
conduction  band  of  a  GaAs-AlGaAs  quantum  well  with  an  applied  electric  field  taking 
into  account  the  field  dependent  linewidih.  Our  analysis  is  based  on  the  one  electron 
density  matrix  formulation  with  incrasubband  relaxation  processes  due  to  polar  optical 
phonon  scattering  and  tunneling  of  electrons.  We  show  that  (a)  for  an  increasing 
electric  field  the  absorption  peak  corresponding  to  the  transition  of  states  1  — >  2  is 
shifted  higher  in  energy  and  (b)  the  peak  amplitude  increases  if  the  Fermi  level  is  fixed 
and  decreases  if  the  electron  density  in  the  well  is  fixed  when  an  increasing  electric  field 
is  applied.  The  linewidth  broadening  also  reduces  the  peak  absorption  amplitude  . 


'  K." 


1.  Introduction 

The  existence  of  discrete  energy  levels  by  quantum 
confinement  of  carriers  in  semiconductor  quantum  well 
systems  results  in  resonance  optical  transitions  of  infrared 
radiation  of  suitable  polarizations  between  quantized 
subbands.  The  electric-field  dependence  of  the 
intersubband  optical  absorption  is  of  growing  interest 
because  of  potential  device  applications  including  high 
speed  infrared  optical  modulators  and  photodetectors.  6-11 
The  Stark  shifts  of  intersubband  optical  absorption  spectra 
observed  experimentally  6,7  agree  with  the  theoretical 
calculations  of  the  Stark  shifts  of  subband  energy  levels  in 
a  quantum  well  structure  8-10.  Previous  calculations  12-15 
show  that  the  energy  level  of  the  ground  state  Ej  decreases 
and  the  first  excited  state  E2  increases  slightly  when  there 
is  a  moderate  applied  electric  field.  Therefore,  the  energy 
difference  E2  -  Ei  corresponding  to  the  peak  energy  of 
intersubband  absorption  is  increased  when  the  electric  field 
is  increased.  For  the  peak  absorption  amplitude  change 
with  increasing  electric  field,  two  different  experimental 
results  were  reported,  however.  Harwit  and  Harris  6 
observed  that  the  peak  amplitude  decreased  with  an 
increasing  electric  field  .  Bajema  et  al.  ^  showed  that  the 
peak  amplitude  increased  with  increasing  electric  field 
which  agreed  with  our  earlier  calculations  8,9  which 
assumed  a  fixed  Fermi  level  Ep  .  As  we  will  show  later, 
the  difference  between  the  results  of  the  two  experiments 
may  be  due  to  the  difference  in  the  experimental 
conditions.  Both  experiments  show  an  increase  in  the 
energy  of  peak  absorption,  which  agrees  with  the 
theoretical  results.  It  is  well  known  that  the  optical 
absorption  coefficient  depends  on  the  dipole  matrix 
element,  the  density  of  states  (or  population  of  the  initial 
and  the  final  states)  ,  and  the  linewidth.  For  intersubband 

transition,  the  dipole  matrix  element  ^  ^  V2  '  ez  ^  Vp  ^  varies 
slowly  with  applied  electric  field,  so  a  major  change  of  the 
absorption  pe^  could  be  attributed  to  the  change  of  the 
population  and  the  linewidth  broadening,  which  may 
strongly  depend  on  the  condition  of  the  experiment  .In  this 
paper,  w'e  calculate  ilie  linear  intersubband  electro-optical 
absorption  of  a  quar.'jm  well  model  for  two  different 
hypothetical  conditions;  (1)  the  Fermi  level  is  assumed  to 
be  fixed  and  (2)  the  electron  density  in  the  well  is  assumed 
to  be  fixed.  We  show  that  (1)  if  the  Fermi  level  is  fixed, 
the  peak  absorption  amplitude  increases  since  more 
electrons  will  occupy  the  ground  state  with  an  applied 
electric  field  (E]  is  low'ered)  and  (2)  if  the  electron  density 
is  fixed,  the  peak  absorption  amplitude  decreases  slightly 
since  the  optical  dipole  maunx  elements  decrease  slightly 
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and  the  electron  population  difference  change  is  negligible 
when  linewidth  broadening  is  not  taken  into  account  The 
linew'idth  broadening  if  included  will  funher  reduce  the 
absorption  coefficient.  Our  analysis  is  based  on  the  one 
electron  density  matrix  fonnulation  with  intrasubband 
relaxations  .  We  improve  our  previous  calculations, 
8-10  5y  taking  into  account  the  electric-field  dependent 
relaxation  processes  due  to  the  polar  optical  phonon 
scattering  and  tiie  tunneling  of  quasi-bound  states. 

2.  Theory 

Let  us  consider  an  electron  with  a  charge  -lei  and  an 
effective  mass  m*  in  a  quantum  well  with  a  well  width  Lf 
and  a  depth  Vq  in  the  presence  of  an  applied  electric  field  F 
,  an  optical  radiation  with  an  angular  frequency  of  w  and 
the  polarization  along  the  direction  of  the  well  z,  electron- 
polar  optical  phonon  interaction  ,  and  tunneling  of  quasi¬ 
bound  states  (Fig.  1)  .  We  choose  the  origin  to  be  at  the 
center  of  the  well.  To  calculate  the  optical  matrix  element 
and  the  polar  optical  phonon  scattering  rates,  we  use  an 
infinite  quantum  well  model  with  an  effective  well  width  L 
chosen  to  give  the  same  ground  state  energy  of  the  finite 
quantum  well.  Tunneling  of  the  quasi-bound  states  of  the 
finite  quantum  well  with  an  applied  electric  field  is  taken 
into  account  in  the  linewidth  broadening  of  the  absorption 
coefficient.  The  intersubband  optical  absorption  coefficient 
cx  for  the  transition  between  the  lowest  two  levels  is 


given  by 


i  |2 

'  (E2-Ei-'hm)V(1i7,b) 


vath  the  matrix  element 


(.L/2 

Lj  =  lei  62  (z)  z  pi  (z)  dz 
•'-L/2 


where  p  is  the  permeability;  nr  is  the  refracitive  index  (for 
GaAs.  nj^3.2);c  is  the  speed  of  light  in  free  space;  V  is  the 
volume;  kt  is  the  wave  vector  of  the  electron  in  the  x-y 
plane;  Ei  and  E2  denote  the  quantized  energy  levels  for  thie 

initial  and  final  states,  respectively;  Oj  and  are  the 
envelope  functions  for  the  first  and  the  second  lowest 
states,  respectively;  a  and  b  ate  the  short  hand  notations  for 
the  total  v,'ave  vectors  of  the  initial  and  final  states. 

respectively;  fi  is  the  P’anck  c'nstant  divided  by  2  n  t'h'fab 


'n-’-.-'s”'- 


is  the  relaxation  rate:,  and  fa  and  fb  are  the  Fermi-Dirac 
distribution  functions  for  the  initial  and  final  states, 
respectively.  The  overall  factor  of  2  in  front  of  eq.  (1) 

denotes  the  two  spins.  The  intrasubband  linewidth  -fiTab 
is  defined  by 

- - T,;;  I  ■— -r—  j  (3 

‘.-ph.ii  "e-ph/'  t,b  / 


Wticic  l/'-c-pii.ii  aiivj  i/',.  lUc  iJi^  i^iai  cic^Uuii-polar 
optical  phonon  scattering  rates  for  the  states  a  and  b, 

respectively,  and  l/^i,a  and  are  the  tunneling  rates 
for  the  states  a  and  b,  respectively.  The  electric  field 
dependences  of  the  polar  optical  nhonon  scattering  rates 
have  been  calculated  previouly  and  the  tunneling  rates 
can  be  obtained  from  the  imaginary'  pans  of  the  complex 
energy  eigenvalues  of  quasi-bound  states  .  If  we 
change  the  summation  over  into  the  integration  over  the 
transverse  energy  Et ,  eq.  (1)  becomes 


coucm*  r“._  IM2,l^(f,-f.)hYjE,) 


n^ErtTi"  •'O  (E2-Ei-1ico)‘-fetiy,b(E,))‘ 


where  Et  is  the  transverse  component  of  the  total  energy 
defined  bv 


Ej=(k^+  ky  )/(2m  )  (5) 

We  calculate  the  linear  intersubband  optical 
absorption  for  two  cases,  (1)  the  Fermi  level  Ep  is  fixed 
and  (2)  the  electron  density  in  the  quantum  well  is  fixed^  S. 
The  intergration  over  E;  in  eq.  (4)  is  done  numerically. 

The  parameters  used  in  our  calculations  are  as 
fellows: 

m'“=0.0665mo,  T=300  K,  Lf=100  A, 

L=  126.5  A,  Vo=340  meV.  (6) 

* 

We  used  an  effective  well  width  L=126.5  A,  chosen  to 
give  the  same  groupd  state  energy  of  a  finite  well  with  a 
width  of  Lf  =  100  A  and  the  barrier  height  Vq  =  340  raeV 
in  the  Ga.4.s-.41Ga.4s  ,s}  siem,  to  calculate  the  optical  dipole 
matrix  element  and  the  electron-polar  optical  phonon 
scattenng  rates.  Here  mo  is  the  free  electron  mass  and  T  is 
the  temperature.  The  electron  density  n  in  the  quantum 
well  and  the  Fenni  level  Ep  are  related  bv 


m  KpT 


where  Kb  is  the  Boltzmann  constant  and  Ei  is  the  i  th 
quantized  subband  energy  level.  In  the  calculation  of  case 
(1),  we  used  the  value  of  the  Fermi  level  Ep  =-40.47  meV 
for  n=3.0xl0i^/cm3  when  electric  field  F  =  0  and  then  fix 
the  Fermi  level  when  a  field  F  is  applied.  If  Tab  were 
assumed  to  be  a  constant^' iO  ,  eq.  (4)  can  be  integrated 
analytically  and  cx  turns  out  to  be  proportional  to  the 
population  difference  in  the  two  quantized  levels,  nj-nj.  If 
Ep  is  fixed  when  F  >  0,  ni-n2  will  increase  since  E]  is 
lowered.  Thus  one  sees  that  the  peak  absorption  ot  will 
increase  and  E2  is  raised.  In  Fig.  2,  we  plot  the  average 

polar  optical  phonon  scattering  rates  <  1/  t  >  g.ph 

(solid  lines)  and  the  tunneling  rates  <  1/x  >  tunneling 
(dashed  lines)  versus  applied  electric  field  F.  It  can  be  seen 
readily  from  the  figure  that  average  polar  optical  phonon 
scattering  rates  rise  verj'  slowly  with  the  electric  field  F 
and  the  tunneling  rate  of  the  second  quantized  state 
increases  rapidly  when  F  is  greater  than  100  kV/cm  and 
becomes  comparable  with  the  polar  optical  phonon 
scattering  rates  when  F=200  kV/cm.  Our  previous 
calculations shows  that  the  electron-phonon  scattering 
rate  is  increased  by  about  10%  and  the  heavy  hole-phonon 
scattering  rate  is  almost  doubled  at  F  =  200  kV/cm 
compared  with  those  of  t^e  zero  field  case-0  in  a  quantum 
well  of  width  L=126.5  A.  In  Fig.  3,  we  plot  the  linear 
iniersubband  optical  absorption  coefficient  ot  versus  the 
incident  photon  energy  for  the  case  (1)  where  the  Fermi 
level  Ep  is  fixed  to  the  zero  field  value.  The  absorption 
coefficients  are  plotted  for  (i)  F=0  (dotted  line  ),  (ii)  F=200 
kV/cm  without  tunneling  (dashed  line),  and  (iii)F=200 
kV/cm  with  tunneling  (solid  line).  It  can  be  seen  that  the 
peak  amplitude  is  increased  in  genera!  with  increasing 
electric  field  F.  In  Fig.  4,  we  plot  the  linear  intersubband 
optica]  absorption  coefficient  for  case  (2)  when  the  electron 
densi'v  in  the  quantum  well  is  fixed  to  3.0  x 
when  an  electric  field  F  is  applied  .  The  absorption 
coefficients  are  plotted  for  (i)  F=0  (dotted  line) ,  (ii)  F=200 
kV/cm  without  tunneling  (dashed  line),  and  (iii)  F=200 
kV/cm  w'ith  tunneling  (solid  line).  The  peak  amplitude  is 
decreased  with  increasing  elecmic  field  F.  In  both  cases, 
the  absorption  peak  is  shifted  higher  in  energy  due  to  the 
quantum  confined  Stark  effect,  which  agrees  with  both 
experiments  in  refs.  6  and  7.  The  full  widths  at  half 
maximum  (FW'HM)  of  the  iniersubband  absorption  arc  3.0 
meV  (F=0).  3.28  meV  (F-200  kV/cm  without  tunnelintt), 


and  5.31  meV  (F=200  kV/cm  with  tunneling).  We  can 
deduce  that  the  linewidth  due  to  the  polar  optical  phonon 
scattering  is  increased  by  10%  when  F=200  kV/cm 
compared  with  zero  field  case  and  tunneling  contributes 
2.03  meV  when  F=200  kV/cm.  If  we  include  other 
linewidth  broadening  mechanisms  such  as  interface 
roughness  and  impurity  scatterings,  we  would  have  further 
broadening  of  FWHM.  From  the  results  of  cases  1)  and 
2)  and  experimental  obser\'ations  ^•'7,  it  seems  that 
experimental  conditions  may  affect  the  population  change 
when  an  external  electrip  field  ic  applied.  In  ref.  7,  the  well 
is  relatively  large  (264A)  and  the  field  is  relatively  small 
(<12.5kV/cm),  thus  the  tunneling  effect  and  other  field 
induced  broadenings  may  be  ignored.  In  ref  6,  however, 
the  well  is  narrow'er  and  the  field  is  larger;  thus  the 
linewidth  broadenings  due  to  the  electric  field  may  be 
important,  thus  reduce  the  peak  absorption  amplitude. (In 
both  experiment  .however, the  linewidth  broadening  due  to 
the  tunneling  is  negligible,  so  other  mechanisms  may  be 
responsible  for  the  linewidth  broadening  of  ref6.).  The 
pesik  amplitude  is  mainly  determined  by  the  linewidth  and 
the  electron  population  difference  which  need  funher 
study. 


3.  Conclusion 

Electric-field  dependence  of  the  linear  intersubband 
optical  aborption  is  calculated  for  two  cases  :  (1)  fixed 
Fermi  level  and  (2)  fixed  electron  density  in  the  well.  Polar 
optical  phonon  scatterings  and  tunneling  of  quasi-bound 
states  are  considered  for  linewidth  broadening  .  It  is 
shown  (a)  the  absorption  peak  is  shifted  higher  in  energy 
with  increasing  electric  field,  (b)  the  absorption  pe^ 
amplitude  increases  with  increasing  electric  field  for  case 

(1)  and  decreases  with  increasing  electric  field  for  case 

(2) .  Further  study  on  other  linewidth  broadening 
mechanisms  such  as  the  interface  roughness  scatterings, 
the  web  width  inhomogeneity,  and  the  impurity  scatterings 
are  necessary  for  the  interjmetation  of  the  experimental 
results. 
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Figure  Captions 


Fig.  1 .  Linewidth  broadenings  of  the  intersubband 
absorption  due  to 

(a)  electron-phonon  scattering  and 

(b)  electron  tunneling  through  the 
barrier  in  the  presence  of  an  applied 
electric  field. 

Fig.  2.  Average  polar  optical  phonon  scatter¬ 
ing  rates  (solid  lines)  and  tunneling 
rates  (dashed  line)  are  plotted  for  the 
lowest  two  subbands  versus  the  electric 
field  F. 

Fig.  3.  The  absorption  coefficients  arc  plotted  for 

(i)  F=0  (dotted  line),  (ii)F=200  kV/cm  without 
tunneling  (dashed  line),  and  (iii) 

F=200  kV/cm  with  tunneling(solid  line) 
vs.  the  photon  energy  for  the  case  of  the 
fixed  Fermi  level  An  effective  well 
widthL=126A  A,  T=300K,  and  Ef= 

-40.47  meV  (for  n=3xl0l6/cm  3  when  F=0)  ai 
used. 

Fig.  4.  The  absorption  coefficients  are  plotted  for 
(i)F=0(dotted  line),  (ii)_F=2(X)  Wcm 
without  tunneling(dashed  line),  and  (iii) 

F=200  kV/cm  with  tunneling(solid  line) 
vs  photon  energy  for  the  case  of  the  fixed 
electron  density  of  3xl0^6/cm3 
at  T=300K  for  a  quantum  welj,  with  an 
effective  well  width  L=  126.5  A. 
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TUNNELING  EFFECTS  ON  THE  INTERSUBBAND  ELECTRO-OPTICAL 
ABSORPTION  IN  A  QUANTUM  WELL 
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Abstract 


The  electric-field  dependence  of  the  Intersubband  optical  absorption  in  a 
quantum  well  with  llnewidCh  broadened  by  the  tunneling  of  electrons  has  been 
studied  theoretically.  We  show  that  for  an  increasing  electric  field  the  line 
width  broadening  is  enhanced  by  the  tunneling  of  quasi-bound  state  electrons, 
which  reduces  the  absorption. 


'PACS  index:  78.20.Dj,  73.40.Gk,  78.20.Jq) 
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The  optical  transition  between  quantized  subbands  (Intersubband  absorption) 

1-3 

vjii;  >jf  the  most  remarkable  properties  of  quantum-well  structures 

tereir  )  ■  ti  e  .ileci  ilc  field  dependence  of  the  Intersubband  optical  absorption 

4  5 
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experimentally  agree  with  the  trend  of  the  theoretical  calculations  of  the 

electric-field  dependence  of  subband  energy  levels  In  a  quantum-well  struc- 
6'“9  8  9 

ture  .  Previous  calculations  ’  show  that  the  energy  level  of  the  ground 
state  decreases,  and  the  first  excited  state  E2  increases  slightly  when  there 
is  an  Intermediate  applied  electric  field.  Thus,  the  energy  difference  E2  -  Ej 
corresponding  to  the  peak  energy  for  optical  absorption  is  Increased  when  the 
electric  field  F  is  increased. 


For  the  peak  amplitude  variation  with  the  electric  field,  two  different 

4 

experimental  results  were  reported.  Harwlt  and  Karris  observed  that  the  peak 
amplitude  is  decreased  while  the  optical  absorption  llnewidth  is  significantly 
broadened  with  increasing  electric  field  and  Bajema  et  al.^  showed  that  the  peak 
amplitude  '.ncreases  while  the  linewidth  is  not  significantly  affected  by  the 
applied  electric  field.  In  the  work  by  Bajema  et  al.,  wide  quantum  well  (L  = 

264  A)  and  low  field  (F  <  12.5  kV/cm)  are  used;  thus  the  tunneling  effect  and 
field  induced  linewidth  broadening  are  not  important.  However,  in  Harwit  and 
Harris’s  work,  they  use  relatively  narrow  well  (L  =  120  A)  af  intermediate  field 
range  (F  -  36  kV/cm),  where  linewidth  broadenlngs  appear  to  be  significant 
(tunneling  is  still  not  important  for  this  well  width).  The  purpose  of  this 
paper  is  to  study  the  effect  of  tunneling  especially  for  narrow  quantum  wells. 
The  dipole  matrix  el  jnent  |<>J’2  j  decreases  only  slightly  with  increasing 
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electric  field  since  the  wave  functions  are  pushed  to  the  same  side  of  the  quan¬ 
tum  well  (while  the  electron  and  hole  wave  functions  are  pushed  to  the  opposite 
lues  of  the  well  when  1 nter hand  transition  was  considered  previously^^) .  Thus, 
i  !  ,  -  i.ignltule  of  the  1  nt  o  rsiibl-und  ahsoiprlon  coefficient  would  decrease  slightly 
i'  ^  f  ri  ■  ill.  •:  of  electroiis  is  .issamed  to  be 

constant  in  the  well.  The  llnewidth  broadening  due  to  the  tunneling  effect  and 
other  scattering  mechanisms  would  further  decrease  the  absorption  peak  ampli¬ 
tude.  On  the  other  hand,  if  the  Fermi  level  is  assumed  to  be  fixed  and  the 
llnewidth  broadening  is  neglected  when  an  electric  field  is  applied,  the  popula¬ 
tion  difference  of  the  first  two  subbands  will  increase,  thus  increasing  the 


absorption  coefficient 


6,7 


For  a  finite  quantum-well  structure  with  an  applied 


electric  field,  it  is  well  known  that  the  system  does  not,  strictly  speaking, 

11  12 

have  true  bound  states  ’  .  Instead,  it  has  quasi-bound  states  with  a  complex 

energy  E  »  Eg  -  ir/2,  where  E^  and  F  correspond  to  the  quasi-bound  state  energy 
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level  and  the  resonance  width,  respectively  Our  previous  calculation  showed 
that  the  resonance  width,  corresponding  to  the  tunneling,  rises  sharply  with 
increasing  electric  field.  Since  the  tunneling  also  affects  the  llnewidth  for 
the  optical  absorption,  the  llnewidth  of  the  intersubband  optical  absorption  may 
be  broadened  considerably  by  the  electric  field.  Thus,  in  a  narrow  quantum  well 
for  which  tunneling  becomes  more  important  when  an  electric  field  is  applied, 
the  magnitude  of  the  intersubband  optical  absorption  coefficient  is  expected  to 
be  reduced  with  increasing  electric  field. 

In  this  paper,  we  present  theoretical  calculations  for  the  electric-field 
dependence  of  the  intersubband  optical  absorption  taking  into  account  the  tun¬ 
neling  of  quasi-bound  states.  To  calculate  the  optical  matrix  element,  we  use 
an  infinite  quantum-well  model  with  an  effective  well  width  that  gives  the  same 
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ground  state  energy  of  the  finite  quantum  well  .  Tunneling  of  quasi-bound 
states  of  the  finite  quantum  well  with  an  applied  electric  field  Is  then  taken 
into  account  in  the  llnewidth  broadening  of  the  absorption  coefficient. 

The  elect rl c-f 1  eld  dependence  of  other  llnewidth  broadening  mechanisms, 

such  as  electron-polar  optical  phonon  scattering,  Is  neglected  for  Intersubband 

transition,  because  Its  dependence  on  the  field  is  small  for  conduction  band 
13 

electrons  .  For  1 nterband  transition,  however,  it  is  found  that  heavy  hole- 

polar  optical  phonon  scattering  gives  the  same  order  of  transition  rate  as 

13 

heavy-hole  tunneling,  which  is  different  from  the  Intersubband  case 
Inclusion  of  other  scattering  mechanisms  would  give  further  broadening  of  the 
llnewidth.  Theoretical  analysis  of  the  effects  of  tunneling  on  the  intersubband 
electro-optical  absorption  Is  important  to  fully  understand  the  principles  of 
recently  proposed  far-infrared  photodetectors^^  based  on  the  Intersubband 
absorption  with  tunneling-assisted  impact  Ionization. 

Let  us  consider  an  electron  with  charge  -!e|  and  effective  mass  m*  in  an 
infinite  quantum  well  with  an  effective  well  width  L  (which  gives  the  same  ground 
state  energy  of  a  finite  well  of  width  L^  and  depth  V^)  in  the  presence  of  a 
constant  electric  field  F  along  the  positive  direction  of  the  well  z.  We  choose 
the  origin  to  be  at  the  center  of  the  well.  We  consider  a  quantum  well  with  an 
applied  electric  field  In  the  presence  of  an  optical  radiation  with  an  angular 
frequency  id  and  the  polarization  along  the  direction  of  the  well  z.  For  the 
arbitrary  polarization  case,  there  is  an  additional  factor  cos  6  in  the  matrix 
element,  where  0  is  the  angle  between  the  z-axis  and  the  direction  of  the  polari¬ 
zation^  for  intersubband  transitions.  The  intersubband  optical  absorption  coef- 

•^6  7 

flcient  a  for  the  transition  between  the  lowest  two  levels  is  given  by"’  ’ 


■Hi 


(fioi  -  E2  -  Ej)*-  +  (r/2) 


where 


m*K  T  E  -  E 

n  = - Y~  (~-K-T~')] 


and  the  matrix  element 


|e|  / 


(z)z  4(,  (z)  dz  , 


where  E^  and  E^  denote  the  quantized  energy  levels  for  the  initial  state  and  the 
final  state,  respectively,  p  is  the  permeability,  c  is  the  speed  of  light  in 
free  space.  Kg  is  Boltzmann's  constant,  T  is  the  temperature,  n^  is  the  refrac¬ 
tive  index,  mQ  is  the  free-space  electron  mass,  E^  is  the  Fermi  energy,  fi  is  the 
Planck  constant  divided  by  2ir ,  41^  and  <^2  ^re  the  envelope  functions  for  the 
first  and  second  lowest  states,  respectively,  and  P  is  the  llnewidth.  It  is 
seen  clearly  from  Eq.  (1)  that  the  absorption  coefficient  due  to  transitions 
between  the  first  two  subbands  is  determined  by  (1)  the  dipole  matrix  element 

[M2j|,  (11)  the  electron  population  difference  (n^  -  02)  of  the  two  subbands  and 

3 

(ill)  the  llnewidth.  The  linewldth  P  can  be  written  phenomenologically 


2  =  -T  "  I  "  7;^ 


where  Pg  is  the  llnewidth  in  the  absence  of  the  electric  field,  Xj  and  T2  are 
the  tunneling  life  times  of  the  first  and  the  second  quasi-bound  states,  respec¬ 
tively.  The  tunneling  rates  Tx/  x,  and  Ti/  x-  can  be  obtained  from  the  imaginary 


parts  of  the  complex  energy  eigenvalues  of  quasi-bound  states  .  The  llnewldth 

broadening  due  to  tunneling  can  also  be  understood  by  the  broadening  of  the  den- 
19 

slty  of  states  near  the  energy  of  the  original  bound  state  when  there  is  no 
leakage  of  electrons.  In  the  presence  of  an  applied  electric  field,  there  will 
be  new  available  states  due  to  the  leakage  of  the  electrons  or  the  coupling  of 
the  free  states  on  one  side  of  the  well  with  the  bound  state  in  the  well. 


The  parameters  used  in  our  calculations  are  as  follows; 

m*  =  0.0665  m^,  T  -  77  K,  L  =  101.27  A,  n  =  5.0  x  10^^/cm^,  and  Tq  =  10  meV 

We  used  an  effective  well  width  101.27  A,  chosen  to  give  the  same  ground  state 
energy  of  a  finite  well  with  a  width  of  =  65  A  and  the  barrier  height  = 
245  meV  in  the  GaAs-AlGaAs  system.  The  electron  density  n  and  the  Fermi  level 
Ep  are  related  by 

n  =  I  "i 
1 


m*KBT 


-  E  '  E  “I 

I  in  i  +  exp  , 

i  L  J 


and  we  assume  that  only  the  lowest  two  subbands  are  occupied  and  the  electron 
concentration  n  is  fixed  when  there  is  an  applied  electric  field.  (If  we  fix 
the  Fermi  level  Ej,  at  the  value  of  the  zero  field,  then  due  to '  the  changes  of  E 

and  E^,  the  logarithmic  function  in  Eq.  (1)  will  increase  with  increasing 
electric  field,  as  will  the  absorption  coefficient^.  This  result  agrees  with 
that  of  the  experiment  in  ref.  5,  where  the  llnewldth  broadening  effects  are 
negligible  since  the  well  width  is  relatively  wide  and  the  electric  field  is 


.  V 


A.  v.  <■.  ir.  rf-_  . 
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small.  For  that  cas';,  the  carrier  concentration  n  will  be  changed.)  The  zero- 

field  llnewldth  Tq  =  10  cieV  is  taken  from  the  experimental  results  .n  ref.  1. 

In  Figure  1,  we  plot  the  tunneling  life  times  (in  picoseconds)  of  the  electrons 

for  the  ground  state  n  =  1  (solid  line)  and  the  first  excited  state  n  =  2 

(dashed  line)  versus  the  electric  field.  Tunneling  life  times  have  been  calcu- 

1  2 

lated  from  the  complex  energy  of  the  quasi-bound  states  .  From  this  figure  we 

can  see  that  electrons  in  the  second  quantized  level  tunnel  through  the  barrier 

more  rapidly  than  the  electrons  in  the  ground  state  because  the  barrier  for  the 

excited  electrons  Is  lower  than  that  of  the  ground  state  electrons.  In  Figure 

2,  we  plot  the  half  llnewldth  r/2  given  In  Eq.  (3)  versus  the  electric  field  F, 

The  llnev^idth  is  doubled  when  F  =  60  kV/cm  for  this  particular  quantum  structure 

and  is  increasing  almost  linearly  after  F  =  50  kV/cro.  In  Figure  3,  we  plot  the 

absorption  coefficient  a  versus  incident  photon  energy  fiw  (eV)  for  both  electric 

field  F  =  0  (dashed  line)  and  60  kV/cm  (solid  line).  The  magnitude  of  the  peak 

absorption  coefficient  is  reduced  by  half  when  F  =  60  kV/cm.  It  can  be  seen 

that  both  Stark  shifts  of  the  peak  absorption  energy  and  the  linewidth 

A 

broadening  show  the  same  trend  observed  by  Harwit  and  Harris  .  One  should  note 
that  the  total  linewidth  in  the  experiment  would  include  not  only  the  tunneling 
effect  and  the  scattering  effect  but  also  other  mechanisms  such  as  the  Interface 
roughness  scatterings  and  the  well  width  Inhoraogenei ty  effect.  Our  calculation 
for  L  =  120  A  used  in  ref.  h  shows  that  tunneling  is  negligible  for  that  case 
since  the  well  is  still  considerably  wide.  For  the  potential  device  applica¬ 
tions  such  as  infrared  photodetectors^^  if  the  applied  electric  field  is  too 
strong  and  the  well  is  narrow,  the  optical  absorption  will  be  reduced  and  the 
efficiency  of  the  device  will  be  affected. 


’tjr.r.f.t 


In  conclusion,  we  have  studied  Che  tunneling  effects  on  the  InCersubband 
electro-optical  absorption  In  a  finite  quantum  well.  It  Is  shown  that  the 
magnitude  of  Che  peak  absorption  coefficient  is  reduced  by  half  when  the  field 
Is  60  kV/cm  for  the  well  width,  L  =  101.27  A  (or  finite  well  of  =  65  A  and 
Vq  =  245  A)  due  to  tunneling  .^ndnced  linewldth  broadening. 

We  are  grateful  for  usefu.  discussions  with  K.  Hess.  This  work  was  sup¬ 
ported  by  Che  Air  Force  Contract  F33615-84-K-1557 .  One  of  the  authors  (D.  Ahn) 
would  also  like  Co  thank  GTE  for  providing  support. 
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FIGURE  CAPTIONS 


Figure  1.  Tunneling  life  times  for  the  ground  state  n  =  1  Csollc  line)  and  the 
first  excited  state  n  =  2  (dashed  line)  are  plotted  versus  the 
electric  field  F. 

Figure  2.  The  half  llnewldth  r/2  for  optical  Intersubband  absorption  Is 
plotted  versus  the  electric  field  F. 

Figure  3.  The  electro-optical  absorption  coefficients  a  for  an  electric  field 

F  =  0  (dashed  line)  and  60  kV/cro  (solid  line)  are  plotted  versus  the 

incident  photon  energy  hw  (eV)  for  a  well  of  width  L  =  101.27  A  with 
16 

5.0  X  10  electrons/cm^  at  77  K. 
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A  FIELD  EFFECT  QUANTUM-WELL  LASER 


WITH  LATERAL  CURRENT  INJECTION 


D.  Ahn  and  S.  L.  Chuang 


Departm<inC  of  Electrical  and  Computer  Engineering 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  IL  61801 


Abstract 


Polarization  dependent  gain  switching  in  a  field  effect  quantum-well  laser 
with  lateral  current  injection  is  studied.  We  use  the  k 'p  method  for  the  opti¬ 
cal  dipole  matrix  and  take  into  account  the  intraband  relaxation.  Gain 
switching  is  achieved  by  the  field  effect  from  a  gate  terminal  in  a  lateral 
current  injection  quantum-well  laser  structure.  It  is  shown  that  the  peak  gain 
exhibits  a  red  shift  and  the  peak  gain  amplitude  decreases  with  an  increasing 
electric  field  for  the  TE  mode.  The  TM  mode  is  less  affected  by  the  external 
field. 
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QuaiLum-well  lasers  are  of  growing  interest  because  of  their  superior 
characteristics,  such  as  low  threshold  current  and  narrow  gain  spectrum,  com- 
pared  with  those  for  conventional  semiconductor  lasers  .  Recently,  an 
AlGaAs/GaAs  lateral  current  injection  quantum-well  laser  was  proposed  in  which 
the  structure  is  planar,  suitable  for  integration  with  other  optical  devices  and 
exhibits  built-in  index  guiding  and  a  very  low  stray  capacitance^.  The  possibil¬ 
ity  of  electric-field-induced  photoluminescence  or  gain  switching  in  a  similar 
structure  with  a  third  terminal  has  been  suggested  ’  .  Field-induced  gain 
switching  has  an  advantage  over  the  conventional  method  by  controlling  injection 

current  in  that  the  former  method  is  not  limited  by  the  carrier  lifetime  which 

6  7  8 

limits  the  switching  speed  '  '  .  Recent  experimental  results  by  T.  Takeoka  et 
3 

al,  show  that  a  short  optical  pulse  as  narrow  as  140  ps  full  width  at  half 
maximum  is  generated  by  field-induced  gain  switching  in  an  optically  pumped 
quantum-well  laser. 

9 

Arakawa  et  al .  showed  active  Q  switching  achieved  in  another  coupled  sec¬ 
tion  adjacent  to  the  active  region  of  the  conventional  quantum-well  laser 
through  the  quantum  confined  stark  effect.  There  is  no  applied  voltage  directly 
across  the  active  region  of  the  quantum-well  laser  in  the  structure  and  as  a 
result,  there  is  no  direct  control  of  the  intrinsic  gain  in  the  active  region. 

In  this  paper,  we  report  theoretical  results  for  gain  switching  in  a 
lateral  injection  quantum-well  laser  with  a  controlling  gate  similar  to  that  of 
a  field  effect  transistor  (Fig.  1).  We  call  it  a  field  effect  quantum-well 
laser  because  (i)  the  electrons  and  holes  are  injected  laterally  by  the  p  and  n 
junctions  and  they  conduct  parallel  to  the  quantum-well  interfaces,  and  (ii)  the 
gain  switching  is  achieved  by  an  applied  voltage  at  the  gate,  thus  creating  an 
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internal  electric  field  perpendicular  tt>  the  quantum-well  interfaces.  The 

7 

region  below  the  gate  can  be  an  insulating  AZ^Gaj^_^As  or  AMs  region  .  The 
electric  field  will  push  the  electron  and  hole  wavefunct ions  to  opposite  sides 
of  the  well  in  the  direction  perpendicular  to  the  well,  therefore  decreasing  the 
overlap  integral  tor  optical  transitioi.s  in  the  laser  structure  (Fig.  2).  Thus, 
the  gain  is  decreased  by  the  applied  perpendicular  field.  We  show  a  model 
calculation  of  the  electric-field  dependence  of  the  linear  gain  in  a  quantum 
well  wher?  carriers  are  injected  parallel  to  the  well.  Our  analysis  is  based  on 

the  density  matrix  formulation  with  intraband  relaxation ^  \  We  use  the  lc*p 

12  .1113 

perturbation  method  and  envelope  function  approximation  *  for  the  optical 

matrix  element  and  transition  energy. 

Let  us  consider  the  electrons  and  holes  in  a  quantum  well  with  a  well  width 
L  in  the  presence  of  an  applied  electric  field  F  along  the  direction  of  the  well 
z.  We  use  an  infinite  well  model  to  calculate  the  wave  functions  when  there  is 
an  applied  electric  field.  We  choose  the  origin  to  be  at  the  center  of  the 
well.  For  simplicity,  we  assume  isotropic  and  parabolic  bands,  and  will  con¬ 
sider  only  the  heavy  hole  in  the  valence  band.  The  effect  due  to  the  light  hole 
band  is  smaller  because  of  the  smaller  density  of  states  since  its  effective 
mass  is  smaller.  Effects  of  more  complicated  band  structures  such  as  band 
mixing  will  be  a  subject  of  future  study.  In  a  one-electron  model,  the  linear 
gain  in  a  quantum-well  laser  is  given  as 


¥%%% 


^  '■‘J  k.  (E-u  -  hu)  +  (h/i. 


\  -r.  y. 


where  u  is  the  angular  frequency  of  photons,  p  is  the  pertneab i I ity ,  n^  is  the 

refractive  index,  c  is  the  speed  of  light  in  free  space,  j)''  and  are  the 

envelope  functions  for  electrons  and  holes  in  the  i-^  and  the  subbands, 

respectively,  V  is  the  volume,  is  the  wave  vector  in  the  x-y  plane,  e  is  the 

polarization  vector,  is  the  optical  dipole  moment  matrix  element  between  the 

i.“h  subband  in  Lne  ^uuUuctLOu  Land  jiid  lue  j-*’  subband  in  the  heavy  hole  band 

where  i  and  j  dependences  of  the  matrix  elements  M  come  from  the  orientational 

properties  of  the  Bloch  functions  at  the  zone  center,  is  the  transition 

energy  between  the  i“^  subband  in  the  conduction  band  and  the  subband  in  the 

valence  band,  f''  and  f;^  are  the  Fermi  functions  for  the  i^^  subband  in  the  con- 
c  n 

duct  ion  band  and  the  subband  ir.  the  valence  band,  respectively,  and  is 

the  intraband  relaxation  time.  The  factor  of  2  comes  from  two  spin  states. 

When  F  “  0  the  optical  transition  is  possible  only  when  i  *  j  by  the  selection 

rule;  then,  Eq.  (1)  reduces  to  the  conventional  expression  derived  in  Refs.  (Q 

and  l2.  However,  when  F  *  0,  electron  and  hole  wavef unct ions  are  pushed  to  the 

opposite  sides  of  the  well  (Fig.  2)  and  the  quantum  well  becomes  asymmetric. 

Therefore,  the  overlap  integral  for  i  #  j  is  reduced  and  the  optical  transition 

between  the  subbands  i  and  j  where  i  *  j  becomes  possible.  Previous  calcula- 
13  IS 

tions  ’  show  that  the  energy  level  of  the  ground  state  decreases  for  both 
electrons  and  holes  (quantum  confined  Stark  effects).  Therefore,  the  peak  gain 
exhibits  a  red  shift  in  energy  and  a  decrease  of  the  magnitude  with  an 
increasing  electric  field  F.  The  optical  dipole  matrix  elements  are  evaluated 
using  the  J«p  method  in  the  four-band  approximation 


( I  +  cos  9^ ) 


j  )  ^  c  i  r,^ 


sin  9. 

i 


for  TE 


for  TM 


where  9^  is  Lhe  angle  between  the  k  vector  and  the  z  axis  for  the  i*-^  subband, 
!>!*■  J  I  ^  is  given  by 


Eg  (Eg  +  A ) 


3  a) 


where  e  is  the  electron  charge,  m  is  the  effective  mass  of  an  electron,  Eg  is 

c 

the  band  gap,  and  A  is  the  spin-orbit  splitting  in  the  valence  band.  The  envelope 

functions  and  are  given  by  a  linear  combination  of  Airy  funct ions^ ^ ^ ^ .  As 

a  numerical  example,  we  have  calculated  the  linear  gain  g^^^u)  in  a  GaAs-A£GaAs 

quantum-well  structure  with  an  external  electric  field.  Here  we  assume  the 

electric  field  F  is  applied  perpendicular  to  the  well  through  the  gate  and  the 

current  is  injected  parallel  to  the  quantum  well  (lateral  injection)  as  shown  in 

Fig.l.  The  screening  effect  and  the  electric  field  inhomogeneity  are  neglected. 

The  screening  effect  would  reduce  the  effective  electric  field  strength  in  the 

well  and  create  the  spatial  dependence  of  the  field^ .  We  have  assumed  that  the 

intraband  relaxation  time  t-  is  the  same  for  all  subbands  and  is  independent  of 

in 

the  electric  field.  Calculated  results  of  the  gain  spectra  for  TE  are  shown  in 

Fig.  3  for  the  well  width  L  =  100  A  with  the  intraband  relaxation  time  t.  =  1  x 

in 

- 1 3 

10  (sec)  at  T  *  300  K  for  the  cases:  F  »  0  (dotted  line),  F  “  50  kV/cm 
(dashed  line),  and  F  «  100  kV/cm  (solid  line).  We  have  assumed  that  the  den- 

18  3 

sities  of  electrons  and  holes  remain  at  2.0  x  10  /cm  for  all  F.  From  the 
figure,  it  can  easily  be  seen  that  the  peak  gain  spectra  show  a  red  shift  and 
the  peak  gain  amplitude  decreases  with  an  increasing  electric  field  F.  Figure  4 
shows  the  gain  spectra  for  the  TM  mode  for  the  well  with  the  same  parameters  as 
in  Fig.  3.  It  can  be  seen  that  the  sain  f-r  the  TM  mode  is  niuch  smaller  and 
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less  affected  by  the  field  than  that  for  the  TE  mode.  Results  in  Fi?.  3  show 
that  effective  gain  change  could  be  achieved  at  a  relatively  low  electric  field 
(<100  kV/cm).  The  field-induced  gain  change  has  an  advantage  over  the  gain 
change  by  injection  current  in  a  conventional  diode  laser  in  that  the  former  is 
not  limited  by  the  carrier  lifetime.  Therefore,  electric-field-induced  gain 
switching  can  be  used  in  ultrafast  optical  signal  processing.  The  physical 
nature  of  the  field-induced  gain  switching  is  attributed  to  the  quantum-confined 
Stark  effects. 


In  conclusion,  we  have  calculated  the  electric-field  dependence  of  the 
linear  gain  in  a  GaAs-MGaAs  field  effect  quantum-well  laser  with  lateral 
current  injection.  It  is  shown  that  the  peak  gain  shows  a  red  shift  and  the 
peak  gain  amplitude  decreases  with  an  increasing  electric  field  for  TE  modes. 
Gain  spectra  for  TM  modes  are  less  affected  by  an  applied  electric  field. 
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Figure  Captions 

Fig.  1.  SLructure  of  a  field  effect  quantum-well  laser  with  lateral  current 
inject  ion . 

Fig.  2.  The  band  diagram  for  a  quantum  well  with  an  applied  electric  field. 

F  and  F  are  the  quasi-Fermi  levels  for  the  electrons  and  holes, 
n  p  ’ 

respect ive ly . 

Fig.  3.  Linear  gain  spectra  calculated  for  the  TE  mode  for  the  cases:  F  =  0 

(dotted  line),  F  <•  50  kV/cm  (dashed  line),  and  F  =  100  kV/cm  (solid 

f  R 

line).  We  use  T  =  300  K,  the  carrier  concentration  N  =  2.0  *  10  cm 
and  the  well  width  L  =  100  A. 

Fig.  4.  Linear  gain  spectra  calculated  for  the  TM  mode  for  the  cases:  P  =  0 

(dotted  line),  F  »  50  kV/cin  (dashed  line),  and  F  =  100  kV/cm  (solid 

1  ft 

line).  We  use  T  »  300  K,  the  carrier  concentration  PJ  =  2.0  x  10^  cm 


and  the  well  width  L  =  100  A. 
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Abstract  We  present  numerical  calculations  of  the  second  harmonic  suscep¬ 
tibility  due  to  intersubband  transition  within  the  conduction  band  of  a  quantum 
well  in  an  external  applied  electric  field.  The  asymmetry  of  the  quantum  well 
due  to  Che  electric  field  accounts  for  the  nonvanishing  of  the  second  order 
susceptibilities.  It  is  shown  Chat  for  moderate  values  of  an  applied  electric 
field  of  10  kV/cm  to  70  kV/cra,  the  second  harmonic  susceptibility  is  generally 
10  to  100  times  larger  than  chat  of  bulk  GaAs.  Furthermore  this  procedure  of 
second  harmonic  generation  can  be  controlled  by  an  external  modulating  voltage. 


Quantum  confinement  of  carriers  In  a  semiconductor  quantum  well  leads  to  the 

formation  of  discrete  energy  levels  and  the  drastic  change  of  optical  absorption 

spectra^.  Optical  transitions  between  the  size-quantized  subbands  are  feasible. 

2  -3 

These  have  been  studied  experimentally  for  cases  without  an  electric  field 
and  with  an  electric  field  .  Large  dipole  strength  and  a  narrow  bandwidth  were 
observed.  This  suggests  large  optical  nonlinearities  associated  with  intersub¬ 
band  transitions^  The  linear  and  nonlinear  optical  absorption  coefficients 

for  general  quantum  well  systems  with  an  applied  electric  field  have  been  calcu- 
7  ”8 

lated  •  In  this  paper  we  study  the  electric  field  control  of  second  harmonic 
generation  In  a  quantum  well.  From  symmetry  considerations,  it  is  clear  that 
the  second  order  electric  susceptibility  is  only  nonzero  if  the  quantum  well  is 

9 

asymmetric.  Such  asymmetry  can  be  created  by  compositional  grading  or  by  an 
external  applied  electric  field.  In  the  following  we  present  calculations  of 
the  second  harmonic  susceptibility  due  to  intersubband  transitions  of  a  quantum 
well  in  an  applied  electric  field.  It  Is  shown  that  the  susceptibility  can  be 
10  to  100  times  that  of  bulk  GaAs  for  moderate  values  of  electric  field  strength. 
An  additional  advantage  is  that  this  procedure  of  second  harmonic  generation  can 
be  controlled  by  an  external  modulation  voltage. 

Consider  an  asymmetric  quantum  well  in  the  presence  of  optical  radiation 
E(t)  with  the  polarization  along  the  direction  of  the  well  z  (fig-  !)•  The 
asymmetry  can  also  be  introduced  by  applying  an  electric  fle\,d  to  a  symmetric 
quantum  well.  Cases  for  the  arbitrary  polarization  will  introduce  an  additional 
factor  cos  0  in  the  matrix  element  where  6  is  the  angle  between  the  z-axis  and 


the  direction  of  polarization  for  intersubband  transitions. 
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Let  p  be  the  one  electron  density  matrix  for  this  system,  Hq  be  the  unper¬ 
turbed  Hamiltonian  for  this  system  with  an  asymmetric  potential  energy,  and  M  be 
a  dipole  operator.  Since  we  are  considering  an  asymmetric  quantum-well  system, 
the  dipole  operator  has  nonvanishing  diagonal  elements  which  contribute  directly 
to  the  nonlinear  susceptibility.  The  one-electron  density  matrix  equation  with 


Intraband  relaxations  becomes 
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|t  “  [^0  ~  p]  -  '2-  [r(p-p^*^^)  +  (p-p^*^^)r] 


where  [  ,  ]  is  the  quantum  mechanical  commutator,  fi  is  Planck's  constant  divided 

by  2ti,  p^*^^  is  the  unperturbed  density  matrix,  and  T  is  the  phenomenological 

operator  responsible  for  the  damping  due  to  the  electron-phonon  interaction, 

collisions  among  electrons,  etc.  We  assume  that  T  is  a  diagonal  matrix  and  its 

elements  f  are  the  inverse  of  the  relaxation  time  for  the  state  ]u>.  We  use 
uu 

shorthand  notations  such  that  |u>  =  In  ,  ^  >,  where  'k  is  the  wave  vector  for 

'  '  u  tu  cu 

electrons  in  the  x-y  plane  and  n^  =  1,  2,  3,  ...  denote  the  size  quantized 

energy  eigenstates  of  the  subband.  The  Hamiltonian  H^  is  diagonal  in  Iu>  with 

F  =  V  ^  *  * 

energy  k^^/2m  with  ra  being  the  effective  mass  of  electrons  near 

the  conduction  band  edge.  The  diagonal  elements  of  the  operator  T  are  <u  |rl  u> 

=  1/t  =ry  :  Y  Y  “  (1/t  1/t  )/2.  Equation  (1)  can  be  solved  bv  pertur- 

u  uu  uv  vu  u  V  ' 

.  10-11 
ba  1 1  on 

p  =  y  p^"^  (2) 


with  the  unperturbed  density  matrix  p  assumed  to  have  only  diagonal  terms. 
Let  "he  Incident  optical  field  be  represented  as 

follows 


with  (Uj ,  j  »  1,  2,  N  as  the  angular  frequencies  of  the  optical  radiation. 


Then  the  second-order  soljtlon  of  the  density  matrix  is 
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where  W  =  (E„  -  Ey)/fi  and  M  =  <u  ]m  1  v>.  The  first  subscript  d  in 
uv  uv 

(2) 

£  refers  to  the  definition  of  (5)  for  di  f  f  erent  frequencies  u.  and 

Because  of  conservation  of  momentum  in  the  xy  plane,  the  k^'s  of  lu>  and  |v> 

must  be  the  same  for  the  matrix  element  M  to  be  nonzero.  The  problem  of  a 

uv 

\  7  8  1 

quantum  well  with  an  applied  electric  field  F  has  been  solved  ’  ’  “ .  The 
solution  is  (r)  =  <r|U>  =  <r  ]n  ,  k  >  =  A  U  (r)  exp(ik  'r^)  (}i  (z),  where 

U  t  C  tU  t  Oy 

A  is  the  area  of  the  well,  r  is  the  position  vector  in  the  x-v  plane,  and  U  i 

t  c 

the  cell  periodic  cell  function.  The  wave  function  41  (z)  can  be  expressed  in 

7  8 

terms  of  the  Airy  functions  ’  . 
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The  matrix  element  is 
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with  k  ”  k  .  The  second  order  electronic  polarization  Is 
tv  tu 
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In  the  following,  we  illustrate  the  results  for  a  two  level  system  |a> 


1,  k  >  and  lb>  =  j2,  k  >.  We  assume  Ferral-Dlrac  statistics  for  and 

t  t  aa  bb 


Also  let  the  optical  radiation  be  of  angular  frequency  w  so  that  in  (8),  Uj 


=  Hi  and  E  -  E  =•  E  . 

*■  ~  ~K 


Then,  the  second  harmonic  susceptibility  is,  from  (b 


and  (9), 
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(2u)  +  W.  +  I-y  .  )(a)  +  W  +  if 
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where  -  N2  ”  AN  Is  Che  electron  population  difference  of  the  two  levels  and  Is 


r  h'  h 

.*k,T  '  *  '"P  (  k.T  ) 

1  +  exp  ( - ^ - ) 


with  E  being  the  Fermi  level.  A  simplified  result  can  be  obtained  from  (10)  by 
F 

retaining  only  the  near-resonant  term  at  2u)  =  Wv.*  Then 


<"l  -  "2>  l"bal  <"bb  -  "ka> 

(2«  -  .  n^) 


We  note  frora  (12)  chat  d  Is  zero  if  diagonal  matrix  eJements  of  Che 
dipole  flioraenC  operator  M  are  zero,  as  in  a  symmetric  quantum  well.  In  the  pres- 
ent  case,  d^  Is  zero  when  the  applied  electric  field  F  is  zero. 


moderate  applied  electric  field  with  F  *  10  kV/cra  and  F  «  70  kV/cra,  The  values 
of  relaxation  times  are  chosen  to  be  "  0.14  ps.  We  note 


aa 


^  ab 


that  the  normalized  )d  j  can  reach  values  of  10  to  5  *  10  for  moderate 
values  of  applied  electric  field.  In  Figure  3  and  Table  I  the  calculated  peak 
values  of  |d  '^  |  in  (1/9  x  10  MKS  units)  are  illustrated  as  a  function  of 
applied  electric  field  F.  For  example,  at  F  =  50  kV/cm,  the  peak  value  Is 

5398.48  (1/9  *  10  MKS  units).  This  Is  much  higher  than  the  bulk  GaAs  value^^ 

C  2u) ) 

of  72.  Another  Interesting  feature  of  the  result  is  Chat  d  is  proportional 


i 


to  the  difference  between  M,  ,  and  M  .  Thus  a  possible  further  enhancement  of 

bb  aa 

the  second  harmonic  susceptibility  is  by  combining  nonsymmetric  compositional 

9 

grading  of  the  quantum  well  with  Che  applied  electric  field.  For  example,  the 
compositional  grading  can  be  such  chat  decreases  while  increases  when 
the  electric  field  is  applied.  In  this  paper,  we  considered  the  intersubband 
effects  't  infrared  frequencies.  The  same  idea  of  electric  control  can  also  be 
extended  to  optical  frequencies  by  considering  transitions  between  conduction 
and  valence  bands. 
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Fig •  1  PoCenClal~energy  profile  for  an  infinite  quantum  well  with  width  L 
subject  to  an  external  electric  field  F  in  tl.e  presence  of  incoming 
radiation  with  angular  frequency  Ruj 


Fig ■  2  The  normalized  second  harmonic  susceptibility  (d  as  a  function  of 

fiu)  in  eV  for  L  ■  126.5  A  T  »  77  K,  n  =  3.0  x  10^^/cm^.  1/y  “ 

aa 

1/y.  ,  *  1/y  .  “  0.14  ps,  and  F  »  10  kV/cm  and  F  «  70  kV/cm. 

00  ao 


Fig.  3  The  second  harmonic  susceptibility  |d  “|  (peak  value)  in  (1/9  x 
“22 

10  MKS  units)  is  a  function  of  electric  field  strength  F  in  kV/cm 


TABLE  I 


2(13 

A  Cable  for  comparison  of  ]d  j  in  bulk  GaAs  and 
quantum  well  with  an  applied  electric  field  F. 


macert  al 

Bulk  GaAs 

72 

GaAs 
quantum 
well  with 
an  applied 

F  = 

O.OI  kV/cm 

1.14 

field  F 

5  kV/cm 

571.22 

10  kV/cra 

1140.45 

20  kV/cra 

2265.1 3 

50  kV/cm 

5398.48 

70  kV/cm 


7 167.65 


